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ABSTRACT 


A symbolic  processor,  ( pronounced  "five  star"  ),  to  assist 
in  the  generation  of  stiffness  matrices  for  finite  elements,  based  on 
a recently  developed  symbolic  processor,  is  presented.  Operations  are 
performed  upon  element  characteristics  and  material  properties  in 
symbolic  form  to  produce  a "matrix  template,"  consisting  of  the 
algebraic  expression®  generated  for  the  stiffness  coefficients  as 
functions  of  the  problem  parameters  in  literal  form.  The  template  may 
be  evaluated  for  a given  element  by  binding  these  symbolic  forms  to 
the  numerical  values  associated  uith  a specific  element.  The 
evaluation  process  is  further  facilitated  by  permitting  specification 
of  a variety  of  output  formats  for  the  resulting  matrix  template. 
Required  input  is  minimized  by  automatically  synthesizing  the 
constituent  matrices  of  the  formulation  from  user-supplied 
specifications  of  shape  functions,  material  properties  and 
stress-strain  relationships,  all  in  symbolic  notation. 

The  processor,  written  in  RACSYI1A,  is  highly  interactive  providing 
prompts  for  user  input,  enumeration  of  available  program  options,  and 
extensive  on-line  assistance.  The  user  may  input  a "?"  in  place  of  a 
prompted  input  to  request  instructional  text.  The  file  handling 
capabilities  of  flACSYtIA  are  utilized  to  retain  a complete  record  of 
each  program  run.  These  records  facilitate  the  handling  of 
diagnostics,  assist  in  further  processing  and  permit  the  generation  of 
statistics  valuable  for  system  development.  Error  checking  is 
accomplished  through  semantic  checks  built  into  the  program  functions 
and  syntactic  checks  performed  uithin  the  MACSYflA  operating 
environment. 


A partial  list  of  user  input  includes: 

1)  Method  Selection  - Isoparametric  or  generalized  coordinate 
formulations.  2)  Element  Parameters  - Number  of  nodes,  number  of 
degrees  of  freedom  per  node  and  related  terms.  3)  Material  Properties 
- This  matrix  misy  be  selected  from  a library  of  standard  forms  (e.g. 
plane  stress,  plane  strain)  or  supplied  by  the  user.  4)  Strain 
Specification  - Components  are  entered  in  a user-oriented  calculus 
notation  (e.g.  du/dx  is  input  as  0(u,x)  ).  5)  Shape  Functions  - Shape 
functions  may  contain  trigonometric  functions  and  a large  class  of 
intrinsic  functions  as  uell  as  polynomial  terms.  6)  Output  Control 
Specification  - A description  of  the  output  format  of  the  generated 
matrix  template. 

Possible  output  forms  include  a tabular  display  of  the  matrix 
coefficients  In  symbolic  form  and  the  coefflcents  In  FORTRAN  card 
image  format. 


Background  material  includes;  The  objective  of  this  study;  The 
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(derivation  of  the  stiffness  matrices;  A summary  of  previous  research; 
A brief  (description  of  MACSYMA. 

□(itails  of  the  implementation  of  ****»•»  cover:  The  design 
objoctivQBj  Octal  lo  of  the  algorithms  used  and  how  they  were 
Implemented;  A description  of  *iViV**  and  its  limitations. 

Sample  runs  include  the  formulation,  using  both  tlie  isoparametric 
and  generalized  coordinate  methods,  of  the  stiffness  matrices  for  a; 
Bar  element  with  constant  cross-sectional  area;  Bar  element  with 
linearly  varying  cross-sectional  area;  Constant  Strain  Triangle  with 
uniform  thickness;  Four  Node  Quadrilateral. 

Conclusions  arc  drawn  and  recommendations  for  future  work  are 
made.  Appendix  I contains  notes  on  operating  and  accessing  the 
processor. 
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Chapter  1 

INTRODUCTION 


1.1  Objective 

The  finite  element  method  involves  two  processes:  the  generation 
of  elements,  where  the  computation  effort  is  linear  in  the  number  of 
elements,  and  the  solution  of  the  discrete  system  equations,  where  the 
effort  increases  as  some  power  of  the  number  of  degrees  of  freedom, 
and  thus  of  the  number  of  elements.  Improvements  over  the  past  decade 
in  decomposition  and  solution  techniques  have  reached  the  point  where 
in  many  problems  the  element  generation  effort  exceeds  that  for  system 
solution.  Thus,  from  a practical  standpoint  there  is  a great 
incentive  to  attempt  to  drastically  reduce  the  computational  effort  in 
element  generation,  the  bulk  of  which  is  taken  up  in  the  numerical 
quadrature  Involved.  With  the  large  variety  of  new  elements  being 
developed  or  investigated,  there  is  a similar  incentive  to  reduce  the 
amount  of  manual  algebraic  manipulations  required  before  the  stiffness 
matrices  for  a new  element  can  be  cast  in  a form  suitable  for 
processing. 

This  research  project  was  undertaken  to  develop  a mechanism  that 
would  facilitate  the  generation  of  element  sti ffness  matrices  by 
eliminating  the  numerical  quadrature  process  and  by  minimizing  the 
amount  of  manual  algebraic  manipulation  required. 

A symbolic  processor,  VoWrfnV,  to  assist  in  the  generation  of 
stiffness  matrices  for  finite  elements,  based  on  a recently  developed 
symbolic  processor,  is  presented.  Operations  are  performed  upon 
element  characteristics  and  material  properties  in  symbolic  form  to 
produce  a "matrix  template,"  consisting  of  the  algebraic  expressions 
for  the  stiffness  coefficients  as  functions  of  the  problem  parameters 
in  literal  form.  The  template  may  be  evaluated  for  a given  element  by 
binding  these  symbolic  forms  to  the  numerical  values  associated  with  a 
specific  element.  The  evaluation  process  is  further  facilitated  by 
permitting  specification  of  a variety  of  output  formats  for  the 
resulting  matrix  template.  Required  input  is  minimized  by 
automatically  synthesizing  the  constituent  matrices  of  the  formulation 
from  user-supplied  specifications  of  shape  functions,  material 
properties  and  stress-strain  relationships,  all  in  symbolic  notation. 

1.2  Derivation  of  the  Stiffness  flatrix 

This  section  introduces  the  nomenclature,  details  the  required 
operations  and  contains  the  derivation  of  stiffness  matrices  for  both 
the  isoparametric  and  generalized  coordinate  formulations. 

The  stiffness  matrix  for  any  finite  element  is  given,  in  general. 
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li'ij  {'.ife  rii.’tiCii  ■:in(i  Ab(;M3l.  Zi enk  i eui cz  HHi  ) : 

. / CK:]  = j fBI  IC]  (Bl  tlV  (1) 

The  termr,  i ii  equation  (1)  tire  cii  ocuoapcl  beloii. 

1.2.1  I so|;tnr orne  trie  F or  mu  I a I i on 

For  isoparametric  elements  equation  (1)  is  rei-irittcn  in  terms  of 
natural  coordinates.  In  three  dimensions  this  produces 


IK] 


III 


(Bl  [Cl  IBI  detdJl)  dr  ds  dt 


(2) 


in  i-ihich  IK]  is  the  stiffness  matrix  and  r,  s and  t are  natural 
coordinates  (Desai  and  AljellSI).  The  derivation  of  tlie  constituent 
motrices  of  cciuotion  (2)  is  given  helou. 

In  the  isoparametric  formulation.  tlie  functional  relationship 
describing  the  element  geometrg  and  the  element  displacement  are  the 
sohm;: 

Ix)  » I N(r,s,t)  I IXnl  (3) 

ful  » { N(r,s,  t)  I Iql  (A) 

Here  Ix)  represents  the  cartesian  coordinates  of  the  element;  IXn) 
is  the  vector  of  nodal  coordinates  in  the  global  coordinate  sgsteni; 
lul  stands  for  the  values  of  the  displacements  interior  to  the 
clement;  fql  represents  the  nodal  displacements  and  IW(r,s,t)l  are  tlva 
interpolation  functions  in  terms  of  the  non-dimensional  natural 

coordinates.  (Although  the  present  development  assumes  a global 

cartesian  system,  a similar  derivation  can  be  uritten  for  other 
coordinate  systems).  The  dimensionality  of  the  above  relationship  is 
clarified  in  the  follouing  expansion  of  equation  (4).  Letting  m,  = the 
number  of  nodes, 
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(5) 


The  next  ste|)  is  to  determine  the  strains  uhich  arc  derivatives  of 
the  displacements: 


le)  - IB)  Iq) 


(6) 


in  uhich  ie)  represents  the  strain  components  defined  in  the  global 
coordinate  system;  and  (Bl  represents  the  derivatives  of  the 
interpolation  functions  in  equation  (A)  uith  respt’ct  to  th(?  global 
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coordinoles. 

Since  these  i n t er|30 1 a t i on  functions  (and  hence,  tlio  di  tp  I aceinont  o) 
are  functions  of  tlio  natural  coordinates,  differentiation  must  he 
performed  by  the  chain  rule.  Employing  the  derivative  of  eqiiotion 


with  respect  to  each 
chain  rule  produces': 


natural  coordinate  in  the  expression  for 


(3) 

the 


dN 

3r 

dW 

aN 

at 


in  which  the  Jacobian,  [J] 

'aw 
ar 
aivi 
as 
aN 
at 


[j] 

3x3 


3xm 

is  defined  by 


an 

ax 

aN 

aij 


(7) 


3xm 


[J] 


I (Xnl  (Ynl  IZnI  I 


111x3 


(S) 


3xffl 


As  the  quantity  on  the  right  side,  of  equation  (7)  is  required  in 


the  formulation  of 


(Bt 


he 


i over  ted. 


the  Jacobian  must 
Preniu  I t i p I y i ng  both  sides  of  equation  (7)  by  tlie  inverse  of  the 
Jacobian  produces  an  expression  for  the  derivatives  of  the  sltape 
functions  with  respect  to  the  global  coordinates.  These  derivatives 
are  assembled  into  the  (Bi  matrix  and  ordered  according  to  the 
specifications  given  by  the  strain  component  vector,  lei. 


The  [CJ  matrix  contains  the  stress-strain  relationships: 

(si  - (C]  (el 


(9) 


Finally,  when  the  volume  integral  is  converted  from  global  to 
natural  coordinates,  the  differential  volume  becomes 


dx  dy  dz 


det([J))  dr  ds  dt 


(10) 
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in  which  IJ]  is  the  Jacobian  os  dofinecJ  in  etiuation  (8). 


1 . 2. 2 Cenera I i zed  Coordinate  Formulation 


The  generalized  coordinate  formulation 
expressing  lul  , the  displacements  internal 
function  of  a set  of  yt.  to  be  determined  generalized 
represented  by  the  column  vector,  lai : 


begins  with  a relation 
to  the  element  as  a 
coordinates 


(u! 


{S(x)l  lai 


(11) 


The  shape  functions,  (Si,  are  polynomials  in 
coordinates  (denoted  by  x)  and  are  chosen  to  conform  t 
requirements.  In  general,  the  order  of  the  polynomials 
tlie  number  of  generalized  coordinates  is  equal  to  the  to 
degrees  of  freedom  of  the  element.  Utilizing  this  pr 
displacements  interior  to  the  element  may  be  expressed  in 
nodal  displacements  as  follows.  Substituing  the  nodal 
into  the  shape  functions  produces  the 
matrix,  (A),  which  relates  the  nodal 
generalized  coordinates: 


di sp I acement 
di sp I acement s. 


lc|i 


Solving  for  (ai 


(A)  (ai 


-1 


the  global 
o convergence 
is  such  that 
tal  number  of 
inciple,  the, 
terms  of  the 
coord i nates 
ransf ormat i on 
(qi  , to  the 


(12) 


(ai  = (A)  (qi 


(13) 


and  substituting  into  equation  (11)  produces  the  desired 
relationship  between  nodal  and  element  displacements 


(ui 


(S(x)i  (A)  (qi 


(14) 


Thio  form  Is  fundamentally  the  same  as  equation  (4)  for  the 
isoparametric  case. 


For  a specified  set  of  strain  components,  the  (Bi  matrix  for  the 
generalized  coordinate  approach,  (Ba)  , may  be  determined  by 
appropriate  differentiation  of  equation  (11). 

(el  - (Bai  (ai  (15) 

Since  the  str.iins  are  derivatives  of  the  displacements  with 
respect  to  the  global  coordinates  and  the  shape  functions  are 

polynomials  in  these  coordinates,  no  coordinate  system  transformation, 
as  in  equation  (7)  for  the  isoparametric  case,  is  required. 

The  material  properties  matrix,  (Cl,  is  the  same  as  in  the 

isoparametric  method. 
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SubBtituting  the  constituent  matrices  into  equation  (1)  and 
integrating  with  respect  to  the  differential  volume,  dV,  in  terms  of 
the  global  coordinates,  produces  a stiffness  matrix,  [Kal , associated 
with  tfie  generalized  coordinates. 

IKaJ  - J IBa)^  [C]  IBal  dV  (IG) 

Equation  (13)  gives  the  transformation  necessary  to  express  the 
stiffness  matrix  with  respect  to  the  nodal  dicplacemcnt  quantities:. 

-1  t -1 

IK]  « [A]  (Ka)  (A)  (17) 

Equations  (2)  and  (17)  are  the  expressions  for  the  stiffness 
matrix  employed  in  the  computations  perfomed  by  the  processor. 


1.2.3  Formulation  of  Element  Loads,  Mass  and  Stresses 

In  addition  to  the  stiffness  matrix,  several  other  element 
matrices  are  relevant  to  finite  element  analysis.  These  matrices 
include  the  element  mass  matrix,  the  stress  - displacement 
relationships  and  the  load  vectors  representing  (a)  body  forces,  (b) 
surface  tractions  and  (c)  initial  strains. 

Using  the  nomenclature  introduced  in  section  1.2.1  and  lettit'g  (L) 
represent  the  shape  functions  for  either  the  isoparametric  or 
generalized  coordinate  methods,  the  formulation  of  the  above  matrices 
fol lows: 


The  element  body  force  vector,  (QDl  is  given  by 


(QBI 


(fbi  dV 


(18) 


in  which  (fbi  represents  the  specified  body  forces.  If  (L)  contains 
the  generalized  coordinate  shape  functions,  (QB)  must  be 
pre-mu  1 1 i p I i ed  by  the  inverse  - tranpose  of  the  [A]  matrix  of  equation 
(12). 


The  element  surface  traction  vector,  IQTI , is  given  by 


IQTl 


Ift)  dV 


(19) 


in  which  (ftl  represents  the  specified  surface  tractions  and 
integration  is  performed  over  the  surface  of  the  element.  The 
transformation  involving  the  (A)  matrix  would  be  applied  as  in  the 
formulation  of  IQBI  . 


The  element  initial  force  vector,  IQI) , is  given  by 


1.3  Motivation  for  Symbolic  Processing 

Numerical  integration  techniques  require  the  evalvation  of  the 
integrand  of  equation  (1)  at  specified  points,  producing  a square 
matrix  of  numerical  values  of  order  equal  to  the  number  of  degrees  of 
freedom.  As  constituent  terms  in  this  equation  are  functions  of  the 
nodal  coordinates  (or  are  geometric  values  which  must  be  evaluated  at 
these  points)  the  evaluation  must  be  performed  on  each  element 
individually.  If  a large  number  of  elements  is  to  be  employed  in  the 
finite  element  model,  a significant  computational  overhead  may  be 
1 ncurred. 

Similarly,  numeric  evaluation  of  the  inverse  and  determinant  of 
the  Jacobian  (equation  (7))  and  the  inverse  of  the  displacement 
transformation  matrix  (equation  (13))  also  requires  nodal  coordinates 
to  be  bound  to  numeric  values  and  thus  must  be  performed  upon  elements 
individual ly. 

Using  symbolic  manipulation  techniques,  nodal  coordinates  may  be 
retained  and  operated  upon  in  literal  form  throughout  the  computation. 
The  coordinates  are  not  bound  to  numerical  values  and  thus  may 
represent  any  set  of  actual  (numerical)  coordinates.  The  stiffness 
matrix  produced  is  expressed  as  a template  in  terms  of  these  unbound 
values.  It  may  be  readily  evaluated  for  a given  set  of  actual  nodal 
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coordinates  during  execution.  Ihe  computations  detailed  in  the 
previous  section  need  only  performed  once  to  generate  this  general 
template. 

/ 


1.4  Summary  of  Previous  Uork 


Several  researchers  have 
symbolic  computing  languages 
matrices  directly  in  literal 
the  specific  numerical  values 


proposed  or  investigated  the  use  of 
for  generating  finite  element  stiffness 
form,  to  be  subsequently  evaluated  for 
of  a given  element. 


Luftlll  proposes  individual  special  purpose  routines  with 
limitations  as  to  problem  type,  element  shape,  and  displacement 
function  specification.  Two  example  programs  are  provided.  The 
discussion  for  the  Tirst,  a processor  for  rectangular  elements  using 
generalized  coordinates,  indicates  that  exact  integration  is  being 
performed  because  the  element  boundaries  correspond  to  constant  values 
of  X and  Y.  A program  employing  an  isoparametric  formulation 
comprises  the  second  example.  Since  it  utilizes  Gauss  Quadrature  and 
the  numerical  evaluation  of  the  Jacobian,  it  must  be  executed  for  each 
element  individually.  Contributions  detailed  include:  (a)  the 
introduction  of  'intrinsic  matrices'  (the  integral  of  the  product  of  a 
form  of  the  interpolation  functions  with  its  transpose)  for  minimizing 
and  organizing  intermediate  calculations  and;  (b)  the  detailed 
specification  of  a polynomial  manipulator  for  performing  the  requisite 
operat i ons. 

The  processor  descibed  in  Gunderson (2)  employs  a modified 
generalized  coordinate  approach  and  requires  the  user  to  specify  the 
displacement  transformation  matrix  (i.e.  the  inverse  of  the  matrix 
relating  the  generalized  coordinates  to  the  nodal  displacements).  The 
program  makes  use  of  a sophisticated  scheme  for  data  organization  by 
representing  the  polynomials  as  multi-dimensional  integer  arrays  and 
by  defining  the  matrix  operations  accordingly. 


Taigl3) 
coef  f i c i ents 
the  cases  of 


has  presented  programs  for  evaluating  the  stiffness 
for  quadrilateral  plate  elements  with  in-plane  forces  for 
rectangular  and  trapezoidal  panels. 


Anderson  and  Noor  141  and  Anderson  and  Bowen (5)  demonstrate  the  use 
of  HACSYflA's  symbolic  integration  facilities  in  the  development  of 
shallow-shell  finite  elements.  In  both  papers  the  required  integrals 
are  divided  into  classes  using  group  theoretic  techniques.  In 
reference  (5),  the  symbolic  expressions  are  incorporated  into  a 
FORTRAN  computer  program.  A study  is  made  of  computation  time  and 
memory  requirements. 

Wong  IG)  suggests  the  use  of  symbolic  computation  in  the 
formulation  of  finite  element  stiffness  matrices  but  admits  that 
CONFORM,  the  system  that  he  co-authored.  Is  insufficiently  flexible 
for  the  required  handling  of  data  structures. 
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The  survey  in  Jensen  113)  of  symbolic  computing  languages  and  their 
applications  in  mechanics  suggests  that  only  flACSYUA  possess 
sufficient  flexibility  and  pouer  to  carry  out  the  objectives. 

In  summary  the  uork  to  date  has  been  directed  touarcle  the 
development  of  eitheri  (a)  programG  that  take  advantage  of  properties 
associated  with  a very  specific  element  type  ( 13),  14),  15)  )j  or  (b) 
processors  which  synthesize  matrices  for  any  element  but  cannot 
automate  the  entire  process  ( 111,  12)  ). 


1.5  Description  of  HACSYHA 


Thi  s 
MACSYMA, 
m.iter  i a I 


section  contains  a brief  description  of  the  capabilities  of 
the  base  language  of  jVjVjVjV*,  and  wi  I I serve  as  background 
for  subsequent  dicussion  of  the  implementation  of  the  system. 


MACSYtlA  (pronounced  "maxima").  Project  flAC’ s SYmbolic  Manipulation 
system,  is  a large  computer  program  written  in  LISP  devoted  to  the 
manipulation  of  algebraic  expressions.  I1ACSYMA  runs  under  the  ITS 
timesharing  system  (originally  developed  at  the  M. I.T.  Artificial 
Intelligence  Laboratory),  on  the  Mathlab  PDP-10  computer  at  M. I.T. 
With  a syntax  resembling  ALGOL  B0,  MACSYMA  has  capabilities  for 
manipulating  algebraic  expressions  involving  constants,  variables  and 
functions  til. 


Provisions  and  attributes  of  MACSYMA  of  interest  to  the  computer 
language  designer  include: 


(a)  The  base  language  is  a specially  designed,  enriched  version  of 
LISP  called  MACLISP; 

(ta)  The  system  contains  a large  body  of  intrinsic  functions; 

(c)  ALGOL-like  control  structure,  compound  statements  and  block 
structure  permit  the  incorporation  of  the  intrinsics  into  user-defined 
functions.  As  in  LISP,  all  functions,  user-defined  or  intrinsic, 
return  a value; 

(d)  A set  of  commands  may  be  pre-stored  on  a disk  file  and  executed  by 
means  of  the  ’BATCH’ * command; 

(e)  An  editor,  modeled  after  TECO,  can  be  Invoked  to  edit  input  or  to 
correct  syntax  errors; 

(f)  Data  types  include  atomic  variables,  lists,  arrays,  matrices,  and 
strings.  Numeric  constants  m,iy  be  integers,  ratioral  numbec,  floating 
point  numbers  or  "bigfloats"  (floating  point  va'.ues  of  essentially 
arbitrary  precision); 

(g)  Debugging  aids  and  trace  functions  are  provided; 

(h)  The  user  may  deciare  and  manipulate  properties  of  atoms; 
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(i)  Pattern  matching  facilities  exist.  These  include  type  testing  and 
general  pattern  matching  functions  which  permit  the  user  to  test 
expressions  for  combinations  of  syntactic  and  semantic  patterns  and  to 
automatically  have  variables  set  to  parts  of  the  expressions  which  fit 
the  patterns  IGl ; 

(j)  It  is  possible  to  store  expressions,  values  and  functions  on  disk 
f i I es. 

Mathematical  and  computational  functions  include: 

(a)  Evaluation  and  simplification; 

(b)  Differentiation  and  integration; 

(c)  Part  selction  and  substitution; 

(d)  Solving  for  roots  of  an  equation; 

(e)  Matrix  functions  including  transpost i on,  multiplication,  inversion 
and  evaluation  of  determinants. 

(f)  Manipulation  of  rational  expressions  (expressions  which  are  the 
quotient  of  two  polynomials); 

(g)  Taylor  series  and  power  series; 

(h)  Graphing; 

(i)  Poisson  series; 

(j)  Tensor  manipulation; 

(k)  Laplace  transforms; 

(l)  Finding  the  limit  of  an  expression  as  a constituent  variable 
approaches  a value  from  a given  direction. 

Command  lines  to  MACSYMA  are  strings  of  characters  representing 
mathematical  expressions  involving  equations,  arrays,  functions,  and 
programs.  Extra  spaces,  tabs,  and  all  carriage  returns  are  ignored 
(except  when  these  occur  in  quoted  strings)  (21. 


Comm,ind  lines  are  terminated  by 


or  "S"  (dollar  sign).  A "; 


causes  the  command  line  to  be  evaluated  and  the  result  di  spayed  . The 
terminator  "8"  causes  the  command  line  to  be  evaluated  but  the  result 
is  not  displayed  (3).  As  *>Wn‘nV  supervises  the  display  of  output, 
both  terminators  have  the  same  effect. 

The  command  (input)  lines  are  indexed  by  labels  of  the  form  "(Ci)" 
where  i is  incremented  with  each  new  command  typed  by  the  user. 
Similarly,  the  results  of  computations  ore  also  indexed  by  a label  of 
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the  form  "(Ul)'*;  thuo,  uoual  lij  the  i th  input-output  pair  Mill 
be(Ci)-(Di)  15].  Intermediate  results  (if  any)  are  tagged  uith  line 
latsels  of  the  form  "(Ei)".  Line  labels  may  also  be  used  to  reference 
associated  expressions. 

Bogen  171  contains  a description  of  the  major  features  and 
functions  of  nACSYMA.  flathlab  IS)  provides  details  on  logging  into  the 
system  and  contains  a script. 

An  Important  characteristic  is  that  flACSYMA,  In  its  present, 
experimental  version,  provides  a significantly  limited  amount  of 
user  storage  space  for  symbolic  expressions. 


l.G  Organization  of  This  Report 

The  remainder  of  this  report  is  organized  as  follows; 

Chapter  2 describes  the  system  implementation,  including  the 
Design  Consideration  and  Constraints;  Details  of  the  implementation  of 
the  computation  of  the  stiffness  matrix;  and  constituent  matrices; 
Specification  of  algorithms;  Logic  Hierarchy  and  Data  Flow;  Unique 
Input  and  Output  Facilities;  Error  Recovery;  and  System  Limitations. 

Chapter  3 contains  a number  of  illustrative  examples  of  runs  made 
on  the  system  and  the  verification  of  results  for  several  test  cases. 

Chapter  4 consists  of  a summary  and  conclusions  giving  a final 
assessment  of  the  work  performed  and  suggestions  for  future  research. 

Appendix  1 is  a brief  user’s  manual  to  flACSYflA  and  Vo'n'o’nV, 
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Chapter  2 

SYSTEM  IMPLEMENTATION 


2.1  Implementation  Considerations 

The  considerations  in  the  design  of  ****»•#  uere  Imposed  by  the 
operating  environment  of  MACSYMA  and  a set  of  initial  design 
decisions.  The  major  design  decisions  were  as  follows: 

(a)  The  system  had  to  be  easily  usable  by  engineers  possessing  a 
minimum  knowledge  of  its  operation. 

(b)  Ease  of  maintenance  and  modification  were  be  emphasized.  The 
incorporation  of  improved  algorithms  and  additional  methods  of 
analysis  would  thus  be  facilitated. 

(c)  The  modularization  of  the  system  was  intended  to  resemble  the 
steps  in  the  problem  formulation. 

(d)  Consistent  with  MACSYMA’ s interactive  mode  of  operation,  the 
system  would  operate  as  an  interactive  process.  Techniques  which 
provide  selective  ouput  were  necessary  to  override  MACSYMA’s  procedure 
of  echoing  commands  and  output. 

(e)  Since  computation  time  Increases  with  the  number  and  complexity  of 
the  expressions,  in  general,  there  is  no  "compute  versus  store" 
tradeoff  in  MACSYMA.  Processes  may  also  become  storage  bound  with 
respect  to  invoked  program  units,  both  user-defined  and  intrinsic. 
Whereas  user  program  modules  can  be  deleted  in  order  to  free  storage, 
MACSYMA  program  segments,  once  loaded,  remain  for  the  duration  of  the 
job.  A constraint  was  also  placed  upon  computation  time  because  the 
system  is  to  be  accessed  interactively. 

(f)  All  input  had  to  be  format-free. 

(g)  Semantic  and  syntactic  checks  on  input  had  to  performed. 
Integration  of  MACSYMA’s  syxtax  checker  and  editor  was  deemed 
appropr i ate. 

(h)  It  was  decided  that  a copy  should  be  retained  of  each  session 
during  which  the  system  is  used.  Such  records  would  aid  in 
diagnostics  and  development. 

(i)  As  the  system  is  to  run  interactively,  a help  facility  had  to  be 
provided  to  supplement  error  messages  generated  by  semantic  checks  as 
well  as  to  assist  in  entering  input. 

As  an  alternative  to  the  general  method  of  computation,  which 
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produces  the  complete  matrix  template  by  performing  all  integrations 
and  matrix  operations  symbolically,  a hybrid  symbolic-numeric  scheme 
was  also  considered.  In  this  approach,  symbolic  operations  are  used 
only  to  produce  the  template  for  the  triple  product  Integrand  of 
equation  (1)  (section  1.2),  in  terms  of  unbound  values  of  the 
coordinates  at  the  numeric  quadrature  points.  At  execution  time  of 
the  analysis  program  invoking  it,  the  template  is  numerically 
evaluated  at  each  quadrc7ture  point  to  produce  the  integral.  The 
hybrid  technique  would  use  the  same  computation  sequence  as  either  the 
isoparametric  or  generalized  coordinate  formulations,  up  to  the 
symbolic  integration  processing. 


2.2  Computation  Implementation 

Algorithms  employed  to  perform  the  comivutat  i ons  described  in 
section  1.2  are  presented. 

2.2.1  Reformulation  of  the  Stiffness  Matrix 

To  minimize  computational  effort,  the  formulation  of  the  stiffness 
matrices  given  in  equation  (2)  for  the  isoparametric  method  and 
equation  (17)  for  the  generalized  coordinate  approach  may  be 
reorganized  as  follows. 

For  the  isoparametric  formulation,  expressing  the  Inverse  of  the 
Jacobian,  (J) , of  equation  (8)  as  the  adjoint  divided  by  the 
determinant,  the  latter  may  be  factored  from  (Bl  and  the  transpose  of 
(Bl  . Combining  these  terms  with  the  determinant  from  equation  10  and 
using  (BJl  to  denote  the  reduced-  (Bl  matrix,  yields  the  following 
expression  for  the  stiffness  matrix 


[K]  - / IBJ)  (Cl  IBJ)  dr  ds  dt  (23) 


^ det((J]) 

This  formulation  produces  a greatly  simplified  (Bl  matrix  and 
reduces  the  number  of  required  divisions  in  the  finished  matrix 
template, 

A similar  procedure  is  employed  to  factor  the  [A]  matrix  in  the 
generalized  coordinate  method.  Letting  (AOJ  represent  the  determinant 
of  the  submatrices  comprising  the  (A)  matrix  and  (AJl  stand  for  the 
form  of  the  lA)  matrix  produced  by  replacing  these  submatrices  with 
their  adjoint,  equation  (17)  may  be  rewritten  as 


1 t 

IKl  [AJ]  (Kal  [AJ]  (24) 


de  t ( lAD) ) iWr2 

Once  again  the  computations  are  greatly  simplified. 
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2.2.2  Reduction  of  Dependent  Natural  Coordinates 

The  numljer  of  natural  coordinates  specified  by  the  user  must  be 
equal  to  or,  at  moot,  one  greater  than  the  number  of  global 
coordinates  specified.  In  the  latter  case,  the  natural  coordinates 
are  not  independent. 

To  facilitate  computation  the  last  coordinate  is  treated  as  the 
dependent  one  and  is  replaced  by  an  expression  equal  to  unity  minus 
the  sum  of  the  remaining  coordinates. 


2.2.3  Generation  of  the  Jacobian 

The  Jacobian,  as  defined  in  equation  (8),  is  formed  using  a 
refined  version  of  the  following  algorithm. 

Let 

NDOF  - the  number  of  degrees  of  freedom  per  node, 

NN  » the  number  of  nodes, 

G[i]  - the  name  of  the  i th  global  coordinate, 

L[i3  - the  name  of  the  1th  natural  coordinate, 

S ■ the  vector  of  shape  functions, 

□IFF,  GC  and  JVAL  are  local  values. 

Then 

1 for  i«l  to  NDOF  do 

2 DlFF-(the  derivative  of  S with  respect  to  L[i3) 

3 for  j»l  to  NDOF  do 

4 GC-GtJ] 

5 JVAL»0 

G for  k-1  to  NN  do 

7 JVAL-JVAL  + DlFF(k],-<GC[k] 

8 JACOBIANIi,  j]-JVAL 

Note  that  GC  is  assigned  the  name  of  a global  variable  (say  X),  in 
line  4 and,  in  line  7,  assumes  symbolic  values  of  associated  nodal 
coordinates  ( X(k)  for  k»l  to  NN  ) simply  by  appending  the  subscript, 
"[k]". 


2.2.4  Formulation  of  the  Inverse  of  the  lA]  Matrix 

The  inverse  of  the  tA]  matrix  of  equations  (12),  (13),  (14)  and 

(17)  is  determined  using  the  following; 

Expanding  equation  (12),  for  the  case  of  3D  cartesian  coordinates 
produces 
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in  uiiich  'in'  represents  the  number  of  nodes  and  (Lin),  (Vnl  , (l-lnl  are 
the  nod,Tl  values  of  the  displacements.  The  vectors  (Sll,  (S2)  and 
(531  are  the  shape  functions  for  each  of  the  degrees  of  freedom  ( U, 
V,  U ) , evaluated  at  the  nodal  values.  That  is,  the  i th  row  of  (Sjl 
is  the  value  of  the  shape  function  for  the  jth  degree  of  freedom 
evaluated  at  the  ith  node. 

For  tlie  case  in  which  each  degree  of  freedom  is  governed  hg  tlie 
same  shape  function,  all  ISjl  sul^matriccs  are  ec|ual.  The  inverse  of 
tl'ie  [AJ  matrix,  which  is  of  order  ( 3m  >:  3'ii  I,  can  l.■>p  obtained  by 
inverting  one  of  the  ISjl  matrices  of  order  I m x m I and  assembling  N 
copies  (tihere  N is  tlie  number  of  degrees  of  freedom  per  node)  along 
the  main  diagonal  of  a previously  zeroed  matrix.  In  addition,  the 
determinant  of  (Sjl  m,iy  be  factored  from  the  inverse. 

For  unique  values  of  the  ISjl.  each  submatrix  would  have  to  Ije 
inverted.  In  either  case,  the  computational  savings  is  appreciable: 
Rather  than  inverting  a matrix  of  order  ( Nm  x Nm  ),  only  K inverses 
of  order  ( m x m ) need  be  taken  where  K is  equal  to  unity  if  all 
shape  functions  are  the  same  or  N if  all  are  unique. 


2.2.5  Processing  of  Strain  Components 

MACSYIIA’s  General  Pattern  Hatching  Functions  are  used  to  convert 
user  specifications  of  strain  components  into  a database  from  which 
the  (Bl  matrix  for  the  isoparametric  formulation  and  the  fBal  matrix 
for  the  generalized  coordinate  method  are  created. 

Th«j  procedure  consists  of  associating  predicates  with  pattern 
variables  and  defining  functions  of  forms,  containing  these  variables, 
against  which  input  m.oy  be  tested. 

This  process  permits  components  to  l5e  specified  in  a calculus 
notation.  For  example,  3u/c)x  is  input  as  Dfu.x)  and  cl*u/dxc)y  may  be 
entered  as  D(u,x,y).  Currently  permissilile  forms  for  components  are 

a * 0(u,x) 

a * niu.x)  + b v<  Dlu,x) 
u/a 

ubere  a,  b represent  scalars, 

u represents  a displacement  variolate, 

X represents  a global  coordinate  variable. 


The  present  implementation  will  parse  terms  specifying  second 
derivatives  but  will  output  a warning  that  the  incorporation  of  such 
derivatives  in  the  IBI  or  IBal  matrices  is  beyond  present 
capabi I i t i es. 

The  database  is  in  matrix  form  where  each  row  contains  the 
information  associated  with  a term  in  the  specified  strain  component. 
A singie  component  may  have  more  than  one  algebraic  term.  Letting 
DB[i,j]  denote  a particular  element  in  the  database,  the  meaning  of 
the  entries  in  the  i th  row  is: 

□Bli.l]  » the  number  of  the  component,  which  is  the  same  as  the  row 
that  the  term  will  occupy  in  the  IBl  or  iBal  matrix, 

□B[i,2]  a the  scalar  multiple, 

DB[I,3]  - index  of  the  displacement  variable  with  respect  to  the  list 
of  displacement  variables, 

DB[i,4]  - index  of  the  first  global  coordinate  variable  in  the 

derivative  with  respect  to  the  list  of  coordinate  variables,  taqged  as 
the  number  of  degrees  of  freedom  plus  unity  for  components  of  the 
form,  "u/a". 

□B[i,5]  - Index  as  in  column  4 but  for  the  second  coordinate  variable 
( if  any  J in  the  derivative,  set  to  zero  since  second  derivatives  are 
not  yet  implemented. 

As  an  example,  the  database  for  the  strain  components, 

{ du/3r,  u/r,  du/dz,  du/dz  + du/dr  ) 


given  the  list  of  displacement  and  coordinate  variables  as  t r,  z ] 
and  [ u,  w ] respectively  Is 

11110 

2 1/r  1 3 0 

3 12  2 0 

4 112  0 

4 12  10 

The  strain  components  would  be  entered  as 

{ D(U,R),  U/R,  D(U,Z),  D(U,Z)+0(U,R)  I 

The  user  is  given  the  option  of  either  specifying  the  strain 
components  using  the  above  notation  or  by  making  a selection  from  a 
library  of  pre-stored  databases. 

The  current  library  options  are: 


(1)  User-Supplied  Values 
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(2)  One  Oi mnnsi ona I Elasticity 

(3)  Plane  Stress 

(4)  Plane  Strain 

(5)  Ax  i syinme  tr  i c 

(G)  Linear  Isotropic  Elasticity  - 30 

Options  (3)  and  (4)  reference  identical  values  and  are  listed  as 
two  separate  options  only  for  conformabi  I i ty  with  the  Material 
Properties  library  options  (see  section  2.7).  Option  (1)  is  provided 
in  the  event  that  the  user  has  designated  library  specification  but 

the  provided  options  are  not  appropriate. 

« 

Examples  of  both  user  and  library  specifications  are  presented  in 
Chapter  3. 


2.2.6  Generation  of  the  IBl  Matrices  ] 

The  synthesis  of  the  IBl  or  the  (Bal  matrix  requires  the  I 

generation  of  the  appropriate  derivatives  of  the  shape  functions  and  • j 

the  assembly  of  these  derivatives  according  to  the  specification  of  j 

the  database  presented  in  section  2.2.5. 

. ] 

For  the  generalized  coordinate  formulation,  the  algorithm  for  the  1 

synthesis  of  the  (Bal  matrix  may  be  abstracted  asj  j 

Let  I 

NBC  » the  total  number  of  strain  components,  i 

NST  ■ the  total  number  of  strain  terms,  | 

NN  « the  number  of  element  nodes, 

NDOF  » the  number  of  degrees  of  freedom  per  node,  ! 

S ■»  the  vector  representing  the  shape  function,  | 

D,  FACTOR,  ROU,  START,  and  COL  are  local  variables.  j 

1 

then  j 

1 Zero  a matrix  of  order  NSC  by  (NN  NOOF)  ] 

2 for  i»l  to  N5T  do 

3 START  - (0BIi,3]  - 1)  * NN  I 

4 if  DB ( i , 4] -NOOF+1  I 

5 then  D ■ S . j 

B else  0 ■ (the  derivative  of  S with  respect  to  the  i 

coordinate  variable  specified  in  DBli,4])  j 

7 FACTOR  - DB(i.2]  | 

8 ROU  - OBIi.l] 

9 for  j«l  to  NN  do 

10  COL  - START  + J j 

11  Ba (ROU, COL]  -Ba  (ROU, COL]  + FACTORED 


The  algor  itm  for  the  generation  of  the  (Bl 
isoparametric  formulation  is  presented  in  two  parts: 
the  derivatives  and  the  assembly  process. 


17 


SYSTEM  IIIPLEHENTATION 


To  cjenerate  a matrix,  DTEMP,  the  rowo  of  which  contain  the 
derivatives  of  the  shape  functions,  S,  with  respect  to  the  global 
coordmates,  perform  the  following: 

Let 

NDOF  ■ the  number  of  degrees  of  freedom  per  node, 

NN  » the  number  of  nodes, 

S  ••  the  vector  of  shape  functions, 

D,  JVAL  are  local  variables. 


Then 


1 Generate  ADJJ,  the  adjoint  of  the  Jacobian 

2 for  j=l  to  NDOF  do 

3 D » (derivative  of  S with  respect  to 

the  Jth  natural  coordinate) 

4 for  i=l  to  NDOF  do 

5 JVAL  - ADJJti.J] 

G for  k.«l  to  NN  do 

7 DTEn(Mi,k]  =DTEI1P[i.k]  + JVAL*D(k) 

8 Append  a copy  of  S (for  the  0th  derivative)  as  the  last  ( i.e. 
NDOF+lst  row  ) to  facilitate  processing  of  terms  of  the  form 
"u/a“  . 

The  process  of  assembling  the  derivatives  to  form  the  (Bi  matrix 
is  abstracted  in  the  following  algorithm  which  closely  resembles  its 
counterpart  in  the  generalized  coordinate  approach. 

Let 

NN,  NDOF,  DB,  DTEMP,  NSC,  NST  « as  previouly  defined 
STARTl,  START2,  ROW,  FACTOR,  ROW,  COL  are  local  variables 

then 

1 Zero  a matrix,  IB)  of  order  NSC  by  (NN  * NDOF) 

2 for  i«l  to  NST  do 

3 STARTl  - (DB[i,3]  - 1)  * NN 

4 START2  - DB[i,4] 

5 FACTOR  - DB(i,2] 

6 ROW  - OB[i,l] 

7 for  j«l  to  NN  do 

8 COL  - STARTl  + j 

3 B(R0W,C0L1  - B[R0W,C0L1  + FACTORvfOTEnP  [START2,  j] 

Whereas  the  algorithm  for  the  generalized  coordinate  approach 
generates  the  derivatives  ’on-the-f ly’ , the  derivatives  referenced  in 
the  algoritm  above  are  prestored  in  the  DTEMP  array. 


To  identify  the  occurrence  of  the  "u/a"  forms,  DB[i,4]  and  hence 
START2  (line  4)  were  set  to  NDOF+1.  The  algorithm  generating  DTEMP 
appends  the  0th  derivative  of  the  shape  function  os  the  NDOF+lst  row. 
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2.2.7  Material  Properties  Specification 

the  user  has  two  options  in  the  specification  of  the  Material 
Properties  Matrin  of  equation  (9)  (see  section  1.2.1): 

For  the  "User-Supplied"  option,  the  user  may  enter  the  upper 
triangular  portion  of  the  material  properties  matrix  in  row  major 
format.  Specification  of  a constant  scalar  multiple  is  also  accepted. 
Editing  functions  which  permit  the  display  or  modification  of  either 
the  matrix  or  scalar  multiple,  are  provided; 

The  "Library"  option  permits  access  to  a set  of  pre-stored 
material  properties  matrices.  The  library  options  coincide  with  those 
for  the  strain  component  specification: 

(1)  User-Supplied  Values 

(2)  One  Dimensional  Elasticity 

(3)  Plane  Stress 

(4)  Plane  Strain 

(5)  Ax  i syinme  tr  i c 

(G)  Linear  Isotropic  Elasticity  - 30 

Option  (1)  is  provided  in  the  event  that  the  user  has  designated 
library  specification  but  the  provided  options  are  not  appropriate. 

The  library  matrices  use  the  symbols,  "E"  and  "NU"  to  represent 
Young’s  Modulus  and  Poisson's  Ratio  respectively.  Until  an  edit 
function,  which  would  permit  the  renaming  of  these  values,  is  created, 
the  user  should  avoid  using  these  names  in  the  specification  of  other 
var iables. 


2.2.8  Element  Volume  Modification 

The  nominal  elemental  volume  is  calculated  as  the  product  of  the 
differential  forms  of  the  global  coordinates  for  the  generalized 
coordinate  approach  or  of  the  natural  coordinates  for  the 
isoparametric  formulation. 

The  user  may  modify  these  forms  by  inputting  appropriate  scalar 
multiples.  For  example,  in  axi  syinmo  tr  i c problems,  the  user  might 
enter  the  factor,  2nr,  as  a multiple  of  the  products  of  differentials? 
dr*dz.  The  modification  terms,  also  referred  to  as  Auxiliary  Terms, 
may  contain  any  syntacticai  ly  correct  arithmetic  expressions. 
Examples  of  the  use  of  element  volume  modification  terms  are  presented 
in  Chapter  3. 


These  values  are  placed  on  a stack  upon  input.  The  system  uses 
the  same  stack  to  retain  auxiliary  terms  encountered  during 
computation  (for  example,  the  determinant  of  (A)  or  of  the  Jacobian). 
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2.2.3  IJti  term  i iv)  I i on  of  the  Limit'!  of  I nto'jr.i  I i on 

In  tl'Kj  c)c?nero  I i ^ecl  coordinate  formulation,  numeric  values  foi-  the 
limits  of  integration  cannot  l>e  estal)  I i sliecl  until  tlie  element  is 
located  in  the  finite  element  menh:  onlij  tlitn  caii  maxiriu.'tn  and  minimum 
coordinate  valuer,  he  determined.  To  c i rc:uinve'nt  this  (M'ohlem  the 
limits  are  taken  as  identifiers  formed  hy  concatenat i ng  the  coordincite 
varial'jlp  n.'imcr.  uith  the  strings  "I1IN"  and  "I'lAX"  (for  example  Xf1lN, 
XMAX,  YTUN  and  Yl'IAX)  . The  matrix  template  pi'odur.ed  is  tlius  a function 
of  these  values.  Numerical  values  for*  tlie  -f1IN  and  -I1AX  variables  may 
be  established  at  execution  time  by  a|:>plying  the  FORTRAN  intrinsics 
AMIN.l  and  AUAXl  to  the  approp.  iate  nodal  coordinate  vector. 

The  cle  t er  m i na  t i on  of  the  limits  of  ititegration  for'  the 
i so|:rarame  ti' i c fot'mulation  is  conditional  upon  whether-  the  natural 
coord  i ria  les  are  independent  or  dependent. 

If  the  coordinates  are  independent,  the  user-  r.|;'ec:  i f i cs  limits  from 
tire  sri't:  I -I,  (1,  1 I uith  the  stipulation  that  the.  lower  limit  is 

less  than  the  upper  limit. 

if  the  coordinates  are  dependent,  the  loner  limit  may  be  specified 
but  the  upper  limit  is  prescribed  by  tlie  following  relationship: 

Let  N(i)  denote  the  i th  natural  coordinate.  Then 


upper  limit  of  N ( i ) 


1 if  i = 1 

i-1 

1 - 5]  NIj)  otherwise 

j = l 


2.,?. IB  Integration  Processing 


The  steps  in  performing  the  integration  for  both  the  generalised 
coordinate  and  i sopor am, itr i c methods,  as  given  in  equations  (2)  and 
(IG)  respectively,  are  comparable: 


(a)  Form  the  quadratic  form  involving  the  (Cl  matrix  and  either  (Bi  or 
IBal  . 


(b)  tlultiply  the  result  by  any  auxiliary  terms  or  element  volume 
modifier  terms  which  are  functions  of  the  variables  of  integation. 

(c)  Integrate  the  terms  in  the  upper  triangular  portion  witli  respect 
to  each  coordinate  and  evaluate  at  the  limits  of  integration. 


I 

J 
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2.3.1  Logic  Hierarchy 

I'oWnV*  is  divided  into  program  units  call  'phases.*  Each  phase 
performs  a specific  computational  or  system  booKkeeping  task. 

A phase  consists  of  a command  file  and  a function  file.  The 
command  file  is  a collection  of  MACSYMA  commands  and  invocations  of 
system  functions  and  is  executed  via  the  MACSYMA  'BATCH*  command. 
Phases  may,  in  turn,  execute  other  phases  by  issuing  the  appropriate 
BATCH  command.  The  function  file  contains  the  definition  of  the 
system  functions  invoked  in  the  command  file  and  all  of  their  external 
references.  To  minimize  the  amount  of  storage,  the  function  file  is 
loaded  at  the  start  of  the  phase  and  all  functions  are  deleted  before 
the  phase  is  exited.  In  addition,  run  configuation  processes  which 
perform  a trace,  dump  and  break  ( see  section  2.B  ) are  enacted. 

The  general  format  of  a command  file  is: 

(a)  Print  the  header  message  giving  the  decription  of  the  phase, 

(b)  If  the  phase  is  to  be  traced,  turn  on  the  trace  facility, 

(c)  Load  the  function  file, 

(d)  Execute  the  system  functions, 

(e)  If  specified,  perform  a dump, 

(f)  if  specified,  perform  a break, 

(g)  Remove  all  functions  loaded  from  the  function  file. 

2. 3. 1.1  Hierarchy  for  the  Isoparametric  Formulation 

This  section  contains  a listing  of  the  phases,  in  the  order  in 
which  they  are  executed,  for  the  isoparametric  formulation.  Several 
of  the  phases  are  common  to  both  formulations. 

FSTAR  - System  Entry  Point 

The  user  initiates  execution  of  the  system  by  issuing  the  command, 
'BATCH!  [FSTAR, CM0,0SK,AklGl, ON):'.  This  phase  also  opens  the  disk 
file  which  is  to  contain  a record  of  the  execution  and  loads  a set  of 
global  system  functions.  The  phase  then  executes  phases  SYINIT, 
SELECT  and  TERM IN. 

SYINIT  - System  Initialization 

Set  the  operating  environment,  print  opening  messages  and  input  the 
user  identification. 

SELECT  - Method  Selection 

Input  the  user  specification  of  the  formulation  (Isoparametric  or 
generalized  coordinate  ) and  initiate  execution  of  the  appropriate 
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formulation  executive. 


ISOEXC  - Isoparametric  Formulation  Executive 

Execute  phases  associated  uith  the  isoparametric  formulation.  This 
phase  requires  no  user  input. 

ISOEIN  - Isoparametric  Formulation  Initialization 

Process  user  specifications  for  the  run  configuration  (TRACE,  DUMP 
and  BREAK  settings). 


INPISO  - Problem  Parameter  Specification 

Input  and  process  user  specifications  of  the  number  of  element  nodes, 
number  of  degrees  of  freedom  per  node,  number  of  natural  coordinates 
and  the  names  of  the  global  coordinates,  natural  coordinates  and 
displacement  variables. 


SFNISO  - Shape  Function  Processor 

Input  the  shape  f unci  ions  and,  if  necessarg,  express  them  in  terms  of 
an  independent  set  of  natural  coordinates. 


BtIDATA  - B Matrix  Database  Generation 

Process  strain  component  specifications  to  synthesize  the  database  to 
be  used  in  the  generation  of  the  IBl  and  IBal  matrices. 
Specifications  may  be  either  by  library  selection  or  manual  input. 

MATERL  - Material  Properties  Selection 

Process  specification  of  material  properties  by  either  library 
selection  or  manual  input. 


AUXTER  - Auxiliary  Term  Processor 

Input  and  store  the  factors  modifying  the  elemental  volume. 


LIMTIS  - Integration  Limits 

Check  if  the  natural  coordinates  are  dependent  or  independent.  Input 
and  process  limits  accordingly. 

JACOBN  - Jacobian  Generation 

Form  the  Jacobian  matrix  and  its  determinant.  This  phase  requires  no 
input. 

BMXISO  - (BI  Matrix  Generation 

Form  the  (BI  matrix  as  specified  by  the  database.  No  user  input  is 
required. 

INTISO  - Integration  Processor 

See  section  2.2.10  Integration  Processing.  No  user  input  is  required. 
DISIPP  - Display  Pre-processor 

Algebraically  simplify  the  results  of  the  integration  phase  to 
facilitate  evaluation.  User  input  is  not  required. 

DISPLY  - Display  Processor 

Input  user  specification  of  the  output  format  of  the  stiffness 
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matrix,  Proccas  accordingly. 

TERMIN  - Syatem  Termination 

Pr  i nt  ,'C  I os  i ng  messages,  close  the  file  containing  the  record  of  the 
execution  and  process  the  specifications  for  terminating  the  run. 

2.3. 1.2  Hierarchy  for  the  Generalized  Coordinate  Formulation 

This  section  contains  a listing  of  the  phases,  in  the  order  in 
which  they  are  executed,  for  the  generalized  coordinate  formulation. 
For  phases  which  are  common  to  both  formulations  reference  is  made  to 
the  description  in  the  previous  section.  ' 

FSTAR  - System  Entry  Point 
(See  previous  section) 

SYINIT  - System  Initialization 
(See  previous  section) 


SELECT  - Method  Selection 
(See  previous  section) 


GECEXC  - Generalized  Coordinate  Formulation  Executive 

Execute  phases  associated  with  the  generalized  coordinate 

formulation.  No  user  input  is  required. 

GECEIN  - Generalized  Coordinate  Formulation  Initialization 

Process  user  specifications  foe  the  run  configuration  (TRACE,  DUMP 

and  BREAK  settings). 


INPGEC  - Problem  Parameter  Specification 

Input  and  process  the  user  specifications  for  the  number  of  nodes, 
number  of  degrees  of  freedom  per  node  and  the  names  of  the 
displacement  and  global  coordinate  variables. 


SFNGEC  - Shape  Function  Processor 
Input  and  process  shape  functions. 


BMDATA  - IBI  Matrix  Database  Generation 
(See  previous  section) 

MATERL  - Material  Properties  Selection 
(See  previous  section) 

AUKTER  - Auxiliary  Term  Processor 
(See  previous  section) 

LIMTGC  - Integration  Limits 
Determine  the  limits  of  integration. 

No  user  input  is  required. 

AMXINV  - Inversion  of  the  (A)  Matrix 

See  the  algorithm  in  section  2.2.4.  No  user  input  is  required. 
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BI1XGEC  - IBa)  tiatrix  Generation 

. Form  the  IBa)  matrix  as  specified  by  the  database.  No  user  input  is 

o required. 

' INTGEC  - Integration  Processor 

See  the  description  in  section  2.2.10.  No  user  input  is  required. 

□ISGPP  - Display  Pre-processor 

Algebraically  simplify  the  form  of  the  stiffness  matrix  to  facilitate 
{ evaluation.  No  user  input  is  required. 

« i 

DISPLY  - Display  Processor 
(See  the  previous  section) 

TERMIN  - System  Termination 
(See  the  previous  section) 

2. 3. 1.3  Execution  Hierarchy 

Upon  completion  of  a phase,  control  returns  to  the  module  which 
issued  the  BATCH  command  invoking  it.  For  phases  ISOEIN  through 
DISPLY  for  the  isoparametric  method,  the  invoking  module  is  ISOEXC 
which  then  executes  the  BATCH  command  for  the  next  phase  in  that 
sequence.  An  identical  relationship  exists  between  the  sequence  of 
phases  GECEIN  through  DISPLY  with  respect  to  GECEIN  for  the 
generalized  coordinate  formulation.  When  execution  of  these 
executive  phases  is  completed,  control  reverts  to  SELECT.  Phase  FSTAR 
executes  phases  SYINIT,  SELECT  and  TERIIIN. 

This  hierarchy  is  advantageous  in  that  it  facilitates  the 
reorganization  of  existing  phases  and  the  addition  of  new  or 
a I ternat i ve  ones. 


2.3.2  Data  Flow 

Data  Flow  refers  to  the  transfer  of  data  values  among  the 
different  program  segments  which  comprise  the  system.  The 
designation,  ’system  value’,  will  refer  to  those  values  which 
represent  a major  term  in  the  computation  and  which  are  not  local  to 
the  phase  in  which  they  are  created. 

HACSYMA’s  rules  of  scope  are  similar  to  those  implemented  in  ALGDL 
60.  Thus,  system  values  are  automatically  placed  in  a global  pool, 
the  flACSYflA  ’VALUES’  list,  simply  by  not  declaring  then  local  to  any 
function.  Parameters  local  to  a function  appear  in  brackets  at  the 
beginning  of  the  function  body  as  prescribed  by  flACSYflA  syntax. 

To  increase  programming  clarity,  the  following  design  decisions 
relating  to  data  flow  were  made: 


(a)  System  values  are  returned  only  through  functions  explicitly 
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invoked  In  commend  files.  MALSYMA  functions,  like  LISP  functions, 
can  only  return  a single  value.  If  multiple  values  are  to  bo 
generated  by  the  same  function,  they  are  returned  as  elements  In  a 
simple  list  and  unpacked  in  the  command  file. 

(b)  Parameters  of  a function  must  appear  in  the  argument  list. 
Though  the  system  values  are  global  in  scope  and  thus  referencable  in 
the  body  of  any  function,  this  stipulation  uas  made  in  keejjing  with 
good  program  modularity. 

User  specified  values  of  element  volume  modification  factors  and 
other  auxiliary  terms  are  retained  on  a simple  stack.  This  stack  is 
also  used  to  retain  values  resulting  from  internal  computation  (for 
example,  the  determinant  of  the  Jacobian  and  of  the  [A]  matrix). 

The  integration  processor  examines  the  stack  for  expressions  which 
are  functions  of  the  variables  of  integration  and  multiplies  them  into 
the  quadratic  form  (see  section  2.2.10).  The  display  processors  will 
identify  those  elements  in  the  stack  which  are  not  functions  of  these 
integration  variables  and  collect  them  into  a single  expression.  This 
expression  constitutes  a common  factor  of  the  stiffness  matrix. 

A list  of  system  values,  the  phase  in  which  they  are  created 
(designated  in  parentheses)  and  a brief  definition  follows: 

A INVERSE  - (AMXINV) 

The  factored  (A)  matrix,  referred  to  as  (AJ)  in  section  2.2.1 

BMATRIX  - (311X150,  BMXGEC) 
the  (BI  or  IBal  matrix 

BI1ATRIXDATABASE  - (BMDATA) 

Database  containing  the  specification  for  the  synthesis  of  the  (BI  of 
IBa)  matrices  (see  section  2.2.6) 

BREAXPOINTS  - (ISOEIN,  GECEIN) 

List  of  the  phases  for  which  a BREAK  is  requested 
□ELAYSUITCH  - (SYSINT) 

Batch  file  processing  switch,  initialized  to  *0N’,  necessary  to 
overcome  the  Inability  of  MACLISP,  the  base  language  for  MACSYMA,  to 
have  more  than  one  file  open  at  a time. 

DETAMATRIX  - (AMXINV) 

Determinant  of  the  tA)  matrix 

DETJACOBIAN  - (JACOBN) 

Determinant  of  the  Jacobian  matrix 

DUMPPOINTS  - (I see IN,  GECEIN) 

List  of  the  phases  for  which  a DUMP  is  requested 


DISPLACEMENTS  - (INPISO,  INPGEC) 


m 
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List  of  the  nainee  of  the  di  o|3 1 ocemonts 

FACTOR  - (DISIPP,  DISGPP) 

Common  scalar  multiple  of  FUNCTIONAL 

FUNCTIONAL  - (INTISO,  INTGEC) 

Matrix  of  stiffness  coefficients  with  common  multiples  factored 
GLOBALCOORDS  - (INPISO,  INPGEC) 

List  of  the  names  of  the  global  (geometry)  coordinates 
HELP  - (SYINIT) 

String  entered  to  request  the  help  text,  initialized  to  ’?’ 

JACOBIAN  - (JACOBN) 

The  Jacobian  matrix 

LIMITS  - (LIMTIS,  LIMTGC) 

Matrix,  the  rows  of  which  contain  the  upper  and  lower  limits  of 
integration 

MATER lALSMATR IX  - (MATERL) 

Material  properties  matrix  with  common  multiple  factored 

MATERIALSMATRIXMULTIPLIER  - (MATERL) 

Common  multiple  of  MATERIALSMATRIX 

METHOD  - (SELECT) 

Integer  corresponding  to  the  chosen  formulation:  1» i soparamatr i c, 
2=general ized  coordinate 

NATURALCOORDS  - (INPISO,  INPGEC) 

List  of  the  names  of  the  natural  coordinates 

MDOF  - (INPISO,  INPGEC) 

Number  of  degrees  of  freedom  per  node 

NUMNATURAL  - (INPISO,  INPGEC) 

Number  of  natural  coordinates 

NUMNOOES  - (INPISO.  INI’GEC) 

Number  of  element  nodes 

NUMSTRA INCOMPONENTS  - (BMDATA) 

Number  of  strain  components 

NUMSTRAINTERMS  - (BMDATA) 

Total  number  of  strain  terms;  a strain  component  may  contain  more 
than  one  algebraic  term 
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SHAPEFUNCTI0N5  - (SFNISO.  SFNGEC) 

Mntrix,  the  nous  of  uhich  contain  the  shape  functions 
STACK/-  (SYIWIT,  AUKTER) 

Push-down  stack  containing  the  element  volume  modification  factors, 
and  auKiliory  terms 

STACKPOl NITER  - (SYINIT) 

Pointer  to  just  below  the  top  element  in  STACK 
TRACEPOINTS  - (ISOEIN,  GECEIN) 

List  of  the  phases  for  which  a TRACE  is  requested 
USERNAME  - (SYINIT) 

User  identification;  used  to  name  the  Record  File 


2.4  Help  Text 

Anticipating  that  diagnostic  messages  and  prompts  may  not  provide 
sufficient  information  about  a particular  input  value,  a mechanism 
providing  on-line  assistance  is  incorporated. 

The  user  may  answer  any  prompt  with  "?"  and  a brief  help  text  will 
be  printed  at  the  terminal.  The  text  contains  a description  of  the 
requested  input  and,  if  a selection  Is  to  be  made,  an  enumeration  of 
the  a I ternat i ves. 

Upon  exiting  the  help  text,  control  returns  to  the  prompt  for  the 
i nput. 


2.5  Input  and  Output  Functions 


2. 5. 1 Record  File 

To  assist  in  diagnostics  and  system  development,  a record  of  each 
execution  of  the  system  is  made.  The  record  is  written  to  a disk  file 
and  contains  a copy  of  all  user  inputs  command  executions  and  system 
dumps. 


2.5.2  Input  and  Output  Facilities 

In  the  default  mode,  MACSYMA  commands  are  echoed  at  the  terminal 
as  they  are  typed  in  if  entered  interactively  or  as  they  are  executed 
if  prestored  on  a disk  file.  A command  terminated  with  a semicolon 
will  display  the  returned  value,  while  if  terminated  with  a dollar 
sign,  will  not.  All  output  to  the  terminal,  excluding  the  echo  of 
typed  input  is  suppressed  if  the  switch,  ’TTYOFF’  is  set  to  TRUE. 
This  condition  prevails  until  the  switch  is  reset  to  FALSE. 
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As  in'ritin’t  is  intended  to  be  an  interactive  system,  it  is  necessary 
to  selectively  suppress  output  to  the  user’s  terminal:  the  echo  of 
commands  and  the  automatic  display  of  the  results  retrieved  by 
functions  is  prevented  but  prompts,  help  text,  error  messages  and 
requested  values  are  displayed.  In  addition,  a complete  record  of  user 
input,  all  command  executions  and  values  returned  by  command  level 
functions  is  placed  on  the  file  containing  the  record  of  the 
execution.  To  accomplish-  this,  the  following  steps  were  implemented: 

(a)  All  function  calls  in  the  command  files  are  terminated  with  a 
dollar  sign  thus  suppressing  the  automatic  display  of  returned  values. 

(b)  The  TTYOFF  switch  is  set  to  TRUE  in  phase  FSTAR,  System  output 
functions  locally  set  TTYOFF  to  FALSE,  execute  flACSYriA  display 
commands  and  then  reset  TTYOFF.  Thus  all  output  is  suppressed  except^ 
that  which  is  channeled  through  the  system  output  functions. 

(c)  Since  TTYOFF  only  affects  output  to  the  user's  terminal,  the  file 
containing  the  record  of  the  execution  is  unaffected. 

2.5.3  Input  Characteristics 

The  use  of  the  MACSYUA  ’REAO’  function  presents  some  difficulties. 
The  function  will  read  in  and  evaluate  one  expression.  If  the  input 
corresponds  to  the  name  of  a variable  local  to  the  function  containing 
the  READ  or  one  in  the  global  pool,  the  evaluation  step  causes  the 
input  expression  to  be  replaced  by  the  value  of  that  variable.  The 
resolution  was  provided,  upon  request  from  HACSYMA  system  programmers 
in  the  form  of  the  heretofore  undocumented  ’READIN’  function.  The 
specification  for  the  READIN  function  is  the  same  as  that  for  READ 
with  the  exception  that  the  input  expression  is  not  evaluated. 

All  input  is  format  free  with  the  exception  that  expressions  must 
be  terminated  with  either  a semicolon  or  a dollar  sign.  It  is  not 
necessary  to  press  the  RETURN  key.  Either  upper  or  lower  case  may  be 
used. 


2.6  Run  Configuration 

Run  Configuration  refers  to  the  state  of  the  TRACE,  BREAK  and  DUMP 
facilities  which  may  be  set  for  individual  phases.  These  features  are 
intended  primarily  for  system  maintenance. 


The 

command 

TRACE  facility  causes  an 
file  as  they  are  executed. 

echoing  of 

the  commands 

I n 

the 

If 

the  BREAK  switch  is  set 

for  a phase. 

execut i on 

will 

be 

suspended  just  prior  to  the  phase  termination  and  control  will  revert 
to  MACSYMA  command  level.  Execution  resumes  upon  entry  of  "EXIT;”. 


Ulhereas  TRACE  and  BREAK  processes  are  intended  purely  for  system 
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mn i nton.incc),  the  IJUflP  feature  may  be  valuable  to  the  user.  If  UUnP  Is 
set,  the  system  values  generated  during  a phase  are  displayed  at  the 
user’s  console.  Independent  of  this  setting,  these  values  are  always 
dumped  to  the  file  containing  the  record  of  the  execution.  As  certain 
dumps  are  quite  lenghty,  this  feature  should  be  used  judiciously. 

The  run  configuration  may  be  set  in  the  phases  ISOEIN  and  GECEIN. 


2.7  Error  Recovery 

System  error  recovery  emcompasses'  three  processes; 

(1)  Semantic  errors  are  trapped  during  error  checking,  by  system 
functions.  A diagnostic  is  printed  and  the  user  is  prompted  to 
re-enter.  If  the  diagnostic  offers  insufficient  explanation,  the  user 
may  • answer  the  prompt  with  which  is  the  request  for  the  help 
text. 

(2)  Syntactic  errors  are  trapped  by  MACSYMA.  The  erroneous  input  is 
echoed  along  with  a pointer  indicating  the  subexpression  in  error.  A 
diagnostic  is  printed  and  the  message,  "Please  rephrase  or  edit,"  is 
displayed.  The  user  may  either  retype  the  input  or  press  the  escape 
key  to  enter  the  HACSYIIA  editor.  The  editor,  a derivative  of  TECO,  is 
documented  in  Bogenf7I. 

(3)  Errors  which  cannot  be  trapped  by  either  of  the  above  procedures 
are  associated  with  a system  failure  and  are  fatal.  The  system  will 
attempt  to  shutdown  as  neatly  as  possible  by  preserving  the  global 
pool  of  values  and  closing  the  file  containing  the  record  of  execution 
so  that  diagnostics  may  be  run.  A message  directing  the  user  to 
sources  of  assistance  is  displayed  and  the  user  Is  queried  for  the 
system  termination  procedure  as  in  a normal  run. 


2.8  Limi tations 


The  limitations  of  the  present  system  are; 

(a)  Each  node  must  have  the  same  degrees  of  freedom; 

(b)  Each  degree  of  freedom  must  be  represented  by  the  same  shape 
function; 

(c)  Allowable  forms  for  the  specification  of  strain  components  are; 

a Vi  D(u,x) 

a Vi  D(u,x)  + b Vr  Dlu,x) 
u/a 

where  ’a’  and  ’b’  are  scalars,  'u*  represents  a displacement  variable 
and  ’x’  represents  a geometry  (global  coordinate)  variable, 

(d)  Problem  size  is  limited  by  available  list  space  on  the  computer 
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suppor  t i ng  MACSYMA . 

Minor  modification  uill  permit  the  processing  of  second 
derivatives  andprov i s i ons  in  the  program  have  been  incorporated  to 
facilitate  the  elimination  of  restrictions  (a)  and  (b) . 


2.3  I mp I emontat i on  of  Element  Loads,  Mass  and  Stresses 

The  present  version  of  the  processor  mag  be  extended  to  permit  the 
implementation  of  the  quantities  defined  in  section  1.2,3. 

I 

Many  of  the  constituent  matrices  are  formulated  as  part  of  the 
synthesis  of  the  stiffness  matrix.  These  include  the  IBI  and  tCl 
matrices  and  the  matrix  of  shape  functions.  The  remaining  constituent 
matrices  may  be  taken  from  user  specification  of  [ml,  (fbl , Iftl,  lei) 
and  the  names  of  the  nodal  displacements. 

Integration  functions  can  be  used  to  evaluate  the  required  volume 
and  surface  integrals. 
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diopter  3 

ILLUSTRATIVt  EXAUPLES 


This  chapter  presents  a number  of  illustrative  examples  processed 
on  iv-.v^nvvr.  The  bulk  of  the  presentation  consists  of  the  execution 
record.  Comments  are  placed  in  angle  brackets  ( ">"  ) and  mere 

not  part  of  the  execution. 

The  internal  storage  scheme  for  the  stiffness  matrix  is  a vector, 
"FUNCTIOIOAl. ",  containing  the  elements  of  the  louer  triangular  portion 
in  roil  major  format  and  a variable,  "FACTOR",  containing  the  common 
multi.ple  of  each  element.  In  the  verification  sections  for  each 
example,  these  c\uantities  are  retrieved  and  MACSYUA  functions  are  used 
to  recast  the  results  into  alternate  forms  and  to  establish 
i dent i ties. 

HACSYMA  functions  emploijed  in  the  verification  sections  include 
EV,  uhich  evaluates  its  first  argument  in  the  environment  specified  liy 
the  remaning  arguments,  and  FACTOR,  uhich  factors  its  argument  into 
factors  irreduciljle  over  integers,  (The  function,  FACTOR,  should  not 
be  confused  with  the  -itirMn’t  variable,  FACTOR.  They  are  I'eadily 
di  st i ngui sabl e in  that  the  function  requires  an  argument  in 
parenthesis).  Details  of  these  and  the  other  functions  may  be  found 
in  reference  171. 

The  examples  ore  presented  as  follows: 

Section  3.1  - Bar  Element  with  Uniform  Cross-Sectional  Area 
Part  3.1.1  - Genera  I i::ed  Coordinate  Formulation 
Part  3.1.2  - Isoparametric  Formulation 
Part  3.1,3  - Comparison  and  Verification 

Section  3.2  - Bar  Element  with  Linear  Variation  of  Cross-Sectional 
Area 

Port  3.2.1  - Generalized  Coordinate  Formulation 
Part  3.2.2  - Isoparametric  Formulation 
Part  3.2.3  - Comparison  and  Verification 

Section  3.3  --  Constant  Strain  Triangle  (CST)  with  Uniform  Thickness 
Part  3.3.1  - Generalized  Coordinate  Formulation 
Part  3.3.2  - Isoparametric  Formulation 
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I he  fine'll  section  is  devoted  to  a discussion  of  the  a I yebra  i c form 
of  the  results. 


I 
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3.1  Bar  Eirmcr.  with  Uniform  Crosr.-Scct  1 ona I Area 

<Thq  /bar  einment  has  one  cJccirne  of  freedom,  in  the  axial  direction, 
and  i B modeled  usintj  a linear  dioplocement  function. > 

3.1.1  Generalized  Coordinate  Forniuiation 

<I  mmedi ate  I y subsequent  to  the  LOGIN  procedure,  a copy  of  flACSYMA  is 
loaded  ( via  ";A")  and  the  BATCH  command  is  issued  to  initiate 
execut  ion  of  »ViVjnV*.  > 

*:  A 

This  is  riACSYnA  26B 

F1X26G  8 D5K  IIACSYMA  being  loaded 

I oad i ng  done 

(Cl)  BATCH ( [FSTAR.CnO,DSK,AKlG]  ,0N); 

(C2)  TTY0FF:TRUE11 

* SYSTEft  INITIALIZATION  * 

IJiiLCOflli  TO  *****  VERSION  1.0 

It  is  now  THURSDAY  DECEUBER  1,1977  18:2:10 

Tlie  current  file  is  ISYINIT,  FCNI 

The  current  device  and  username  is  [OSK,  AKIGI 

Report  problems  to 
At. AN  R.  KORNCOFF 
DEPT.  OF  CIVIL  ENGINEERING 
CARNEG I E-MELLON  UNIVERSITY 
CnU-10A,  AlClG 

Terminate  all  input  with  a SEMICOLON  - 
Input  for  HELP 


Input  your  LOGIN  NAME 
? 

AKIG; 

Is  AKIG  correct  ? (YES;  or  NO;) 

? 

YES; 

GREETINGS  AKIG 


* METHOD  SELECTION  * 


ll.I.lJ'URAdVH  EXAIin.ES 

enter  tlie  number  of  the  method  choeen  (1  OR  2). 

? 

GENERALIZED  COORDINATE  FORHULATION 
* GENERALIZED  COORDINATE  FORfUJLATION  EXECUTIVE  -.v 
Vf  GENERALIZED  COORD 1 NATE  FORMULATION  INITIALIZATION  * 


Do  you  wish  to  set  DUMP,  BREAK  or  TRACE  POINTS 
? 

YES; 


1.  GECEIN 

2.  INPGEC 

3.  SFNGEC 

4. 

BMDATA 

S.  MATERl. 

G.  AUXTER 

7.  LIMTGC 

S. 

AMXINV 

9.  BMKGCC 

10.  INTGEC 

11.  DI5GPP 

12.  DISPl.Y 

SET  DUMP  POINTS 

Des i gnate 

selected  PHASES 

by  entering 

ttie  associated  integer  INDEX  or  'ALL'  for 

a 1 

1 pliases. 

Type  ’END’ 
? 

-7 , 

to  TERMINATE. 

? 

4; 

? 

5; 

? 

END; 

1.  GECEIN 

2.  INPGEC 

3.  SFNGEC 

4. 

DI1DATA 

5.  MATERL 

G.  AUXTER 

7.  LIMTGC 

8. 

AMXINV 

3.  BI1KGEC 

10.  INTGEC 

11.  DISGPP 

12.  DISPLY 

SET  BREAKPOINTS 

Dos i gnate 

selected  PHASES  by  entering 

llie  associated  integer  INDEX  or  ’ALL’  for 

a 1 

1 phases. 

Type  ’END’ 

9 

to  TERMINATE. 

END; 

1.  GECEIN 

2.  INPGEC 

3.  SFNGEC 

4. 

RMDATA 

5.  MATERL 

G.  AUXTER 

7.  LIMTGC 

8. 

AMXINV 

3.  BMKGEC 

10.  INTGEC 

11.  DISGPP 

12.  DISPLY 

SET  TRACE 

POINTS 

Des i gnate 

selected  PHASES 

by  entering 

the  associated  integer  INOEX  or  ’ALL’  for  all  phases. 

Type  ’END’  to  TERMINATE. 

? 

END; 

>v  PROBLEM  PARAMETER  SPECIFICATION  - GENERAL  I ZEO  COORDINATE  * 


J 

J 


Input  the  NUMBER  OF  ELEMENT  NODES 
? 
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Input  tb(*  NUMIJER  OF  OCGREES  OF  FREEDOfI  PER  NODE 
? 

1;  ' ' 

Input  thR  vector  of  tltn  NAMES  OF  THE  GLOBAl.  COORDINATES 
Ihere  shoulti  be  1 element 6 
ELEMENT  ].  « 

X; 

Input  the  vector  of  the  NAMES 'OF  THE  DISPLACEMENT  VARIABLES 
There  should  l.ie  1 elements 
ELEMENT  1 = 

U; 

DUMP-  FOR  I NPGEC 

NIJMNODES  « 2 
NOOF  = 1 

GLOBAL COORDS  = IX] 
DISPLACEMENTS  » [U] 


••v  SHAPE  FUNCTION  PROCESSOR  - GENERALIZED  COORDINATE  Vr 

ENTER  the  terms  of  the  SHAPE  FUNCTION  ordered  from 
GENERALIZED  COORDINATE  1 through  coordinate  2. 

ELEMENT  1 » 

1; 

ELEMENT  2 - 
X; 

< This  input  represents  the  displacement  function:  u = al  + aZvrx  , in 
uhich  "al"  and  "a2"  are  the , genera  I i red  coord i nates. > 

SHAPE  FUNCTION  MODIFICATION 
The  OPTIONS  are 

(1)  DISPLAY  THE  SHAPE  FUNCTIONS 

(2)  MODIFY  THE  SHAPE  FUNCTIONS 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

1; 


The  terms  of  the  SHAPE  FUNCTION  arc 
Term  1 
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» 

X 

Enter  thR  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

3; 

* B MATRIX  DATABASE  GENERATION  * 

The  OPTIONS  for  f.|3cci  fijinc)  STRAIN  COMPONENTS  are 

(1)  USER-SUPPLIED  VALUES 

(2)  LIBRARY  VALUES 

Please  ENTER  the  NUMBER  ASSOCIATED  WITH  YOUR  SELECTION 
? 

2; 

The  LIBRARY  OPTIONS  for  SPECIFYING  STRAIN  COMPONENTS  are 

(1)  USER-SUPPLIED  VALUES 

(2)  ONE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 

(4)  PLANE  STRAIN 

(5)  AX I SYMMETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 30 
Please  ENTER  the  NUMBER  ASSOCIATED  WITH  YOUR  CHOICE 
? 

2; 

< The  library  contains  pre-stored  strain  components  specifications  in 
the  database  format  described  in  section  2.2.5.> 

DUMP  FOR  3M0ATA 

BI1ATRIXDATABASE  =[111101 
NUMSTRAINTERMS  - 1 
NUMSTRA INCOMPONENTS  » 1 


* MATERIAL  PROPERTIES  SELECTION  * 

The  options  for  the  selection  of  the  MATERIAL  PROPERTIES  MATRIX  are 

(1)  USER-SUPPLIED  MATRIX 

(2)  LIBRARY  MATRIX 

Please  enter  the  NUMBER  ASSOCIATED  WITH  YOUR  SELECTION 
? 

2; 


The  LIBRARY  OPTIONS  for  SPECIFYING  MATERIAL  PROPERTIES  are 
(1)  USER-SUPPLIED  VALUES 
(?)  ONE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 

(4)  PLANE  STRAIN 

(5)  AXI SYMMETRIC 

(6)  LINEAR  ISOTROPIC  ELASTICITY  - 3D 
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Pln<.!SG  PNITR  tho  NUtlBER  ASSOCIATED  UlTH  YOUR  CHOICE 
? 


DUflf’  FOR  MATERL 

riATERlALSriATRlX  =11] 
nATERlALSriATRlXniJLTl  FLIER  = E 


>v  AUXILIARY  TERM  PROCESSOR  * 

ENTER  ELEMENT  VOLUME  MODIFICATION  FACTORS  OR  AUXILIARY  TERMS 
TYPE  ’END;’  to  TERMINATE 
? 

A; 

? 

END; 

AUXILIARY  TERM  MODIFICATION 
The  OPTIONS  are 

(1)  DISPLAY  THE  AUXILIARY  TERMS 

(2)  MODIFY  AN  AUXILIARY  TERM 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

,v  INTEGRATION  LIMITS  - GENERALIZED  COORDINATE  * 

,v  INVERSION  OF  MATRIX.  A * 

■.V  B MATRIX  GENERATION  - GENERALIZED  COORDINATE  * 

* INTEGRATION  PROCESSOR  - GENERALIZED  COORDINATE  * 

Vf  DISPLAY  PRE-PROCESSOR  - GENERALIZED  COORDINATE  * 

* DISPLAY  PROCESSOR  * 

The  options  for  OUTPUTTING  the  STIFFNESS  MATRIX  are 

(1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARD  IMAGE  format 

ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 

i 


Enter  the  FORTRAN  ARRAY  NAME  for  the  STIFFNESS  MATRIX 
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Hultiply  Ctnch  coefficient  bij 

-A*E/(XI2)-.X(l))*,-f2 

The  STIFPNE5S  I1ATRIX  is 

ST  IFF  (1,1)  = -(XHIN-^XnAX) 
STIFF (2,1)  « XniN-XhAX 
ST  IFF  (3,1)  = -(XIIIN-XnAX) 


* SYSTEM  TERMINATION 

Thfit’c  all  for  nou,  AKIG 
It  is  now  12/1/77  18:13:57 

AccuhuilatecI  CPU  TIME  = 18348  MSEC 

A record  of  this  session  is  recorded  in  file  lAKlG,  >,  D5K, 

The  options  for  cxitinq  the  system  arc 

(1)  TERMINATE  THE  RUN 

(2)  TERMINATE  THE  RUN  AND  THE  JOB 

(3)  TERMINATE  THE  RUN,  THE  JOB  AND  LOGOUT 

Please  enter  the  number  of  the  option  chosen  (1,  2 or  3) 

? 

3: 

Accumulated  CPU  TIME  = 18642  MSEC 

TERMINATE  AND  LOGOUT 

.•LOGOUT 

<The  system  automatically  issued  the  ":L0G0UT"  command, 
the  termination  procedure  are  given  in  the  Appendix. > 

<Option  1 of  the  DISPLAY  PROCESSOR  mould  provitJe  the 
output: > 


The  options  for  OUTPUTTING  the  STIFFNESS  MATRIX  are 

(1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARD  IMAGE  format 

ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 

1; 

Multiply  each  coefficient  by 

A E 


2 

(X  - X ) 

2 1 


USERS  ] 


Details  of 

f o I lowing 


ROIJ  2,  COL  2 


- (XniN  - XHAX) 


3.1.2  Isoparametric  Formulation 


* SYSTEM  INITIALIZATION  * 

UELCOME  TO  **;•.**  VERSION  1.0 

It  is  now  THURSDAY  DECEMBER  1,1977  1:41:25 

The  current  file  is  ISYINIT,  FCNl 

The  current  device  and  username  is  (DSK,  AKIGI 

Report  problems  to 
Al.AN  R.  KORNCOFF 
DEPT.  OF  CIVIL  ENGINEERING 
CARNEGIE-NELLON  UNIVERSI  TY 
CMU-10A.  AKIG 

Terminate  all  input  with  a SEMICOLON  - 
Input  ’?;•  for  HELP 

Input  your  LOGIN  NAME 
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(2)  FHl;  GlINCRALlZrO  COORDINATE  111. THOU 

Plttiise  enter  the  nuwljor  of  the  utothocl  cliojipn  (1  OR  2). 

? 

1 ; 

1 OOPARAni;  TR 1 C FORnUL  AT  I ON 
,v  ISOPAMAUCTRIC  FORflULATlON  EXECUTIVE  * 

* ISOPARAHMETRIC  FORtUJLATION  INITIALIZATION  * 

Oo  you  wish  to  set  DUMP,  BREAK  or  TRACE  POINTS 
7 

NO; 

* PROBLEM  PARAMETER  SPECIFICATION  - ISOPARAMETRIC  * 

Input  the  Nl.IMPER  OF  ELEMENT  NDOES 

7 

Input  the  NUMBER  OF  OEGREES  OF  FREEOOM  PER  NODE 

7 

1; 

Input  the  NUMBER  OF  NATURAL  COORDINATES 

7 

i; 

Input  the  vector  of  the  NAMES  OF  THE  NATURAL  COORDINATES 
There  should  be  1 elements 
ELEMENT  1 « 

L; 

Input  the  vector  of  the  NAMES  OF  THE  GLOBAL  COORDINATES 
There  should  be  1 elements 
ELEMENT  1 » 

X; 

Input  the  vector  of  the  NAMES  OF  THE  DISPLACEMENT  VARIABLES 
There  should  be  1 elements 
ELEMENT  1 - 
U; 

* SHAPE  FUNCTION  PROCESSOR  - ISOPARAMETRIC  >v 

ENTER  the  terms  of  the  SHAPE  FUNCTION  ordered  from 
node  1 through  node  2. 

The  2 elements  will  be  prompted  for 
ELEMENT  1 » 

.5*(I-L); 

ELEMENT  2 - 
.S,vll-*L); 


* 
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< Tliin  input  rttprc’sents  the  d i np  I nccment  ec|(i.it  ion; 
u " 0.5»v(l-l)  + 0.5»v(l  + l)  .> 

SHAPE  FUNCTION  HOD I FI CATION 
The  OPTIONS  are 

<1)  DISPLAY  THE  SHAPE  FUNCTIONS 
(21  f1(iniFY  THE  SHAPE  FUfCTIONS 
(3)  TERtHNATE  THIS  FUNCTION 
Enter  the  NUflBER  ASSOCIATED  with  the  CHOSEN  OPTION 


The  terms  of  the  SHAPE  FUNCTION  are 
Term  1 

0.5*(1  - L) 


Enter  the  NUHBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

3; 

,v  B MATRIX  DATABASE  GENERATION  * 


The  OPTIONS  for  specifying  STRAIN  COMPONENTS  arc 

(1)  USER-SUPPLIED  VALUES 

(2)  LIBRARY  VALUES 

Please  ENTER  the  NUMBER  ASSOCIATED  WITH  YOUR  SELECTION 


The  LIBRARY  OPTIONS  for  SPECIFYING  STRAIN  COMPONENTS  are 

(1)  USER-SUPPLIED  VALUES 

(2)  DNE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 
(A)  PLANE  STRAIN 
(S)  AXI SYMMETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 3D 
Please  ENTER  the  NUMBER  ASSOCIATED  UITH  YOUR  CHOICE 
? 

2; 


* MATERIAL  PROPERTIES  SELECTION  * 


The  options  for  the  selection  of  the  MATERIAL  PROPERTIES  MATRIX  are 

(1)  USER-SUPPLIED  MATRIX 

(2)  LIBRARY  MATRIX 

Please  enter  the  NUMBER  ASSOCIATED  UITH  YOUR  SELECTION 

7 
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lh«  LIDNARY  DPI  IONS  for  SPECIFY! MG  IIATERIAl.  PHOPERTItS  or« 
11)  USER-SUPPLIED  VALUES 

(2) . ONE  DIUHNSIONAl.  ELASTICITY 

(3)  PLANE  STRESS 
(A)  I'^LANE  STRAIN 
(5)  AXISYMI1ETRIC 

(G1  LINEAR  ISOTROPIC  ELASTICITY  - 30 
Pl«?cisc  ENTER  tho  NUMBER  ASSOCIATED  WITH  YOUR  CHOICE 
? 

2; 


* AUXILIARY  TERM  PROCESSOR  * 

ENTER  ELEMENT  VOLUME  MODIFICATION  FACTORS  OR  AUXILIARY  TERMS 
TYPE  ’ENO;’  to  TERMINATE 
? 

A: 

? 

END; 


AUXILIARY  TERM  MODIFICATION 
The  OPTIONS  cire 

(1)  DISPLAY  THE  AUXILIARY  TERMS 

(2)  MODIFY  AN  AUXILIARY  TERM 

(3)  TERMfNATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

3; 

,v  INTEGRATION  LIMITS  - ISOPARAMETRIC  * 

ENTER  the  LIMITS  OF  INTEGRATION  for  NATURAL  COORDINATE,  L 
LOUER  LIMIT  - 
7 

-1? 

UPPER  LIMIT  - 
? 

1; 


.V  JACOBIAN  GENERATION  - ISOPARAMETRIC  * 

* B MATRIX  GENERATION  - ISOPARAMETRIC  * 

* INTEGRATION  PROCESSING  - ISOPARAMETRIC  * 

* DISPLAY  PROCESSOR  * 

The  options  for  OUTPUTTING  the  STIFFNESS  MATRIX  ore 

(1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARO  IMAGE  format 
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ENIFiR  thp  NUriBER  ASSOCIATED  with  YOUR  SEl.I-CUON. 
? 

1; 

/ 

hultiply  each  coefficient  by 

2*A*E 

X - X 
2 1 

ROW  1.  COL  1 

1 

2 

ROW  .2,  COL  1 

1 

2 

ROW  2,  COL  2 

1 

2 

* SYSTEM  TERMINATION  * 

That’s  all  for  now,  AKIG 
It  is  now  12/1/77  1:45:24 


Accumulated  CPU  TIME  - 17040  MSEC 

A rocortl  of  this  Gcsnion  io  rncordeci  in  file  (AKIG,  >,  DSK,  USERS  ] 
<DiBplny  option  2 produces: > 

The  options  for  OUTPUTTING  the  STIFFNESS  MATRIX  are 

(1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARD  IMAGE  format 

ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 


Enter  the  FORTRAN  ARRAY  NAME  for  the  STIFFNESS  MATRIX 
? 

STIFF; 


1 


ll.LUSTKATIVC  LKAIII-LtS 


Hultiplij  eijch  coefficient  by 

, /2*A-,vE/(X(2)-K(1)) 

The  STIFFNESS  FIATRIX  is 

STIFF(l.l)  = 1/2 
STIFF(2,1)  = -1/2 
STIFF(3,1)  = 1/2 


3.1.3  Cowp-nr  i fion  one!  Ver  i f i ciit  i on 

<The  files.  C/IG.voIg)  and  (^ll.vals)  contain  the  value?  generated 
during  tlie  runs  in  sections  3.1.1  and  3.1.2  respectively.  For  the 
purposes  of  comparison,  for  the  gener'aliced  coordinate  formulation, 
tl'ie  cooi'dinate,  "X",  is  taken  to  be  a rnininuim  at  node  1 and  a m.Tximum 
at  node  2 ( i.e.  x(l)<x(2)  ),  > 

liine«  h msec. 

(01)  [OSK.  AKIGI 

(C2)  L0A0FILE(’>ill,VALS)5 

%11  VALS  DSK  AKIG  being  loaded 
loading  done 
ti  tl’ie--  291  msec. 

(02)  DONE 

(C3)  FUICTIONAI.; 
time"  0 msec. 

[ 1 } 

( - ] 

( 2 ] 

( ] 

[ 1 ] 

(03)  [ - - I 

I 2 ] 

[ ] 

[ 1 1 
[ - ) 

[ 2 ) 


(CA)  FTL;FUNCTIONAL*FACTOR;  ! 

tiinr?**  lb  mr^ec.  \ 


II.UJSTr'IATIVfc;  IZXAI1I1ES 


(DA)  [ ] 

[ 2 (X  - X ) ] 

I 2 1] 

, t ] 

[ A E ] 

t ] 

[ X - X ] 
•[21] 


(C5)  FACTOR (%); 
time--  78  msec. 

[ A E ] 

[ — ] 

[ X - X ] 

[2  1 ] 

[ ] 

[ A E ] 

(CIS)  [ ] 

[ X - X ] 
[ 2 1 ] 
t ] 

t A E ] 

[ ] 

[ X ~ X ] 

C 2 1 ] 


< This  is  the  lower  triangular  portion  of  the  stiffness  matrix  for  the 
isoparametric  formulation  expressed  in  row  major  format.  > 

(CB)  LOADFILEt’/lG.VALS); 

%1G  VAl.S  DJ5K  AKIG  being  loaded 
I oad i ng  done 
timR«  258  msec. 

(DG)  DONE 

(C7)  FTL2:FUNCT10NAL*FACT0R; 
time"  15  msec. 

t A E [XniN  - XHAX)  ] 

[ ] 

t 2 ] 

t (X  - X ) ] 

t 2 1] 

t ] 

[ A E (XniN  - XtlAX)  ] 

[ I 

(07)  [ 2 ] 

[ (X  - X ) ] 

[21  ] 

[ ] 

[ A E (XnJN  - XflAX)  ] 

[ ] 

[ 2 ] 


1 
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(X  - K ) 
2 1 


< The,/'  XHIN  and  XMAX  terms  of  the  general  cred  coordinate  stiffness 
matrix  in  expression  (D7)  are  evaluated  at  X(l)  and  X(2)  respectively 
and  the  vector  is  simplified.  > 


(C8)  FACTOR  (EV  (%.  XMIN-X  [1} , XI1AX--.X  (2] ) ) : 
timc=  123  msec. 


( X 
[ 2 


X - X ] 
2 1 ] 
] 

A E ] 


[ X - X ] 

[2  1 ] 


< This  is  the  loner  triangular  portion  of  the  stiffness  matix  for  the 
generalised  coordinate  formulation  in  rou  major  format,  evaluated  at 
the  boundary  conditions.  It  is  identical  to  the  expression  in  (05) 
for  the  isoparametric  formulation.  > 


(C3)  CL0SEFILE(%VER,C1); 


Il.l.Ur^TRATlVC  lIXAIIPLFS 


3.2  Bar  EleniRnt  uith  Linear  Variation  of  Cross-Sect  i ona  I Area 
3. 2.1  Genera  I i zed  Coordinate  Formulation 

< This  run  was  identical  to  that  of  section  3.1.1  with  the  exception 
that  the  cross-sectional  area  of  the  Par  varies  linearly  from  a value 
of  Ad)  at  node  1 to  a value  of  A(2)  at  node  2.  The  input  unique  to 
this  execution  is  ohoun  be  low,  > 


* AUXILIARY  TERM  PROCESSOR  * 

ENTER  ELEMENT  VOLUME  MODIFICATION  FACTORS  OR  AUXILIARY  TERMS 
TYPE  ’END;’  to  TERMINATE 
? 

(X(2)~X)*Ad)/(X(2)-X(l))  + (X-X(1))*A(2)/(X(2)-X(1)); 

? 

END; 

AUXILIARY  TERM  MODIFICATION 
The  OPTIONS  are 

(1)  DISPLAY  THE  AUXILIARY  TERMS 

(2)  MODIFY  AN  AUXILIARY  TERM 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  uith  the  CHOSEN  OPTION 
? 

1; 

The  AUXILIARY  TERMS  are 
Term  1 

A(2)*(X  - X(l))  A(1)*(X(2)  - X) 


X(2)  - X(l)  X(2)  - X(l) 

Enter  the  NUMBER  AFaOCIATED  with  the  CHOSEN  OPTION 
? 

3; 


< Both  options  of  the  display  are:  > 


The  options  for  OUTPUTTING  the  STIFFNESS  MATRIX  ore 

(1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARD  IMAGE  format 

ENTER  the  NUMBER  ASSOCIATED  uith  YOUR  SELECTION. 

? 

1; 

Multiply  each  coefficient  by 

E 
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ILLUSTRATIVl-;  EXAIH-LES 

2 

(X  - X ) 

2 1 


ROU  1.  COL  1 

- (XI1IN  - XHAX)  (A  (2)  XniN  - All)  XIIIN  + A (2)  XMAX  - Ad)  XHAX 

+ 2 Ad)  X - 2 X A(2)  )/(2  (X  - X ) ) 
2 1 2 1 


ROU  2.  COL  1 

(XniN  - Xt1AX)  (A(2)  XniN  - Ad)  XIIIN  + A(2)  K(1AX  - Ad)  XIIAX 

+ 2 Ad)  X - 2 X A(2)  )/(2  (X  - X ) ) 
2 1 2 1 


ROW  2,  COL  2 

- (XIIIN  - XMAX)  (A(2)  XniN  - Ad)  XIIIN  + A(2)  XHAX  - Ad)  XMAX 

♦ 

+ 2 Ad)  X - 2 X A(2))/(2  (X  - X )) 
2 1 2 1 


Tite  options  for  OUTPUTTING  tl'ie  STIFFNESS  MATRIX  are 

(1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARD  IMAGE  format 

ENTER  tl-ie  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 

2; 

Enter  tl.e  FORTRAN  ARRAY  NAME  for  the  STIFFNESS  MATRIX 
? 

STIFF; 

Multiply  each  coefficient  by 
E/(X(2)-X(l))vr*2 
The  STIFFNESS  MATRIX  is 

STI FF  (1 . 1 ) - - (XMI N-XMAX)  ,v (A  (2)  >vXMI N-A  (1 ) *)(MI N+A  (2)  v.XMAX-A  ( 1 ) *XMAX 
1 +2-.VA  ( 1 ) *X  (2)  -2*X  ( 1 ) >vA  (2) ) / (2*  (X  (2)  -X  ( 1 ) ) ) 

STI  FF  (2, 1 ) - (XMI  N-XMAX)  -.v  (A  (2)  *XMI  N-A  ( I ) -.vXMl  N+A  (2)  *XMAX-A  (1 ) *XMAX+ 
1 2*A  ( 1 ) »vX  (2)  -2*X  ( 1 ) *A  (2) ) / (2*  (X  (2)  -X  (1 ) ) ) 

ST  I FF  (3, 1 ) - - (XMI  N-XMAX)  * (A  (2)  *XMIN-A  (1 ) *XM)  N+A  (2)  *XMAX-A  ( 1 ) *XMAX 


m.liSTIMTJVf-:  l-KAnPLES 

1 +2*A  (1 ) *K  (2)  -2*X  ( 1 ) *A (2) ) / (2*  (X  (2)  -X  ( 1) ) ) 


/ 

3.2.2  1 sopor amptr i c Formulation 

< The  input  for  this  case  is  the  same  as  that  for  the  genera  I i zed 
coordinate  formulation.  The  display  options  are:  > 

The  options  for  OUTF’IJTTING  the  STIFFNFSS  flATRlX  are 

(1)  Upper  triangular  portion  in  alychraic  format 

(2)  FORTRAN  CARD  II1AGE  format 

ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 

1; 

Multiply  each  coefficient  by 

2 E 


X - X 
2 1 


ROW  1,  COL  1 

A(2)  + Ad) 


4 


ROW  2,  COL  1 

AI2)  + Ad) 


4 


ROW  2.  COL  2 


A(2)  + Ad) 


4 


The  options  for  OUTPUT  TINS ’the  STIFFNESS  MATRIX  are 
d)  Upper  triangular  portion  in  algebraic  format 
(2)  FORTRAN  CARD  IMAGE  format 
ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 
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2 

(K  - X ) 

2 1 


ROU  1,  COL  1 

- (XniN  - XtlAX)  (A(2)  XniN  - A(l)  XniN  + A(2)  XMAX  - Ad)  XHAX 

+ 2 Ad)  X - 2 X AI2))/(2  (X  - X )) 
' 2 1 2 1 


ROW  2.  COL  1 

(XfllN  - Xt1AX)  (A(2)  XniN  - Ad)  XI11N  + A(2)  KIIAX  - Ad)  XHAX 

+ 2 Ad)  X - 2 X A(2)  )/(2  (X  - X ) ) 

2 1 2 1 


ROU  2,  COL  2 

- (XI1IN  - XMAX)  (A(2)  XHlN  - Ad)  XI1IN  + A(2)  XI1AX  - Ad)  XMAX 

t 

+ 2 Ad)  X - 2 X A(2))/(2  (X  - X )) 
2 1 2 1 


T)ie  options  for  OUTPUTTING  tlie  STIFFNESS  MATRIX  are 
d)  Upper  triangular  portion  in  algebraic  format 
(2)  FORTRAN  CARD  IMAGE  format 
ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 

2; 

Enter  the  FORTRAN  ARRAY  NAME  for  the  STIFFNESS  MATRIX 
? 

STIFF; 

Multiply  each  coefficient  by 
E/(X(2)-Xd))Vf*2 
The  STIFFNESS  MATRIX  Is 

ST  I FF  d , 1 ) - - (XMI N-XMAX)  * f A (2)  *XMI  N-A  (1 ) *XMI  N+A  (2)  v.XMAX-A  ( 1 ) ivXMAX 
1 +2*A  d ) >vX  (2)  -2*X  d ) >vA  (2) ) / (2*  (X  (2)  -X  ( 1 ) ) ) 

STI  FF  (2. 1 ) - (XMI  N-XMAX)  * (A  (2)  ivXMI  N-A  d ) -.vXM I N+A  (2)  >vXMAX-A  ( 1 ) >vXMAX  t- 
1 2*A  (1 ) >vX  12)  -2*X  (1 ) *A  (2) ) / (2*  (X  (2)  -X  (1 ) ) ) 

STIFF  (3, 1 ) - - (XMI N-XMAX) ,v(A (2)  *XMIN-A  (1 ) *XM) N4A  (2)  *XMAX-A  ( 1 ) »vXMAX 


1I.I.USTRAT1  V(:  EXAI1I-1ES 


EiU.;ir  i.h<;  FDRIRAN  ARRAY  NAIIE  for  the  STIFI-NESS  MATRIX 

7 

STIFFS 


flultiply  n^nch  coefficient  by 


2v£/(X(2.l-X(l)) 


The  STIFFNESS  MATRIX  is 

STIFFd.l)  - {A(2)+A(]))/4 
STIFF(2,1)  = -(A(2)-*A(l))/4 
STIFF(3,1)  = (A(2)+A(l))/4 


3.2.3  Contpor  i son  aneJ  Verification 

< The  files,  (%2G,valG)  and  (%2I,valG)  cantairi  values  generated  during 
the  runs  in  sections  3.2.1  and  3.2.2  respectively.  As  in  the 
verification  of  section  3.1.3,  the  niinimuin  and  maxirnuin  values  of  X are 
at  nodes  1 and  2 respectively.  > 

timR=  8 msec. 

(Dl)  mSK,  AKIG] 

(C2)  L0A(]riLEC%2I,VALS); 

%2I  VAL5  D5K  AKIG  being  loaded 
loading  done 
tiinR»  308  msec. 

(LT2)  DONE 

(C3)  FTL:FUNCTIONAL-.vFACTOR; 
time-  14  msec. 

I (A(2)  + Ad))  E ) 

[ ) 

[ 2 (X  - X ) ) 

[ 2 1) 

[ ] 

t (A(2)  + Ad))  E ) 

(03)  f ) 

( 2 (X  - X ) ) 

( 2 13 

[ ] 

[ (A(2)  + Ad))  E I 

[ — — ] 

[ 2 (X  - X ) ) 

t 2 13 


(C4)  L0AnriLEdl2G.VAL5); 


ILl.USTRAllVH  L'XAtIPLtS 


y.?G  VAl.G  Cl'iK  AtflG  bnincj  loaded 
loadtoq  (lone 
t i ino"  377  mncc. 

(LVi)  / DONE 


(C5)  F TL2 :rUNCT I DUAL *f ACTOR; 
t i mo"  18  ntfsec. 

(DO)  IIATRlXd-  E (XniN  - XI1AX)  IA(2)  XHIN  - All)  XIUN  + A(2)  XflAX 


Ad)  XflAX  + 2 Ad)  X - 2 X A(2))/(2  (X  - X ) )], 

2 1 2 1 


(E  (X(1)N  - XflAX)  (A(2)  XnJN  - Ad)  XfllN  + A(2)  XI1AX  - 

3 

■+  2 Ad)  X - 2 X 

A(2))/(2  (X  - X ) )). 

2 1 

2 1 

1-  E (XI11N  - XMAX) 

(A  (2)  XI11N  - Ad)  XdlN  + 

A (2)  XdAX 

3 

+ 2 Ad)  X - 2 X 

A(2))/(2  (X  - X ) ))) 

2 1 

2 1 

(CG)  FACTOR  (EV  (%.  Xni N=X  Cll , XflAX-X  [21 ) ) ; 

time=  711  msec. 

[ (A(2)  + Ad))  E 

] 

[ 

] 

[ 2 (X  - X ) 

1 

[ 2 1 

] 

C 

) 

[ (A(2)  + Ad))  F. 

] 

(GG) 

( 

] 

[ 2 (X  - X ) 

) 

[ 2 1 

1 

[ 

) 

[ (A(2)  + Ad))  E 

) 

[ 

] 

[ 2 (X  - X ) 

1 

( 2 1 

] 

< Note  that  expressions  (DG)  and  107)  are  identical.  > 
(C7)  CL0SEFILE(WER,C2); 
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ii.LU‘,TRATivi:  iixAiirurs 


3.3  C5T  with  Uniform  Tliicknnt.s 

< Thq/constant  strain  triangle  (CST)  has  two  rlecjrecps  of  freedom  at 
^ each  node  and  is  modeled  with  a linear  displacement  function.  > 

3.3.1  Generalized  Coordinate  Formulation 


* SYSTEM  INITIALIZATION  * 

UELCOME  TO  *****  VERSION  1.0 

It  is  now  THURSDAY  OECEMSER  1,1977  2:35:9 
1 he  current  file  is  ISYINIT,  FCN] 

The  current  device  and  username  is  IU5K,  ALIGI 

Report  problems  to 
Al.AN  R.  KORNCOFF 
DEPT.  OF  CIVIL  ENGINEERING 
CARNEG I E -MELLON  UNI  VERS I TY 
CMLI-10A,  AKIG 


Terminate  all  input  with  a SEMICOLON  - 


Input  *?;’  for  HELP 


Input  your  LOGIN  NAME 
? 

AiaG; 

Is  AKIG  correct  ? (YES;  or  NO;) 

? 

YES; 

GREETINGS  AKIG 
* METHOD  SELECTION  * 

The  available  formulations  include 

(1)  THE  ISOPARAMETRIC  METHOD 

(2)  THE  GENERALIZED  COORDINATE  METHOD 

Please  enter  the  number  of  the  method  chosen  (1  OR  2). 
? 

2; 

GENERALIZED  COORDINATE  FORMULATION 


...  rcMroAi  i7cn  rnnnniMATr  cnonin  ATinw  pyrn iti vr  * 


Il.l.USlRATIVi;  EXAUn.fTS 
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NO; 


* PROOLEfI  PARAMETER  SPECIFICATION  - GENERALIZED  COORDINATE  * 

Inpul  tho  NOMDER  OF  ELEMENT  NODES 
? 

3; 

Input  lh(j  NI.IIIDER  OF  OEGREES  OF  FREEDOM  PER  NODE 
? 

: : . 2; 

• 

Inpvjt  the  vector  ot  the  NAMES  OF  THE  GLOBAL  COORDINATES 
There  should  be  2 elements 
ELEMENT  1 - 
X;  ■ 

ELEMENT  2 - 
Y; 


Input  the  vector  of  the  NAMES  OF  THE  DISPLACEMENT  VARIABLES 
There  r.liould  be  2 elements 
ELEMENT  1 » 

U; 

ELEMENT  2 - 
V; 


* SIIAPE  FUNCTION  PROCESSOR  - GENERALIZED  COORDINATE  * 

ENTER  the  terms  of  the  SHAPE  FUNCTION  ordered  from 
GENERALIZED  COORDINATE  1 through  coordinate  3. 

ELEMENT  1 - 
1; 

ELEMENT  2 - 
X; 

ELEMENT  3 - 
Y; 

SHAPE  FUNCTION  MODIFICATION 
The  OPTIONS  are 

(1)  DISPLAY  THE  SHAPE  FUNCTIONS 

(2)  MODIFY  THE  SHAPE  FUNCTIONS 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

3; 

,v  0 MATRIX  DATABASE  GENERATION  * 


The  OPTIONS  for  specifying  STRAIN  COMPONENTS  arc 

(1)  USER-SUPPLIED  VALUES 

(2)  LIBRARY  VALUES 


A 


A 


II  l.U:-.TRATIV(;  F.XAIIPLr-S 

P Iodise  tNIER  thR  NUlWER  ASSOClAlbD  UlTH  YOUR  5IILECTIOI4 
? 

2; 

The  LIBRARY  OPTIONS  for  SPECIFYING  STRAIN  COMPONENTS  are 

(1)  USER-SUPPLIED  VALUES 

(2)  ONE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 

(4)  PLANE  STRAIN 

(5)  AKISYMMETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 30 
Please  ENTER  the  NUMBER  ASSOCIATED  UITH  YOUR  CHOICE 
? 


•.V  MATERIAL  PROPERTIES  SELECTION  * 

The  ojjtions  for  the  selection  of  the  MATERIAL  PROPERTIES  MATRIX  are 

(1)  USER-SUPPLIED  MATRIX 

(2)  LIBRARY  MATRIX 

Please  enter  the  NUMBER  ASSOCIATED  WITH  YOUR  SELECTION 
? 


The  LIBRARY  OPTIONS  for  SPECIFYING  MATERIAL  PROPERTIES  are 

(1)  USER-SUPPLIED  VALUES 

(2)  ONE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 
(41  PLANE  STRAIN 
(5)  AX  I SYMMETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 30 
Please  ENTER  the  NUMBER  AS50CIATE0  UITH  YOUR  CHOICE 
? 

3; 

* AUXILIARY  TERM  PROCESSOR  * 


ENTER  ELEMENT  VOLUME. MODIFICATION  FACTORS  OR  AUXILIARY  TERMS 
TYPE  ’END;*  to  TERMINATE 
? 

T; 


? 


? 


< As  indicated  in  section  2.2.7,  the  nominal  element  volume  is 


calculated  as  the  product  of  the  differentials  of  the  coordinate 
variables  (in  thin  cane  as  dxivdy) . Since  the  area  of  the  triangle  is 
half  this  quantity,  a volume  modification  term,  equal  to  .5  is 
specified.  The  term,  "T",  represents  the  uniform  thickness  of  the 
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AUXILIARY  TERM  MODI  El  CAT  ION 
The  OPTIONS  are  . 

(1)  DISPLAY  THE  AUXILIARY  TERMS 

(2)  /MOniFY  AN  AUXILIARY  TERM 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 


* INTEGRATION  LIMITS  - GENERALIZED  COORDINATE  * 

* INVERSION  OF  MATRIX.  A »v 

,v  D MATRIX  GENERATION  - GENERALIZED  COORDINATE  * 

,v  INTEGRATION  PROCESSOR  - GENERALIZED  COORDINATE  * 
Vf  DISPLAY  PRE-PROCESSOR  - GENERALIZED  COORDINATE  * 

* DISPLAY  PROCESSOR  * 

The  options  for  OUTPUTTING  the  STIFFNESS  MATRIX  are 

(1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARO  IMAGE  format 

ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

? 

i; 


Multiply  each  coefficient  by 

-ET/(2(X  Y - X Y - Y X + Y X + X Y - Y 
2 3 1 3 2 3 1 3 1 2 1 


2 

X ) 

2 


(NU  - 1) 

(NU  + D) 


RON  1,  COL  1 


2 

2 

2 

2 

2 

- (X 

NU  - 2 X X NU  + X NU 

- 2 Y +4 

Y Y 

- X + 2 X X 

- 2 Y 

3 

2 3 2 

3 

2 

- X ) 
2 

2 3 

(XMIN 

3 2 3 

- XMAX)  (YMJN  - 

2 

YMAX)/2 

ROW 

2,  COL  1 

2 

2 

(X 

NU  - X X NU  - X X NU  + X X NU 

-2Y  + 2Y  Y + 2Y  Y 

3 

2 3 13  12 

3 2 3 1 : 

r - T 

t ILl.USTRATIVK  L-XAI1PLE3 


2 

-X+XX+XX-2YY-XX)  (XHIN  - XIIAX)  (YI11N  - Yf1AX)/2 
3/23  13  12  12 


ROU  2.  COL  2 

2 2 2 2 2 
-(X  NU-2X  X NU  + X NU-2Y  + AY  Y - X + 2X  X - 2Y 
3 13  1 3 13  3 13  1 

2 

- X ) (XI1IN  - XMAX)  (YfllN  - YriAX)/2 
1 


RDU  3.  COL  1 

2 

-(XX  NIJ  - X X NU-X  NU  + X X Nu-2Y  Y + 2Y  Y - X X 
23  13  2 12  23  1323 

2 2 

+ X X + 2Y  - 2-Y  Y + X - X X)  (XI11N  - XNAX)  (YH1N  - YMAX)/2 
1 3 2 1 2 2 1 2 


ROW  3,  COL  2 

2 

(X  X NU-X  X NU-X  X NU  + X NU-2Y  Y +2Y  Y - X X 
,2  3 13  12  1 23  1323 

2 2 

+ X X +2Y  Y + X X - 2Y  - X)  (XI1IN  - XI1AX)  (YniN  - YnAX)/2 
13  12  12  11 


ROW  3,  COL  3 

2 2 2 2 2 
-(X  NU-2X  X NU  + X NU-2Y  + 4Y  Y - X +2X  X - 2Y 
2 12  1 2 122  121 

2 

- X ) (XI11N  - XMAX)  (YtllN  - YnAX)/2 
1 


ROW  4,  COL  1 


(X  - X ) (Y  - Y ) (NU  + 1)  (XniN  - XMAX)  (YMIN  - YMAX) 
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3 2 3 2 


2 

/' 


ROW  4, 

COL  2 

(X  Y 

NU  - 2 X Y NU  + X Y NU  h 

■ Y X NU  - 

2 Y X NU  - X Y NU 

3 3 

2 3 

1 3 

2 3 

13  12 

+ 2 Y 

X NU  + X Y 

- X Y - Y 

X 4-  X Y ) 

(XI1IN  - XIIAX) 

12  3 3 

1 3 2 

3 1 2 

(YI1IN 

- YnAX)/2 

ROW  4, 

COL  3 

(X  Y 

NU  - X Y NU  - 

- 2 Y X NU  4 

■ 2 Y X NU 

+ X Y NU  4 X Y NU 

2 3 

1 3 

2 3 

1 3 

2 2 12 

- 2 Y 

X NU  - X Y 

4-  X Y 4-  X 

Y - X Y ) 

(XniN  - XI-1AX) 

12  2 3 

1 3 2 

2 1 2 

(YI1IN 

- YriAX)/2 

ROW  4, 

COL  4 

2 

2 2 

2 2 

- (Y  NU  - 2 Y Y NU 

4-  Y NU  - Y 

4-  2 Y Y - 

2X  4-4X  X - Y 

3 

2 3 

2 3 

2 3 

3 2 3 2 

2 

- 

2 X ) (xniN 

- XNAX)  lYMlN  - YI1AX)/2 

2 

ROW  5. 

COL  1 

(X  Y 

NU  + X Y NU  - 

• 2 X Y NU  - 

2 Y X NU 

4-Y  X NU4-2X  Y NU 

3 3 

2 3 

1 3 

2 3 

13  12 

- Y X NU  + X Y - 

X Y - Y X 

4-  Y X ) (X(1)N  - XflAX) 

1 

2 3 3 

2 3 13 

1 2 

(YHIN 

- Yf1AX)/2 

ROW  5, 

COL  2 

(X  - X ) (Y 

- Y ) (NU  4-  1)  (xniN  - xnAX)  (vniN  - yhax) 

3 1 3 

' 1 

ILLUSTRATIVi:  EXAUPLES 


2 


ROW  5.  COL  3 

- (X  Y NU  - X Y NU  - 2 Y X NIJ  + 2 Y X NU  + 2 X Y NU 

2 3 1 3 -2  3 1 3 1 2 

-Y  X NU-X  Y NU-X  Y + X Y +Y  X - X Y)  (XHIN  - XMAX) 

12  11  23  13  12  11 

i 

(YU IN  - YnAX)/2 


ROW  5,  COL  4 

2 2 2 
(Y  NU  - Y Y NU  - Y Y NU  + Y Y NU  - Y + V Y + V V - 2 X 
3 23  13  12  32313  3 

+ 2X  X +2X  X - Y Y - 2X  X)  (XIUN  - XMAX)  (YdIN  - YriAX)/2 
2 3 13  12  12 


ROW  5,  COL  5 

2 .2  2 2 2 
-(Y  NU-2Y  Y NLI  + Y NU-Y  + 2Y  Y - 2X  + 4X  X - Y 
3 13  13  13  3 131 

2 

- 2 X ) (XI11N  - XdAX)  (YdIN  - YdAX)/2 
1 


ROW  G.  COL  1 


(2 

X 

Y 

NU  - 

2 X Y NU 

- Y X NU  + Y 

X 

NU-X  Y NU 

2 

. 3 

1 3 

2 3 1 

3 

2 2 

+ 2 

X 

Y 

NU  - 

Y X NU  + 

Y X - Y X - 

X 

Y + Y X ) 

1 

2 

1 2 

2 3 13 

2 

2 1 2 

(XniN  - XdAX)  (YdIN  - YdAX)/2 


ROW  G, 

, COL  2 

(2  X 

Y NU  - 

2 X 

Y’  NU  - 

Y X 

NU  + Y 

X NU  + X 

Y NU  - 2 Y X NU 

2 

3 

1 

3 

2 3 

1 

3 1 

2 12 

+ X 

Y NU  + 

Y X 

- Y X 

- X 

Y + X 

Y ) (XdIN 

- XdAX) 

ILLlll.TKAT  I Vt  EKAIIPLtS 
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2 3 13  12  11 


(YI11N  - YnAX)/2 


ROW  G.  COL  3 


(X  - X ) (Y  - Y •)  (NU  + 1)  (XtllN  - XnAX)  (YfllN  - YHAX) 
2 12  1 


ROU  G.  COL  4 


-(Y  Y NU-Y  Y NU-Y  NU  + Y Y NU-Y  Y +Y  Y - 2X  X 
23  13  2 12  2313  23 

2 2 

+ 2X  X + Y - Y Y + 2X  - 2X  X)  (XI11N  - XflAX)  (YIUN  - Yf1AX) 
1 3 2 1 2 2 1 2 


ROU  G,  COL  5 


(Y  Y NU-Y  Y NU-Y  Y NU  + Y NU-Y  Y + Y Y - 2 X X 
23  13  12  1 2313  23 

2 2 

+ 2X  X + Y Y +2X  X - Y - 2X)  (XUIN  - XI1AX)  (YIIIN  - YHAX) 
13  12  12  1 1 


ROU  G,  COL  G 

2 2 2 2 2 
-(Y  NU-2Y  Y NU  + Y NU-Y  +2Y  Y - 2X  + 4X  X - Y 
2 12  1 2 12  2 12  1 


2 X ) (XI11N  - XIIAX)  (YHIN  - Yt1AX)/2 
1 


< Display  option  2 produces;  > 


ILLUlillfATIVi:  LXAI1PLC5 
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The  options  for  OUTPUTTING  the  STIPFNL'SS  I1ATITIX  ore 

(1)  Upper  tricinyular  portion  in  algebraic  format 

(2)  'FORTRAN  CARO  II1AGE  format 

ENTER  the  NUtlBER  ASSOCIATED  with  YOUR  SELECTION. 

? 

Enter  the  FORTRAN  ARRAY  NAHE  for  the  STIFFNESS  MATRIX 
STIFF; 


f1(jlti|jly  each  coefficient  by 

-E*T/  (2*  (X  (2)  *Y  (3)  -X  (1 ) *Y  (3)  -Y  (2)  -.vX  (3)  i Y (1 ) >vX  (3)  +X  (1 ) ivY  (2)  -Y  (1 ) *X  ( 
■1  2))**2-,v(NU-l)*(NU+l)) 

The  STIFFNESS  MATRIX  is 

STI FF (1 . 1 ) =.  - (X (3) •.v.v2*MU-2*X (2) *X (3) .vNU+X (2) **2>vNU-2-.vY (3) *-,v2+A*Y ( 

1 2)  *Y  (3)  -X  (3)  **2-( 2*X  (2)  *X  (3)  -2-,vY  (2)  ,v*2-X  (2)  >Wc2)  ,v (XMl N-XMAX)  »v  ( YMI 

2 N.-YMAX)/2 

STIFF  (2,1)  = (X  (3)  >-nv2*NU-X (2) -.vX (3) -.vNU-X  (1 ) .vX  (3)  *MU^X  (1 ) *X  (2)  ■.vNU-2Vf 

1 Y (3)  **24-2*Y  (2)  :vY  (3)  +2*Y  (1)  -.vY  (3)  -X  (3)  **2+X  (2)  *X  (3)  +X  (1 ) ,vX  (3)  -2*Y 

2 (1 ) *Y  (2)  -X  (1 ) VfX  (2) ) * (XMl N-XMAX)  -.v ( YMI N- YMAX)  /2 

ST  I FF  (3 , 1 ) - - (X  (3 1 **2*NLI-2*X  ( 1 ) ivX  (3)  ivNlJ+X  ( 1 ) !V,v2>vNU-2-.vY  (3)  >v,v2+A*Y  ( 

1 1 ) *Y  ( 3 ) ••  X ( 3 ) **2  (-2*X  ( 1 ) ,vX  ( 3 ) -2*Y  ( 1 ) ,v»2- X ( 1 ) **2 ) * ( XM 1 N-Xf1AX ) * ( YM I 

2 N-YMAX)/2 

ST ! FF  (A , 1 ) = ■-  (X  (2)  *X  (3)  ;VM1J-X  ( 1 ) *X  (3)  .vNU-X  (2)  •.v*2*NU+X  ( 1. ) ,vX  (2)  *NU-2 

1 * Y ( 2 ) *Y  ( 3 ) 1 2*Y  ( 1 ) -.v Y ( 3 ) -X  ( 2 ) *X  ( 3 ) (-X  ( 1 ) ,vX  ( 3 ) 1 2*Y  ( 2 ) **2-2*Y  ( 1 ) >vY  ( 2 

2 ) -t-X  (2)  **2-X  ( 1 ) >vX  (2) ) .V  (XMI N-XMAX)  * (YMI  N- YMAX)  /2 

STIFF  (5, 1 ) = (X  (2)  -.vX  (3)  >vNU-X  (1 ) ,vX  (3)  ivNU-X  (1 ) ,vX  (2)  *NU+X  (1 ) ,v;.-2*NU-2* 

1 Y (2)  *Y  (3)  +2vfY  (1)  *Y  (3)  -X  (2)  *X  (3)  -i-X  (1 ) *X  (3)  +2*Y  (1)  *Y  (2)  +X  (1 ) *X  (2) 

2 -2*Y  (1 ) Vov2-X  (1 ) **2)  * (XMI  N-XMAX)  * (YMI  N-YMAX)  /2 

STIFF(G,1)  = -(X(2)  **2*NIJ-2*X  ( 1 ) *X  (2 ) ,vNU+X  ( 1 ) Vnv2*NU-2*Y  (2 ) iV*2+A  *Y  ( 

1 1 ) .vY  (2)  -X  (2)  ,v*2+2*X  (1)  ,vX  (2)  -2VfY  (1)  ,-nv2-X  (1 ) **2)  Vr  (XMIN-XI1AX)  * (YMI 

2 N-YMAX) /2 

STIFF (7,1)  - -(X(3)-X(2))*(Y(3)~Y(2))>v(NLNl)>v(XI1IN-XMAX)*(YniN~YnA 
1 X)/2 

S T I FF  ( S , 1 ) - ( X ( 3 ) ,v Y ( 3 ) >-fNU-  2*X  ( 2 ) * Y ( 3 ) >vNU+X  ( 1 ) * Y ( 3 ) *NU-v Y ( 2 ) *X  ( 3 ) >vN 

1 U-2*Y  ( 1 ) ,vX  ( 3 ) -.vNU-X  ( 1 ) *Y  ( 2 ) •.vNU+2*Y  ( 1 ) >vX  ( 2 ) -.vNU+X  ( 3 ) *Y  ( 3 ) -X  ( 1 ) >vY  ( 3 

2 ) -Y  (2)  *X  (3)  +X  (1 ) *Y  (2) ) * (XMI  N-XMAX)  -.v  (YMI  N-YMAX)  /2 

STIFF  0,1)  - (X  (2)  >vY  (3)  i-^NU-X  (1)  ,vY  (3)  ,vNU-2*Y  (2)  ,vX  (3)  *NU+2*Y  (1 ) *X  (3) 

1 *NU+X  (2)  >vY  (2)  *NU+X  (1 ) *Y  (2)  *NU-2vfY  (1 ) ,vX  (2)  *NU-X  (2)  *Y  (3)  ^ X ( 1 ) vrY  (3 

2 ) ■»  X ( 2 ) *Y  ( 2 ) - X ( 1 ) *Y  ( 2 ) ) Vf  ( XM  1 N -XM  AX ) * ( YM  1 N - YM  AX ) /2 

ST  I FF  (1 0. 1 ) " - (Y  (3)  iV>v2*NU-2*Y  (2)  >vY  (3)  -.vNU+Y  (2)  »v,v2*NU-Y  (3)  ,v*24.2*Y  (2 

1 ) *Y  (3)  - ^^X  (3)  >-nV?(  AvcX  (2)  *X  (3)  ••  Y (2)  **2-2>vX  (2)  Vnv2)  »v  (XMI  N-XMAX)  * (YM 

2 1 N-YMAX) /2 

ST  I FF  (1 1 , 1 ) - (X  (3)  *Y  (3)  ,vNlJ4X  (2)  *Y  (3)  .•.Nll-2-.vX  ( 1 ) >vY  (3)  v.NU-2*Y  (2)  -.vX  (3 

1 ) v.NU+Y  (1 ) *X  (3)  *NU42*X  (1 ) >vY  (2)  >-.NlJ-Y  (1 ) ,vX  (2)  ^vNl.UX  (3)  *Y  (3)  -X  (2)  *Y  ( 

2 3 ) - Y ( 1 ) .vX  ( 3 ) 4 Y ( 1 ) v.X  ( 2 ) ) iv  ( XM  I f J - XMAX ) iv  ( YM  1 N-YMAX )/  2 

STIFF  (12,1)  - -(X(3I-X(11)*(Y(3I-Y(1)  ),v(NU+l).v(XniN-X'MAX)*(YMIN-YM 
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1 AX)/2 

5T 1 re  ( J 3. 1 ) =>  - (X  (2)  »vY  (3)  »vNlJ-X  ;i ) >vY  (3)  ivMU-2-,vY  (2)  >vX  (3)  >'.Nll  (2-.vY  ( 1 ) ,vX  ( 

1 3)  ,vM|l4  2vvX  ( J. ) ivY  (2)  )VNIl.l-y  (1 ) -.vX  (2)  *N1I-X  (i ) ,vY  (1 ) *MU-X  (2)  -.vY  (3)  -tX  ( J ) -.vY 

Z ( c •(  Y ( ]. ) >vX  ( 2 ) - X ( ] ) ,v Y ( 1 ) ) Vf  ( XI1  i N - XMAX ) ,v  I Ytl ) N-  N'I1  AX ) /2 
STIFF  (J  A ,!)•■•  (Y (3)  >v*2*NLI-Y (2)  *Y (3) >vNU-Y (1 ) ,vY (3)  -.vNll+Y  (i ) ,vY (2)  ivNU-Y 

1 (3)  *:v2-t  Y (2)  *Y  (3)  -t  Y (1 ) ,vY  (3)  ••2-.vX  (3)  >v;;2  (2;vX  (2)  -.vX  (3)  +2*X  ( 1. ) *X  (3)  -Y  ( 

2 1 ) *Y  (2)  -2*K  (1 ) ivX  (2) ) * (XIII N-XI1AX)  -.v  ( YIH  N-YI1AX)  /2 

ST  I FF  (IS.  1 ) - - (Y  (3)  *,v2*NU-2*Y  (1 ) kY  (3)  .vNU+Y  (1 ) .v*2*l>)U-Y  (3)  >vv;2-(-2*Y  ( 1 

1 ) -.vY  (3 ) -2ivX  ( 3 ) iV*2-i  A *X  ( 1 ) VrX  ( 3 ) - Y ( 1 ) ,v*2-2-.vX  ( I ) **2)  * ( X(1 1 N-XI1AX ) ,v  ( YH 

2 II\l-YI1AX)/2 

ST  I FF  (IG,  1 ) - - (2*X  (2)  *Y  (3)  ;vNI.I-2»X  (1 ) -.vY  (3)  -.vNU-V  (2)  ;vX  (3)  V.MLI+Y  (1 ) *X  ( 

1 2)  .vMlJ-X  (2) *Y  (2)  •,vNl.l+2*X  ( ]. ) ivY  (2)  >vNI I-Y  (1 ) -.vX  (2)  ivNU+V  (2)  >vX  (3)  -Y  (1 ) ,vX 

2 (3)  -X  (2) -.vY  (2)  -f-Y  (1 ) ,vX  (2) ) ,v(X(1lN-XnAX)  ,v ( Ylll  N-Yt1AX)  /2 

STI FF  (17, 1 ) = (2*X  (2)  ,vY  (3)  ,vNlJ-2*X  (1 ) *Y  (3)  -.vNlJ-Y  (2)  -.vX  (3)  •.vNlJ+Y  (1 ) ,vX  (3 

1 ) VfNU-fX  (1 ) ,vY  (2)  ,vNU-2,vY  (1 ) *X  (2)  ;vNll4X  ( 1 ) *Y  (1 ) ivNl  I t Y (2)  *X  (3)  -Yd)  -.vX  ( 

2 3 ) -X  ( 1 ) * Y ( 2 ) -i-X  ( 1 ) ivY  ( 1 ) ) ,v  ( X(1 ) N -X(1  AX ) * ( YI1 1 N - YI1  AX ) / 2 

STIFF  (18,1)  = -(X(2)-X(l)),v(Y(2)-Y(l))>v(NlJ4l)>v(XtHN-Xf1AX)-.v(VI1IN-Yf1 
1 AX)/2 

STIFF  (13, 1 ) - - (Y  (2)  >vY(3)  •.vMlJ-Y  (1 ) >vY  (3)  -.vNU- Y (2)  •.v>v2*Nl.l4-Y  (1 ) -.vY  (2)  -.vIMU- 

1 Y ( 2 ) Y ( 3 ) 4 Y ( 1 ) ,v Y ( 3 ) -2-,vX  ( 2 ) -.vX  ( 3 ) 4 2-.vX  ( J ) >vX  ( 3 ) 4 Y ( 2 ) .v:v2 - Y ( 1 ) ,v Y ( 2 ) 4 2 

2 *X  (2)  **2-2*X  ( 1 ) ,vX  (2) ) ,v  (Xtll N-XI1AX)  * ( YI1I N- YI1AX)  /2 

ST  I FF  (20. 1 ) ..  ( Y (2)  ,vY  (3)  .vNlJ-Y  (1 ) ivY  (3)  >vNl  I- V ( J ) -.vY  (2)  »vNl .'4  Yd)  ,v,v2vrNU-Y 

1 (2)  ivY  (3)4  Yd)  ,vY  (3)  ~2*X (2)  >vK  (3)  42*X  d ) ,vX  (3)  4 Y d ) ,vY  (2)  +2*X  (1 ) *X  (2 

2 ) -Y  (1 ) ,-nv2-2*X  (1 ) ,v,v2)  *{Xni  N-XNAX)  -.v  ( YI1I N-Y11AX)  /2 

ST  I FF  (21 . 1 ) = - ( Y (2)  **2-.vNU-2*Y  (1 ) ,vY  (2)  =\NLl4-Y  (1 ) *v.-2-.vf>IU-y  (2)  iV*2+2*Y  d 

1 ) *Y  (2)  -2^'rX  (2)  >v*2+Ai’fX  ( 1 ) jvX  (2)  - Y (1 ) ii'nv2-  2*X  ( 1 ) **2 ) * ( X(1 1 N-X11AX ) * ( Vn 

2 IN-YI1AX)/2 


3.3.2  Isoparametric  Formulation 


* SYSTEM  INITIALIZATION  * 

UELCOME  TO  ***;w.-  VERSION  1.0 

It  is  noi4  THURSDAY  DFCEMBER  1,1977  2:18:4 
The  current  file  is  tSYlNIT,  FCN) 

The  current  device  and  username  is  IDSK,  AXIG) 

Report  problems  to 
ALAN  R.  KORNCOFF 
DEPT.  OF  CIVIL  ENGINEERING 
CARNEG 1 E-MELLON  UNI  VERS  i T Y 
CMU-10A,  AKIG 

Terminate  all  input  with  a SEMICOLON  - 
Input  •?; • for  HELP 


Input  your  LOGIN  NAME 
? 
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ILLUl^TRATIVL  CXAIIPLES 
AKIG; 

!s  AKIG  correct  ? (YES;  or  NO;) 

? 

YES;.  / 

GREETINGS  AKIG 
>v  NETHOO  SELECTION  * 

The  available  formulations  include 

(1)  THE  1 SOP ARAtlE TRIG  HETHOD 

(2)  THE  GENERALIZED  COORDINATE  HETHOD 

Please  enter  the  number  of  the  method  chosen  (1  OR  2). 

? 

1 ; 

I SOPARAHE  TR I C FORHULAT I ON 

* 1 SOPARAHE TR I C FORHULAT I ON  EXECUTIVE  * 

* ISOPARARHETRIC  FORHULATION  INITIALIZATION  * 

Do  you  wish  to  set  DUHP,  BREAK  or  TRACE  POINTS 
? 

NO; 

* PROBLEH  PARAHETER  SPECIFICATION  - ISOPARAHETRIC  * 

Input  the  NUHBER  OF  ELEHENT  NODES 
? 

3; 

Input  the  NUHBER  OF  DEGREES  OF  FREEDOH  PER  NODE 
? 

7. 

Input  the  NUHBER  OF  NATURAL  COORDINATES 
? 

3; 

Input  the  vector  of  the  NAMES  OF  THE  NATURAL  COORDINATES 
There  should  be  3 elements 
ELEHENT  1 - 
LI; 

ELEHENT  2 = 

L2; 

ELEHENT  3 - 
L3; 


Input  the  vector  of  the  NAMES  OF  THE  GLOBAL  COORDINATES 
There  should  bo  2 elements 
ELEHENT  1 » 
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ILl.UL.TIMTIVC  EXAIH-LL'S 
X; 

ELtflliNT  2 » 

Y; 

Input  the*  vector  of  the  NAMES  OF  THE  OISPLACEflENT  VARIABLES 
There  should  be  2 elements 
El.EMENT  1 .. 

U; 

ELEMENT  2 = 

V; 

* SHAPE  FUNCTION  PROCESSOR  - ISOPARAMETRIC  * 

ENTER  the  terms  of  the  SHAPE  FUNCTION  ordered  from 
node  1 through  node  3. 

The  3 elements  wi M.  be  prompted  for 
El-EMI-NT  1 » 

LI; 

ELEMENT  2 = 

L2; 

ELEMENT  3 - 
L3; 

SHAPE  FUNCTION  MODIFICATION 
The  OPTIONS  iire 

(.1)  DISPLAY  THE  SHAPE  FUNCTIONS 

(2)  MODIFY  THE  SHAPE  FUNCTIONS 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 

i ! 

* B MATRIX  DATABASE  GENERATION  * ; 

I 

The  OPTIONS  for  specifying  STRAIN  COMI’ONENTS  ore  I 

(1)  USER-SUPPLIED  VALUES 

(2)  LIBRARY  VALUES 

PIceJse  ENTER  the  NUMBER  ASSOCIATED  WITH  YOUR  SELECTION  ! 

? 

l5 

< Though  the  following  strain  component  specifications  may  be 
automatically  selected  by  chosing  library  option  3 or  A,  the  fol lowing 
demonstrates  the  user  input  option.  > 

I 

Input  the  NUMBER  OF  STRAIN  COMPONENTS 


Input  COMPONENT  1 

? I 

D(U.X);  j 

Input  COMPONENT  2 1, 

1 


J 
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ILl.USTRA  f IVC  IIXAI1l-l.r5 
? 

D(V.Y) ; 

Input  COtIPONIiNT  3 
? 

D(U,Y)  + D(V,X); 


>v  MATERIAL  PROPERTIES  SELECTION  * 

Thtt  options  for  the  selection  of  the  MATERIAL  PROPERTIES  MATRIX  are 

(1)  USER-SUPPLIED  MATRIX 

(2)  LIBRARY  MATRIX 

Please  enter  the  NUMBER  ASSOCIATED  WITH  YOUR  SELECTION 
? 

1;  . 

< T-hout.)h  the  follouincj  material  property  specifications  may  be 
automatically  selected  by  chosiny  library  option  3,  the  following 
demonstrates  the  user  input  option.  > 

Please  enter  the  UPPER  TRIANGULAR  PORTION  of 
the  MATERIAL  PROPERTIES  MATRIX. 

Eacli  element  will  be.  prompted  for  with  its  RON  and  COLUMN  number. 

You  can  also  specify  a CONSTANT  SCALAR  MULTIPLE. 

ERRONEOUS  INPUT  may  be  corrected  after  the  matrix  is  entered. 

Do  you  wish  to  enter  a CONSTANT  SCALAR  MULTIPLE  (YES;  OR  NO;) 

? 

YES; 

Enter  the  SCALAR  MULTIPLE 
? 

E/(,l-NU*Vf2); 

Enter  the  ELEMENTS  of  the  MATERIAL  PROPERTIES  MATRIX 
as  they  arc  prompted  for 
The  matrix  should  be  of  order  3 
RON  1 , COLUMN  1 - 
? 

1; 

ROW  1 , COLUMN  2 - 
? 

MU; 

RON  1 , COLUMN  3 = 

? 

0; 

RON  2 , COLUMN  2 - 
? 

1; 

RON  2 , COLUMN  3 - 
? 

0; 

RON  3 . COLUMN  3 - 
? 
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ll.I.lJf-TRATIVC  EXAI1l-l.ES 


l-Nl.l/2; 


I1ATER1ALS  hATRlX  HOLHFICATJON 
Thf!  OPTIONJ)  ciru 

(1)  DISPLAY  THE  MATERIALS  MATRIX 
MOniPY  THE  MATERIALS  MATRIX 
DISPLAY  THE  MATERIALS  MATRIX  MULTIPLIER 
MODIFY  THE  MATERIALS  MATRIX  MULTIPLIER 
TERMINATE  THIS  FUNCTION 

tin;  WIJMRER  ASSOCIATED  with  the  CHOSEN  OPERATION 


(2) 

(3) 

(4) 

(5) 


Enter 
? 

1; 

THE  MATRIX  IS 


[ 1 
[ 

[ NU 
[ 

[ 


NU 


] 

] 

I 

] 

NU  I 


[ Q 0 1 I 

[ 2 ] 


MULTIPLIED  BY 

E 


2 

1 - NU 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPERATION 
? 

2; 


MODIFY  THE  MATERI ALSMATRIX 

ONLY  THE  ELEMENTS  ABOVE  THE  DIAGONAL  NEED  BE  ENTERED 


ii.i.ur-.TnATiv(-;  t:xAfin.cs 
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* AUXILIARY  TFRU  PRnCFSSOR  * 

FfJ7nR'L‘LL-f1FNT  VDLUflF  nODIFICATJON  FACTORS  DR  AUXILIARY  TERI1S 
TYPE  ’END;’  to  TERMINATE 
? 

T; 

? 

END; 

AUXILIARY  TERM  MODIFICATION 

The  OPTIONS  are  * 

(J)  DISPLAY  THE  AUXILIARY  TERMS 

(2)  MODIFY  AN  AUXILIARY  TERM 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

3; 

INTEGRATION  LIMITS  - ISOPARAMETRIC  * 

Since  the  NATURAL  COORDINATES  Bpecifieci  are  NOT  INDEPENDENT, 
you  need  ENTER  only  the  LONER  LIMITS  OF  INTEGRATION 

NATURAL  COORDINATE.  LI 
UPPER  LIMIT  - 
1 

ENTER  the  LONER  LIMIT 
? 

0; 

NATURAL  COORDINATE.  L2 
UPPER  LIMIT  » 

1 - LI 

ENTER  the  LONER  LIMIT 
? 

1 ; 

THIS  LIMIT  OF  INTEGRATION  MUST  BE  AN  INTEGER.  EITHER  -1  OR  0. 

PLEASE  RE-ENTER. 

? 

0; 

< A sem.intic  check  on  the  input  flagged  an  error  in  the  specification 
of  the  lower  limit  of  integration.  > 

* JACOBIAN  GENERATION  - ISOPARAMETRIC  »v 

* B MATRIX  GENERATION  - ISOPARAMETRIC  * 

* INTEGRATION  PROCESSING  - ISOPARAMETRIC  * 

* DISPLAY  PRE-PROCESSOR  - ISOPARAMETRIC  * 


* DISPLAY  PROCESSOR  * 


ii.i.iJoTRAnvi-:  cxAnri.Ks 
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The  options  for  OUTF’UITING  the  STIF'FfCSS  MATRIX  ore 
(1)'  Upper  triangular  portion  in  algebraic  formal 
12)  F'ORTRAN  CARD  IMAGE  format 
ENTER  the  NUMBER  ASSOCIATED  with  YOUR  SELECTION. 

7 

1; 

Multiply  each  coefficient  bij 

E T 


(X  Y - X Y - y X + Y X + X Y - Y X)(NU-1)(NU  + 1) 
23  13  23  13  12  12 


RDU  1.  COL  1 

2 2 2 2 2 
-(X  NU-2X  X NLI  + X NU-2Y  +4Y  Y - X + 2X  X - 2Y 
3 23  2 3 233  23  2 

2 

- X )/4 
2 


ROIJ  2.  COL  1 

2 2 
(X  NU  - X X NU  - X X NU  + X X NU  - 2 Y + 2 Y Y + 2 Y Y 
3 23  13  12  3 23  13 

2 

-X  + X X +X  X - 2Y  Y - X X)/4 
3 23  13  12  12 


RDU  2.  COL  2 

2 2 2 2 2 
-IX  NU-2X  X NU  + X NU-2Y  +4Y  Y - X +2X  X - 2Y 
3 13  1 3 13  3 13  1 


- X )/4 
1 


ILLUS1RA7IVC-;  tXAIIPLtS 


- (X  X NU  - X X NU  - X NU  + X X NU  - 2 S Y + 1’  Y Y - X X 
23  13  2 12  23  1323 

, /'  2.  2 

+ X X + 2Y  - 2Y  Y + X - X X)/A 
1 3 2 1 2 2 1 2 


ROU  3,  COL  2 

2 

(X  X NU  - X X NU  - X X NU  + X NU  - 2 Y Y 2 Y Y - X X 
23  13  12  1 23  1323 

2 2 

+ X X + 2Y  Y + X X - 2Y  - X)/4 
13  12  12  1 1 


RON  3,  COL  3 

2 2 2 2 2 
-(X  NU-2X  X NU  + X NU-2Y  +4Y  Y - X +2X  X - 2Y 
2 12  1 2 12  2 12  1 

2 

- X )/4 
1 


ROW  4,  COL  1 

(X  - X ) (Y  - Y ) (NU  + 1) 
3 2 3 2 


4 


ROW  4,  COL  2 

(X  Y NU  - 2 X Y NU  + X Y NU  + Y X NU  - 2 Y X NU  - X Y NU 
33  23  13  23  13  12 

+ 2Y  X NU  + X Y -X  Y - Y X +X  Y)/4 
12  33  13  23  12 


RON  4.  COL  3 

(X  Y NU  - X Y NU  - 2 Y X NU  + 2 Y X NU  + X Y NU  + X Y NU 
23  13  23’  13  22  12 

-2Y  X NU-X  Y + X Y + X Y - X Y)/4 


ILUBTHAlIVt  LXAUPLUS 


12  2 2 12  2 2 1 

ROIJ  4./ COL  4 

2 2 2 2 
-(Y  NU-2Y  Y NU  + Y NU-Y  +2Y  Y - 2X  + 4X  X - 
3 23  2 3 23  3 23 

- 2 


ROU  5.  COL  1 

(X  Y NU  + X Y NIJ  - 2 X Y NU  - 2 Y X NIJ  + Y X MU  + 2 X 
33  23  13  23  13  1 

-Y  X NU  + X Y - X Y - Y X +Y 
1 2 3 3 2 3 1 3 1 


ROM  5,  COL  2 

(X  - X j (Y  - Y ) (NU  + 1) 
3 13  1 


4 


ROU  5.  COL  3 

- (X  Y NU  - X Y NU  - 2 Y X NU  + 2 Y X NU  + 2 X Y NU 
2 3 1 3 2 3 1 3 1 2 

-Y  X NU-X  Y NU-X  Y + X Y + Y X - X 
12  11  2 3 13  12  1 


ROU  5.  COL  4 

2 2 
(Y  NU-Y  Y NU-Y  Y NU  + Y Y NU-Y  +Y  Y + Y Y - 2 
3 23  13  12  32313 

+ 2X  X + 2X  X - Y Y - 2X 
2 3 13  12  1 


ROU  5,  COL  5 


2 


2 

Y 

2 

2 

X )/4 
2 


Y NU 
2 

X )/4 
2 


Y )/4 
1 


2 

X 

3 

X )/4 
2 


2 


2 


2 


2 


I LUJ[STRAT  I VC  CXAUrcCS 


- (Y  NU  - 2 Y 

Y 

NU 

+ Y 

NU 

- Y 

+ 

2 Y Y - 2 

X 

+ 4 X 

X 

- Y 

3 1 

3 

1 

3 

1 3 

3 

1 

3 

1 

. / 

2 

- 2 X )/4 

1 

HOIJ  G,  COL  1 

- (2  X Y NU  - 

2 > 

( Y NU 

- ^ 

f X i 

NU 

+ Y X fJU 

- X Y NU 

2 3 

1 

1 

3 

2 3 

1 3 

2 2 

+ 2 

X 

Y 

NU  - 

Y 

X NU 

+ 

Y X - Y 

X 

- X Y 

+ 

Y X ) /4 

1 

■ 2 

1 

2 

2 3 1 

3 

2 2 

1 2 

ROU  G,  COL  2 

(2  X Y NU  - 2 

X 

Y 

NU  - 

Y 

X NU 

+ 

Y X NU  + 

X 

Y NU  - 

2 

Y X NU 

2 3 

1 

3 

2 

3 

1 3 

1 

2 

1 2 

+ 

X 

Y NU 

+ 

Y X - Y 

X 

- X Y 

+ 

X Y )/4 

1 

1 

2 3 1 

3 

1 2 

1 1 

Rf.nj  G,  COL  3 

(X  - 

X 

) (Y 

Y ) (NU  + 1) 

2 

1 

2 

1 

4 


ROU  G,  COL  4 

2 

-(Y  Y NU-Y  Y NU-Y  NU  + Y Y NU-Y  Y + Y Y - 2X  K 
23  13  2 12  2313  23 

2 2 

+ 2X  X + Y - Y Y + 2X  - 2X  X)/4 
1 3 2 1 2 2 1 2 


non  G.  COL  5 

2 

(Y  Y NU-Y  Y NU-Y  Y NU  + Y NU 
2 3 13  12  1 


Y Y + Y 
2 3 1 


+ 2 X 


Y - 2 X 
3 2 

2 


- Y - 


X 

3 


+ 2 X 


X 


Y 


X 


2 

2 X )/4 


Il.l.USTRATIVK  EKAIHIES 
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I 


1 'J  1 2 12  1 1 


ROU  B/  COL  G 

2 2 2 2 2 
-(Y  NU-2Y  Y NU  + Y NU-Y  + 2Y  Y - 2X  +4X  X - Y 
2 12  1 2 12  2 121 

2 

- 2 X )/4 
1 


< Display  option  2 produces:  > 

The  options  for  OUTPUTTING  the  STIFFNLSS  I1ATRIX  are 
(.1)  Upper  triangular  portion  in  algebraic  format 
(2)  FORTRAN  CARD  II1AGE  format 
ENTER  the  NUMBER  ASSOCIATED  uith  YOUR  SELECTION. 

? 

2; 

Enter  the  FORTRAN  ARRAY  NAME  for  the  STIFFNESS  MATRIX 
? 

STIFF; 

1 

Multiply  each  coefficient  by 

-E*T/((X(2)*Y(3)-X(l),vY(3)-Y(2)*X(3)+Y(l)>vX(3)+X(l)*Y(2)-Y(l)*X(2) 

1 )*(NU-l),v(NU+l)) 

The  STIFFNESS  MATRIX  is 

STIFF  (1,1)  = -(X(3)**2*NU-2*X(2)*X(3)*NlJ+X(2)**2*NU-2>,vY(3)-.v*2+4*Y(  ^ 

1 2)  *Y  (3)  -X  (3)  >-nv2'l2*X  12)  »vX  (3)  ••2VfY  (2)  )V>v2"X  (2)  **2)  /4 

STIFF(2,1)  n (X(3) **2*NU-X ( 2 ) *X ( 3 ) ;vNU- X ( 1 ) *X ( 3 ) *NIJ4 X ( 1 ) *X ( 2 ) *NU -2* 

1 Y (3)  :-ov2  i2*Y  (2)  *Y  (3)  +2*Y  (1 ) VrY  (3)  -X  (3)  **24X  (2)  *X  (3)  +X  (1 ) *X  (3)  -2*Y 

2 (l)*Y(2)-X(l)*X(2))/4 

STIFF  (3, 1 ) = - (X  (3)  **2*NU-2*X  (1  )*X  (3)  *Nli+X  (1)  **2*NU-2*Y  (3)  ,-nv2+4*Y  ( 

1 1 ) ,vY  (3)  -X  (3)  **24.2*X  (1 ) ,vX  (3)  -2*Y  (1 ) ,v.v2-X  (1 ) ,v,>2)  /4 

STl  FF  (4 . 1 ) - - (X  (2)  *X  (3)  .-ffJU-X  ( 1 ) *X  (3)  *NU-X  (2)  *v.-2*NU4X  (1 ) ,vX  (2)  *NU-2 

1 *Y  ( 2 ) -.vY  ( 3 ) +2*Y  ( 1 ) *Y  ( 3 ) -X  ( 2 ) *X  ( 3 ) 4 X ( 1 ) *X  ( 3 ) +2* Y ( 2 ) ,-nv2-2* Y ( 1 ) >vY  ( 2 

2 )4X(2)>v*2-XI1)>vX(2))/4 

STI  FF  (5, 1 ) = (X  (2)  >.vX  (3)  >vNU-X  (1 ) -.vX  (3)  -.vNU-X  (1 ) *X  (2)  vrNLNX  (1 ) ,v-v2.vNU-2vf 

1 Y (2)  *Y  (3)  42*Y  (l),vY  (3)  -X  (2)  *X  (3)  +X  (1 ) *X  (3)  42*Y  (1 ) ,vY  (2)  4X  (1 ) *X  (2) 

2 -2*Y  ( 1 ) ,-nv2-X  ( 1 ) **2 ) /4  | 

STIFF  (G.  1)  - -(X(2)**2*NU-2*X(l)*X(2)*NU4X(l),v*2v.NU-2*Y(2)**24-4*Y(  *■ 

1 1 ) ,vY  (2)  -X  (2)  vnv2  (2*X  ( 1 ) *X  (2)  -2*Y  ( 1 ) .•ov2-X  ( 1 ) ,v*2)  /4 

STIFF(7,1)  - -(X(3)-X(2))*(Y(3)-Y(2))*(NU+l)/4 

STIFF(S,1)  - (X(3)*Y(3)VfNU-2*X(2)>vY(3)>vNU+X(l)*Y(3),vNU+Y(2)*X(3)»vN  L 


71 


II.LlJliTRATIVC  rXAIin.FS 

1 ( 1 ) >VK  ( 3 ) )VMl.l-X  ( 1 ) .vY  C:^ ) *Nlj42>vY  ( 1 1 >vX  (L’ ) ,vNlJ+X  ( J ) »vY  ( ) -X  ( 1 ) ,vY  ( 3 

2 )-Y(2)>vX(3)iK(]),vY(2))/A 

ST  1 (T  (3,1)  » (X (2)  *Y  (3) ivMLi-X  ( 1 ) >vY  (3)  >\MI  -2*Y  (2) *X (3)  *MU4  2*Y  ( 1 ) >vX  (3) 
1'  *NU4X  (2)  *Y  (2)  *NIJ4X  (1 ) *Y (2)  >vNl.l-2*Y  (1 ) ,vX  (2)  *NU  -X  (2)  *Y  (3)  4X  (1 ) *Y  (3 
• 2 )4X(2)*Y(2)-X(1),vY(2))/4 

ST  1 FF  (10, 1 ) » - ( Y (3)  •.v*2*Nl.l-2-.vY  (2)  *Y  (3)  •.vNU4Y  (2)  **2*IJU-Y  (3)  •;ov24-2-.vY  (2 
1 ) ,vY  (3)  -2-.vX  (3)  iv*244*X  (2)  -.vX  (3)  -Y  (2)  Vnv2--2*X  (2)  >v*2)  /4 

S T I FF  ( 1 1 , 1 ) - ( X ( 3 1 -.vY  ( 3 ) d JU4X  ( 2 ) :vY  ( 3 ) ,vNl  l-2*X  ( 1 ) *Y  ( 3 ) ,vNlJ-2-.v Y ( 2 ) ,vX  ( 3 

1 ) 4 Y ( 1 ) ,vX  ■(  3 ) -.vNlU  2V.-X  ( 1 ) ,vY  ( 2 ) Y ( i ) *X  ( 2 ) *RU4  X ( 3 ) *Y  ( 3 ) - X ( 2 ) ,v Y ( 

2 3)  - Y (1 ) ,vX  (3)  +Y  (1 ) ,vX  (2) ) /4 

STirF(12,l)  = -(X(3)~X(l))*(Y(3)-Y(l)).v(NU4l)/4 

STIFF (13, 1 ) - - (X (2) >vY (3) VrMU~X (1 ) >vY (3) >vNli-2*Y (2) .vX (3) >vMl.l4  2-,vY (1 ) >vX ( 

1 3)  *NU4  2-,vX  (1 ) *Y  (2)  .vN|l-Y  (1 ) *X  (2)  -.vNtl-X  ( 1 ) *Y  (1 ) ,vNll.  X (2)  ,vY  (3)  4X  (1 ) ,vY 

2 (3)4Y(1)*X(2)-X(1)*Y(1) )/4 

STIFF  (14, 1)  = (Y  (3)  •,vi-.-2*NU-Y  (2)  *Y  (3)  •;.NU-Y  (1)  >vY  (3)  -.vNl^Y  ( 1 ) ;vY  (2)  -.vNU-Y 

1 ( 3 ) •,v*24- Y (2 ) -.vY  ( 3 ) i Y ( 1 ) ,vY  ( 3 ) -2*X  ( 3 ) **24  2*X  ( 2 ) *X  ( 3 ) 42*X  ( 1 ) *X  ( 3 ) - Y ( 

2 1)*Y(2)-2*X(1)*X(2) )/4 

ST  IFF  (IS.  1 ) = - (Y  (3)  **2*M1J-2*Y  (1 ) *Y  (3)  *NIJ4Y  (1 ) **2*Nl.l-Y  (3)  **24-2*Y  ( 1 
1 ) *Y (3) -2*X (3) **244*X ( 1 ) *X (3) -Y (1 ) **2-2*X (1 ) **2) /4 

S T I FF  ( J G . 1 ) - - ( 2*X  ( 2 ) * Y ( 3 ) *f  lU  ••2*X  ( 1 ) * Y ( 3 ) *Nl J ••  Y ( 2 ) *X  ( 3 ) *Ml J + Y ( 1 ) *K  ( 

1 3)  *NU-X  (2)  *Y  (2)  *tni42*X  (1)  *Y  (2)  *N(l-Y  (1)  *X  (2)  *RIJ4Y  (2)  *X  (3)  -Yd)  *X 

2 (3)-X(2)*Y(2)4Y(l)*X(2))/4 

ST  I FF  ( 1 7 , 1 ) - ( 2*X  ( 2 ) *Y  (3 ) *NU- 2*X  ( 1 ) *Y  ( 3 ) *NU-  Y ( 2 ) ^vX  ( 3 ) *NIJ4 Y ( 1 ) *X  ( 3 

1 ) -.vNl.UX  ( 1 ) *Y (2) *NI.I-2*Y (1 ) *X (2) *NU4X (1 ) *Y (1 ) *MJ+Y (2) *X (3)  -Yd) *X ( 

2 3)-X(l)*Y(2)4Xd)*Y(l))/4 

STIFFdS.l)  = -(X(2)-Xd))*(Y(2)-Y(l))*(NU4l)/4 

STIFF  (19. 1 ) - - (Y  (2)  *Y  (3)*NU-Y  (1)  *Y  (3)  *NU-Y  (2)  **2*NU+Y  (1 ) *Y  (2)  *NU- 

1 Y (2)  *Y  (3)  4 Y (1 ) *Y  (3)  -2*X  (2)  *X  (3)  42*X  (1 ) *X  (3)  +Y  (2)  **2-Y  (1 ) *Y  (2)4-2 

2 *X(2)**2-2*X(1)*X(2) )/4 

ST  I FF  ( 20 , 1 ) - ( Y ( 2 ) *Y  ( 3 ) vrfllJ- Y ( 1 ) *Y  ( 3 ) *N(.I-  Y ( 1 ) *Y  ( 2 ) *NU4  Yd)  **2*N1.I - Y 

1 (2)  *Y  (3)  4 Y (1 ) *Y  (3)  •.2*X  (2)  *X  (3)  42*X  ( 1 ) *X  (3)  +Y  (1 ) *Y  (2)  4-2*X  (1 ) *X  (2 

2 ) - Y (1 ) **2-2*X ( 1 ) **2 ) 74 

ST  1 FF  (21 . 1 ) = - ( Y (2)  **2*NIJ-2*Y  ( 1 ) *Y  (2)  *NU4Y  ( 1 ) **2*NU-Y  (2)  **24-2*Y  ( 1 
1 ) *Y  (2)  -2*X  (2)  **24 4*X  (1 ) *X  (2)  -Y  (1 ) i-ov2-2*X  ( 1 ) •.Wf2)  74 


3.3.3  Comporicon  cind  Verification 

< The  verification  for  this  case  consists  of  comparing  the  results  of 
both  formulations  with  a fiducial  vector,  "SUDEV2",  adapted  from 
Rubenste  i n (141  , pages  107-110.  The  element  represented  by  R(JBEV2  is 
located  in  space  such  that:  Xd)=0,  Y(1)“0.  Y(2)'‘0,  Xf1IN“0,  YniN  = 0, 
Xf1AX»>X(2),  YilA>(»Y(3).  For  each  formulation,  the  calculated  vector  is 
evaluated  for  these  boundary  conditions  and  subtracted  from  tlie 
fiducial  vector  thys  producing  a null  vector.  > 

< For  the  generalized  coordinate  formulation:  > 

(Dl)  tOSK,  AKIGI 


(C2)  L0ADFILE('/3G.VALS); 


ILl.UMIVATIVK  I^XAUPLPS 


rSG  VAL5  DGK  AKIG  being  I ceded 
I Oc»H  i tig  done 
time**  2098  mcec. 

(LI2)  /'  D0f€ 

(C3»  FUNCTIONAL tl. 11; 

time-  1 fntsec. 

2 -22  2 
(03)  -(X  NIJ  - 2 X X NIJ  + X NU-2Y  +4Y  Y - X + 2X  X 
3 23  2 3 233  23 

2 2 

- 2 Y - X 1 (XI1IN  - XMAX)  (YIIIN  - Yt1AX) /2 
2 2 

(C4)  FACTOR: 
time’’  0 msec. 

2 

(D4)  - E T/(2  (X  Y - X Y - Y X + Y X ^ X Y - Y X ) 

23  13  23  13  12  12 

(NU  - 1)  (NU  + 1) ) 


(C5)  F7L : FUNCT I ONALivFACTOR# 
t i tne=  255  msec. 

(CG)  FTLtl.ll; 
t i me»  1 msec. 

2 2 2 2 

(DG)  E(X  NU-2X  X NU  + X NU-2Y  + 4Y  Y - X + 2X  X 

3 *23  2 3 233  23 

2 2 

- 2 Y - X ) T (Xm.N  - XOAX)  (YI1IN  - YMAX) 

2 2 

2 

/(4(X  Y - X Y - Y X + Y X + X Y - Y X)  (NU-1) 

23  13  23  13  12  12 

(NU  +1)) 

< Enpression  (DG)  is  element  (1,1)  of  the  stiffness  matrix  generated 
for  the  generalized  coordinate  formulation.  This  matrix  is  stored 
under  the  name,  "FTL".  > 

(C7)  LOAOFILE(CST,FID): 

CST  FID  DSK  AKIG  being  loaded 
loading  done 
time"  1578  msec. 

(D7) 


DONE 
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II.I.USTIVATI  V{-;  LKAIII  lbS 

(C8)  RI.IBF.V2[1.1]  ; 
tiinf!<!  1 mrjec. 

2 

/ (X  - X ) (1  - NU) 

3 2 2 

E (— - + Y ) T 

2 3 

(DS)  r - 

2 

2 X Y (1  - NU  ) 

2 3 

I 

< Expression  (D8)  is  element  (1,1)  of  the  stiffness  matrix  in 
Rubens te i n’ s example.  This  matrix  is  stored  as  ".RUBEV2".  > 

< The  genera  I I lied  coordinate  stiffness  iitatrix  is  evaluated  at  the 
boundary  conditions:  > 

(C3)  FTL:  EV  (FTL,  XlllN-R,  XI1AX=X  (23 . YtllN^R,  YI1AX=Y  (3) , X (1)  =0,  Y (13  «R.  Y (23  =0 
)8 

t i mo=  3771  msec. 

(CIO)  FTLd.l); 
time»  1 msec. 

2 2 2 2 2 
E(X  NU-2X  X NU  + X NU-2Y  - X + 2X  X - X)T 
3 23  2 33  232 

(DIO)  

4 X Y (NU  - 1)  (NU  + 1) 

2 3 

(Cll)  FTL[1.13-RUBEV2(1,138 
time*’  9 msec. 

(C12)  RATSI nr  (’/() ; 
time=  110  msec. 

(D12)  0 

(Cl 3)  FTL-RLIBEV28 
time**  G77  msec. 

(Cl  4)  RATSIf1P(%); 
time**  2879  msec. 

(D14)  MATRIX ((03  . (03  , (03  . (03  , (03  . (03  . (03  . (03  , (03  . (03  , 

(X  + X ) NU  + X - X 


3 2 3 2 

I-  ((X  + X ) Y E NU  + ((X  - X ) Y - 2 Y ( )) 

323  3233  2 


2 

E)  T/(4  X Y NU  - 4X  Y)  3 . (03  , (03  . (03  . (0).  [03  . (03  . (03  . 
2 3 2 3 


ILLU'JTRATIVL  EXAUPLES 
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[(».  [0],  [0]) 

< Expr/etision  (DIA)  contoino  the  term  by  term  difference  Ijctni’en  tlio 
niotrix  generated  by  *iv.v-.v*  and  that  in  Rubenetein’s  example.  For  an 
0 undetermined  reason,  flACSYUA  is  unable  to  recoynire  that  element  11 

of  thin  vector  is  t'C|iial  to  rero.  To  overcomn  this,  tlie  element  is 
presented,  in  the  follouing  linen,  in  a format  tliat  makes  thin 
evaluation  obvious  by  inspection.  > 


(CIS)  Xlll.l]; 
time"  1 msec. 

(015)  - ((X  + X)Y  ENU+((X  - X)Y 
3 2 3 3 2 3 

(X  + X ) NU  + X - X 
3 2 3 2 2 

_ 2 Y ( ))  E)  T/(4  X Y NU  - A X Y ) 

3 2 .2  3 2 3 

(CIG)  NUIK'/C); 
t i mc"  11  msec. 

(DIG)  - ( (X  + X ) Y E NU  + ((X  - X ) Y 
3 2 3 3 2-3 


(X  4 X ) NU  + X - X 
3 2 3 2 

_ 2 Y ( ))  E)  T 

3 2 


(C17)  EV('X.E=1.T-1); 
time"  28  msec. 

(X  + X ) NU  + X - X 
3 2 3 2 

(017)  - (X  4-  X ) Y NU  4-  2 Y ( ) - (X  - X ) Y 

323  3 2 323 

(CIS)  EXPAND  (10; 
time"  25  msec. 

X NU  X NU  X X 
3 2 3 2 

(018)  - X Y NU-X  Y NU4-2Y  ( + 4- ) - X Y 

33  23  32  222  33 


4-  X Y 
2 3 


< Expression  (018)  is  a factor  of  the  numerator  of  the  11th  element. 


After  expanding  the  term  in  parenthesis,  it  can  be  seen  that  the 


fol  lowing  pairs  of  terms  cancel:  1 and  3,  2 and  A,  5 
8.  Thus  element  11  is  indeed  zero  and  the  generated 
is  equivalent  to  the  one  presented  by  Rubenotein.  > 


and  7,  ;ind  G and 
St i f fncss  matr  i x 
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Il.LUr>TRAT|V|-;  RXAIII-'LES 
(r.;i9)  CL05E(-ILE  (%VEn,C3G) ; 

< For  tlie  i !iopnramo  tr  i c fornui  I nt  i on  an  identical  comparison  is  made 
Ijctwroh  the  generated  matrix  and  Rubenstein’ s:  > 

(Dl)  tDSK,  AKIG] 

(CZ)  L0AnFlLE(y,3I.VAl.S);  • 

%3I  VAl.S  D5K  AKIG  being  loaded 
load! tig  done 
time"  1922  msec.  ’ 

(D2)  DOUE 

(C3)  FTL:FUNCTIONAl.*FAGTOR» 
time"  174  msec. 

(C4)  FTLQ.l]! 
tinr*"  1 msec. 

2 2 2 2 
(D4)  E (X  NIJ  - 2 X X NLI  + X NU  - 2 Y + 4 Y Y - X + 2 X X 
3 .23  2 3 233  23 

2 2 

-2Y  - X)T/(4(X  Y - X Y - Y X + Y X + X Y - Y X) 

22  2313231312  12 

(NU  - 1)  (NU  + D) 

< Expression  (D4)  is  element  (1,1)  of  the  stiffness  matrix  generated 
for  the  isoparametric  formulation.  > 

(C5)  LOAUFILE(CST.FID); 

C5T  FID  DSK  AKIG  being  loaded 
I oad i ng  done 
time*'  197G  msec. 

(05)  DONE 


(CG)  RUBEV2I1.1]; 
time"  1 msec. 


106) 


2 

(X  - X ) (1  - NU) 


3 2 2 

E ( + Y ) T 

2 3 


2 X Y (1  - NU  I 
2 3 
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ILLUSTHATIVK  tXAIII-LES 
concJitiono.  > 

(C7)  rTL:EV(FTL,Xll]uO,Ytl]-O.Yl2]«0)a 

You  hove  run  out  of  LIST  space. 

Do  you  wont  more? 

Tupe  Al.L;  NONE:  a level-no.  or  the  name  of  a space. 

LIST; 

time=  3751  msec. 

(CP>)  FTLd.l]; 
tiwe=  1 msec. 

2  • 2 2 2 2 
E(K  NU-2X  X NU  + X N(J-2Y  - X + 2X  X - X)T 

3 23  2 33  232 

(D8)  

4  X Y (NU  - 1)  (NU  + 1) 

2  3 

(C3)  FTL[1,11-RLIBEV2[1,1]; 
time*=  9 msec. 

2 2 2 2 2 
E(X  NU-2X  X NU  + X NU-2Y  - X + 2X  X - X)T 

3 23  2 33  23. 2 

(09)  

4  X Y (NU  - 1)  (NU  + 1) 

2 3 


C. 

(X  - X ) (1  - NU) 


3 2 2 

E ( - + Y ) T 

2 3 


2 

2 X Y (1  - NU  ) 

2 3 

(CIO)  RATS  I (IF  (-4); 
time”  111  msec. 

(DIO)  0 

(Cll)  FTL-RUBEV2« 
time”  GG3  msec. 

(C12)  RATSini’di); 
time”  2822  msec. 

(D12)  MATRIX!  f0),  fOI , tOI . tO] , 101.  101.  101.  10).  tO) , 103, 

(X  + X ) NU  + X - X 
3 2 3 2 


[-((X  +X)Y  ENU+((X  - X)Y  - 2Y  ( )) 


r 
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2 

E)  T/(A  X Y NIJ  - A X Y )].  [0],  [0],  (0],  10],  [01.  [0],  f0] , 

, 2 3 2 3 

[0] . [0] , [0] ) 


< Once  cicjiiin,  the  difference  between  the  generated  and  fiducial 
matrices  is  determined.  Element  11  of  thin  vector  is  handled  as  in 
the  general  i;:ed  coordinate  case.  > 

I 

(C13)  ni2[ll,l]; 
t i ine=  1 msec. 

(013)  - ((X  + X)Y  ENU+((X  - X)Y 


3 2 3 3 2 3 

(X  + X ) NU  + X - X 
3 2 3 2 2 

_ 2 Y (— ^ ))  E)  T/{f(  X Y NU  - 4 X Y ) 

3 2 2 3 2 3 

(C14)  NUH()f); 
time-  11  msec. 

(D14)  -((X+X)YENU+((X-X)Y 
3 2 3 3 2 3 


(X  + X ) NU  + X - X 


3 2 3 2 

- 2 Y ( ))  E)  T 

3 2 


(C15)  EV(%,E-1,T-1) ; 
time-  23  msec. 

(X  + X ) NU  + X - X 
3 2 3 2 

(D15)  -(X  + X)Y  NU  + 2Y  ( ) - (X  - X ) Y 

323  3 2 323 


(Cl 6)  EXPAND  ('/(){ 

time-  25  msec. 


1 
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r 


• ^ 

\ 


t 


3.4  Four  Node  (luadri  lateral 

The  ncxTal  locations  and  the  decrees  of  freedom  for  the  four  node 
quadrilateral  are: 


Z 


J 


3.4.1  Generalized  Coordinate  Formulation 
* SYSTEM  INITIALIZATION  * 


UELCOME  TO  *****  VERSION  1.0 

It  is  now  UEDNESDAY  DECEMBER  7,1 S77  11:28:58 

The  current  file" is  ISYINIT,  FCN] 

The  current  device  and  username  is  IDSK,  AKIGI 


Report  problems  to 
ALAN  R.  KORNCOFF 
DEPT.  OF  CIVIL  ENGINEERING 
CARNEG I E-MELLON  UN I VERS  I TY 
CMU-10A,  AXIG 


4 


Terminate  all  input  with  a SEMICOLON  - 
Input  for  HELP 


Input  your  LOGIN  NAME 
? 

AKIG; 

Is  AKIG  correct  ? (YESj  or  NO;) 
? 

YES; 

GREETINGS  AKIG 


* METHOD  SELECTION  * 

The  available  formulations  include 

(1)  THE  ISOPARAMETRIC  METHOD 

(2)  THE  GENERALIZED  COORDINATE  METHOD 


ILLUSTRATIVE  rxAiin.rs 


Pl(>c)fe  entcjr  tliu  number  of  tlie  inetliod  chor.pn  (J  OR  2). 

? 

2; 

GENLRAIJZED  COORDINATE  FORflULATION 

* GENL-RALIZEO  COORDINATE  FORtIULATION  EXECUTIVE  * 

* GENERALIZED  COORDINATE  FORMULATION  INITIALIZATION 
Do  you  141 sh  to  set  DUMP,  BREAK  or  TRACE  POINTS 

NO; 

* PROBLEM  PARAMETER  SPECIFICATION  - GENERALIZED  COORDINATE 

Input  the  NUMBER  OF  ELEMENT  NODES 
? 

4; 

Input  the  NUMBER  OF  DEGREES  OF  FREEDOM  PER  NODE 
? 

2; 

Input  the  vector  of  the  NAMES  OF  THE  GLOBAL  COORDINATES 
There  should  be  2 elements 
ELEMENT  1 » 

K; 

ELEMENT  2 « 

Y; 

Input  the  vector  of  the  NAMES  OF  THE  DISPLACEMENT  VARIABLES 
There  should  be  2 elements 
ELEMENT  1 ™ 

U; 

ELEMENT  2 » 

V; 


,v  SHAPE  FUNCTION  PROCESSOR  - GENERALIZED  COORDINATE  * 

ENTER  the  terms  of  the  SHAPE  FUNCTION  ordered  from 
GENERALIZED  COORDINATE  1 through  coordinate  4. 

ELEMENT  1 » 

1; 

ELEMENT  2 - 
X; 

ELEMENT  3 - 
Y; 

ELEMENT  4 - 
X,vY; 

< This  specification  represents  the  displacement  functions; 

u » al  + a2*x  + aoVry  + o4*x*y 


ILUJSTRATIVf-:  exauples 


V •»  a5  + ciGVfX  + (j7vfy  + iiSvrx*i,) 
in  uhicli  «!  through  08  ciro  the  generalized  coord  i no  lea.  > 

SHAPE/FUIJCTION  flODiFICATION 
The" OPTIONS  are 

(1)  DISPLAY  THE  SHAPE  FUNCTIONS 

(2)  MODIFY  THE  SHAPE  FUNCTIONS 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  uith  the  CHOSEN  OPTION 
? 

1; 


The  terma  of  the  SHAPE  FUNCTION  ore 
Term  1 

1 


Term  2 


X 


Term  3 


Y 


Term  4 


X,vY 


Enter  the  NUMBER  ASSOCIATED  uith  the  CHOSEN  OPTION 
? 

3{ 


* B MATRIX  DATABASE  GENERATION  * 

The  OPTIONS  for  specifying  STRAIN  COMPONENTS  are 

(1)  USER-SUPPLIED  VALUES 

(2)  LIBRARY  VALUES 

Please  ENTER  the  NUMBER  ASSOCIATED  UITH  YOUR  SELECTION 
? 

2; 

The  LIBRARY  OPTIONS  for  SPECIFYING  STRAIN  COMPONENTS  ore 

(1)  USER -SUPPLIED  VALUES 

(2)  ONE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 

(4)  PLANE  STRAIN 

(5)  AX I SYMMETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 3D 
Please  ENTER  the  NUMBER  ASSOCIATED  UITH  YOUR  CHOICE 
? 

3; 


* MATERIAL  PROPERTIES  SELECTION  * 


ILLUSTRA  I IVC  L'KAIII'LES 


I liR  options  for  tho  selnction  of  thn  IIAILRIAI.  i'KRI-'fcH  1 1 tS  (lAIHlK  arc* 

(1)  UStR-SUPPLIEO  MATRIX 

(2)  LIBRARY  MATRIX 

Plea, so' enter  the  NUMBER  ASSOCIATED  UITH  YOUR  SELECTION 
? 


The  LIBRARY  OPTIONS  for  SPPCIFYING  MATERIAL  PROPERTIES  are 

(1)  USER-SUITLIEO  VALUES 

(2)  ONE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS  ' 

(A)  PLANE  STRAIN 
(S)  AX  I SYMMETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 30 
Please  ENTER  the  NUMBER  ASSOCIATED  UITH  YOUR  CHOICE 
? 


* AUXILIARY  TERM  PROCESSOR  A 

ENTER  ELEMENT  VOLUME  MODIFICATION  FACTORS  OR  AUXILIARY  TERMS 
TYPE  ’END;’  to  TERMINATE 
? 

T; 

? 

END; 

AUXILIARY  TERM  MODIFICATION 
The  OPTIONS  are 

(1)  DISPLAY  THE  AUXILIARY  TERMS 

(2)  MODIFY  AN  AUXILIARY  TERM 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 
? 

3; 


* INTEGRATION  LIMITS  - GENERALIZED  COORDINATE  * 

A INVERSION  OF  MATRIX.  A * 

,v  B MATRIX  GENERATION  - GENERALIZED  COORDINATE  * 

Vf  INTEGRATION  PROCESSOR  - GENERALIZED  COORDINATE  * 

A DISPLAY  PRE-PROCESSOR  - GENERALIZED  COORDINATE  a 

< To  avoid  oxcocrJing  the  list  ntorage  capacity  of  this  experimental 
veroion  of  MACSYMA,  the  expression  simplification  function  in  the 
display  pro-processor  phase  was  bypassed.  > 


A DISPLAY  PROCESSOR  A 


1 
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II.LIXTIVATJVt  LXAni'LCS 

< As  th*f  ou|)ut  is  voluniinuuR,  only  thn  ficst  coofficent  is  ci  i rip  I ayuci 
her«.  > 

The  ,o|/tiono  for  OUTPUITING  the  STIFFNIISS  I1ATR1X  nrn 
« (1)  Upper  triangular  portion  in  algebraic  format 

(2)  FORTRAN  CARO  II1AGE  format 
ENTER  the  NUMBER  ASSOCIATED  uith  YOUR  SELECTION. 

? 

1; 

Multiply  each  coefficient  by 

i - E*T/(EXPT(X  *Y  *X  *Y  - X vrY  *X  *Y  - Y *X  *X  »vY  + Y *X  >vX  »vY 

2344134423A41344 

+ X VfY  *X  >vY  - Y *X  *X  *Y  - X *X  *Y  *Y  + X *X  »vY  *Y  + X *Y  *X  >vY 

1244  1244  2334  1334  2234 

- X *Y  VrX  *Y  - X *X  >vY  *Y  + X *Y  VrY  + Y *K  *Y  ,vX  - Y VrX  ,vY  *X 

11  3 4 1224  1124  2334  1334 

- X *Y  *Y  *X  + X -.vY  v«Y  *X  -f  Y ,vX  *Y  ,vX  - X *Y  *Y  aX  - X aY  aX  aY 

2234  .1134  1224  1124  1233 

+ Y aX  VfX  >vY  + X >vX  aY  v»Y  - X tvY  aX  aY  - Y aX  aY  v*X 

1233  1223  1123  1223 

+ X aY  aV  aX  , 2)a(NU  - DaINLI  ♦ 1)1 
112  3 

< The  MACSYMA  function,  EXPT,  is  usetl  to  represent  a value  cciual  to 
ite  first  argument  raised  to  the  exponent  given  by  ito  second 
argument. > 

RON  1,  COL  1 


1 


(X  a(Y  _ y ) . X aY  - Y aIX  - X ) ♦ Y aX  ) 

24  3 34243  34 

3 3 3 3 

XMINaYMIN  XMAXaYMIN  (1  - NUIaXMIN  aYIUN  (1  - NUIaXMAX  aYMIN 

a( — - +■  - ------------ — - — 

3 3 G 6 

3 3 3 

XMINaYMAX  XMAXaYMAX  (1  - NLDaXMIN  aYIIAX 

+ 

3 3 6 

3 

(1  - NUIaXMAX  aYMAX 

+ ) + 2a(X  a(Y  - Y ) - X aY  - V a(X  - X ) 

G 243  34243 
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Il.l.UlURATI  VH  EXAIIFLCS 

I V *X  ),v(V  *(>;  *Y  - K VrY  ) - X *Y  -.vlY  - Y ) - Y *X  VrY  + X *Y  *Y  ) 
3 A 2 4A  33  22A  3 344  334 

2 2 2 2 
YI1AX  YI11N  YI1AX  YIIIN 

*(XI-1AXVf(-— ) - XniNVf(— )) 

2 2 2 2 

2 

+ (-  X Vf(X  *Y  - X ,vY  ) + X ,vX  >vY  -t  X *Y  *(X  - X ) - X *Y  *X  ) 

2 44'33  344  224  3 334 

(1  - NLDvrYHAX  (1  - NLI)*YI1IN 

,v(XI1AX,v( ) 

2 2 

(1  - NLOvrYHAK  (1  - NLDuYtllN 

- X|11N*( )) 

2 2 

*■  2*(X  vv(Y  - Y > - X *Y  - Y *(X  - X ) + Y *X  ) 

2 4 3 34  24  3 34 

Vf(-  X ,v<X  *Y  - X *Y  ) + X *X  *Y  ♦ X *Y  *(X  - X ) - X *Y  *X  ) 

244  33  344  224  3 334 

2 2 2 

(1  - NLDivXtllN  »vYniN  (1  - NU)  *xnAX  *YniN  (1  - NU)V(XniN  VrYMAX 



4 4 4 

2 

(J  - l^)rtXt1AX  ivYHAX 

+ ) 

4 

• 

2 

+ (Y  *(X  aY  - X *Y  ) - X >vY  *IY  - Y ) - Y *X  *Y  + X *Y  *Y  ) 
244  33  224  3 344  334 

»v(XI1AX-.v(YnAX  - YniN)  - XniN.v(YnAX  - VniN)) 


3.4.2  Calibration  of  the  Generalized  Coordinate  Results 

Since  there  are  no  tabulated,  c I osetl- form  CKpress i ons  against 
iihich  the  generated  coefficients  could  be  compared,  a calibration 
scheme  uas  devised.  The  transverse  deflection  of  the  free  end  of  a 
cantilever  beam,  as  calculated  using  the  generated  stiffness  matrix, 
is  compared  to  the  value  determined  using  numerical  integration  as 
tabulated  in  Cook  I15I  . The  cases  of  a moment  and  a transverse 


1 
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< FACTOR  ond  FUfCTI ONAl.  arc  evaluated  at  the  geometric  boundary 
conditions  and  for  the  specified  material  properties.  > 


(C7)  FACTOR; 
tiine=  R msec. 
(D7) 

(C8)  FUNCTIONAL; 
tiino”  0 msec. 


1.70GGGGBGE-3 


t A32.231G7  ] 

[ ] 

t - 3B9.731G7  ] 

[ ] 

432.231G7 


21G.14GS28 


153.G45S2S 


432.231G7 


153.G45828 


ll.LU'oTIMTI  VI-:  liXAIH-LllS 
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I 

97.GBG2& 

] 

[ 

) 

1 

- 19.53125 

) 

I 

) 

1 

- 97.G5G25 

) 

I 

] 

[ 

19.53125 

] 

[ 

] 

[ 

- 528.G/i58/i 

] 

[ 

] 

t 

505.20834/) 

) 

[ 

] 

1 

1057.291GG 

] 

[ 

] 

I 

19.53125 

] 

[ 

) 

[ 

- 97.G5G25 

] 

I 

I 

[ 

- 19.53125 

1 

[ 

3 

[ 

97.G5625 

3 

I 

3 

[ 

505.208344 

3 

[ 

3 

[ 

- 528.G45S4 

3 

[ 

3 

[ 

- 1033.8541G 

3 

[ 

3 

t 

1057. 291 GG 

3 

< Each  coefficient  is  iniiltiplied  hy  the  common  multiple.  > 

(C9)  FUNCTIONAL; FACTOR, vFLllCT I ONALH 
tiine=  71  msec. 

< The  vector  is  expanclecl  into  a square  matrix,  > 

(C10)  NUf1:Htt 
t i me«»  1 msec. 

(Cll)  FOR  I THRU  8 00  FOR  J THRU  I DO  (NUrirNUIUl.FTL  II , J]  :FTL  [J.  1]  :FUN 

CTIONALINLIII.l]); 

timo»  701  msec. 

(Dll)  DONE 


(C12)  STIFF:GENnATRIX(FTL.S,8); 
time"  37  msec. 

[ 0.7377777SA  ) ( - 0.G3111111B  1 


1 


ILLUGTRATIVt  EXAIH'Lr.S 
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r ' 


[.  d.?G22?.221  ] 


(D12)  Col  1 « [ ] 

( - H.IKGGGGGGG  3 
[ ] 
t - 0.033333333  ] 
[ ] 
[ 0.1GBGG6GGG  ] 
[ 3 

[ 0.033333333  3 

[ - 0.3BSS3SS3  3 
[ 3 

t 0.2G223221  3 

I 3 

[ 0.73777778^.  3 

[ 3 

[ - 0.G31111115  3 
Col  3 = [ 3 Col  4 

[ 0.1GGGGGGGG  3 
t 3 

t 0.033333333  3 
[ 3 

[ - 0.1GGGGGGGG  3 
[ 3 

[ - 0.033333333  3 

[ - 0.1GG6GG66G  3 
[ 3 

[ 0.033333333  3 

[ 3 

[ 0.1GGGGGGGG  3 
[ 3 

t - 0.033333333  3 
Co  1 5 [ 3 Co  I G 

[ 1.80444442  3 

t 3 

[ - 1.7G4 44443  3 

[ 3 

t - 0.00222224  3 
t 3 

[ 0.SG222224  3 


[ 

t 

[ 

[ 

[ 

[ 

t 

Col  7 - [ 
( 
[ 
t 


0.1GGGGGGGG  3 
3 

- 0.033333333  3 

3 

- 0.1GGGGGGG6  3 

3 

0.033333333  3 

3 Col  8 

- 0.90222224  3 

3 

0.8G222224  3 


Col 


[ 

■ [ 
[ 
[ 
[ 
[ 
[ 

= [ 
[ 
[ 
t 
[ 
[ 
[ 
[ 

[ 

[ 

[ 

[ 

[ 

[ 

[ 

[ 

[ 

I 

[ 

t 

I 

( 

t 

t 

[ 

[ 

i 

[ 

[ 

( 

" [ 
t 
[ 
( 


[ - 0.3G88888S  3 

2-1  3 

[ 0.(333333333  3 
[ 3 

I O.IGGGGBGGG  3 
[ 3 

( - 0.033333333  3 
[ 3 

[ - 0.1G&GGGGGG  3 

0.26222221  3 

3 

- 0.3GSSSS88  3 

3 

- 0.G3111111B  3 

3 

0.737777784  3 

3 

- 0.033333333  3 

3 

- 0.1GGGGGGGG  3 

3 

0.033333333  3 

3 

0.1GGGGGGGG  3 

- 0.033333333  3 

3 

0.1GGG6GGGG  3 
3 

0.033333333  3 

3 

- 0.1G6GGGGGG  3 

3 

- 1.7G444443  3 

3 

1.80444442  3 

3 

0.SG2.72224  3 

3 

- 0.30222224  3 

0.033333333  3 

3 

- O.IGGGGGGGG  3 

3 

- 0.033333333  3 

3 

O.IGGGGGGGG  3 
3 

0.8B2.72224  3 

3 

- 0.902222.24  3 


d 
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ll.l.Ut'TRATIVC  FIXAIIPLIIS 


I ] 

[ ] 

t ] 

t - 1.7GV.4/i43  ] 


( 1 

( - 1.7B^iA-Vi/(3  ] 

[ ] 

[ 1.80AA^i^t^2  ] 


< The  miitrix  is  corroctecl  for  the  di  r>|:)  I iiccinent  Ijoundory  conditions  by 
eliminatiny  rous  and  columns  3,  4 and  S.  > 

{ C 1 3 ) S T I FFREDUCED : SUEM A T R 1 X ( 3 . 4 . 8 . S T I FF , 3 , 4 , 8 ) ; 
tiinn^  8 msec. 

t 0.737777784  ] [ - 0.031111115  ] 

I 11  1 

[ - 0.G31111115  1 ( ?i.l2nill?A  ] 

( 1 [ ] 

(D13)  Col  1 - t - 0.1GGGGGGGG  1 Col  2 = f 0.033333333  1 

[ 1 [ 1 

[ - 0,033333333  1 [ O.IGGGGBGGG  1 

[ 1 [ 1 

t 0.1BGGGGGGG  1 [ - 0.033333333  1 


Col  3 


Co  I 5 


1 - 0.1GGGGGGGG 

] 

[ - 0.033333333 

] 

[ 

1 

[ 

1 

[ 0.033333333 

1 

[ O.IGGGGBGGG 

] 

[ 

1 

[ 

) 

[ 1.80444442 

1 Col 

4=1-  1.7G444443 

] 

[ 

1 

[ 

1 

[ - 1.7G444443 

1 

[ 1.80444442 

3 

[ 

] 

[ 

3 

[ - 0.90222224 

3 

[ 0.8G222224 

3 

t O.IGGGGGBGG 

1 

[ 

] 

r - 0.033333333 

1 

[ 

] 

[ - 0.90222224 

1 

t 

1 

[ 0.8G222224 

1 

[ 

] 

[ 1,80444442 

] 

< The  reduced  matrix  is  inverted.  HACSYMA  logs  its  use  of  function 
RAT.  > 


(C14)  RESULT: 
RAT  replaced 
RAT  replaced 
RAT  replaced 
RAT  replaced 
RAT  replaced 
RAT  replaced 
RAT  replaced 
RAT  replaced 
RAT  replaced 


EV  (STlFFREOUCEDtt-1  .NIJIIER:  TRUE) ; 
0.737777784  by  lGG/225  = 0.73777778 
-0.G31111115  by  -142/225  - -0.B3111111 
-0.1GBGG6GGB  by  -1/B  - -0,1BGGGGGBG 
-0.033333333  by  -1/30  - -0.033333333 
O.IGGGGGBGG  by  1/G  - O.IGGGGBGGG 
-0.G31111115  by  -142/225  ■■  -0.63111111 
0.737777784  by  lGG/225  - 0.7377/778 
0.033333333  by  1/30  » 0.033333333 
O.IGGGGGBGG  by  1/B  - O.IGGGGGGGG 


i 


ii.i.U'STnAiivi-:  pxAiirixs 


RAT 

repi 

.•>ced 

RAT 

rep  1 

aced 

RAT 

repI 

aced 

RAT 

rep  1 

.'iced 

RAT 

repi 

need 

RAT 

repi 

aced 

RAT 

repi 

aced 

RAT 

rep  1 

aced 

RAT 

rep  1 

aced 

RAT 

rep  1 

aced 

RAT 

rep  1 

need 

RAT 

rep  1 

aced 

RAT 

rep  1 

aced 

RAT 

rep  1 

aced 

RAT 

rep  1 

aced 

RAT 

r ep  1 

aced 

-0.033333333  by  -1/30  - -0.033333333 
-0.lf.6GF.Ge0G  by  -1/G  " -0.  IGGGGGGGG 
0.033333333  by  1/30  - 0.033333333 
1.80A44A42  by  40G/225  = 1.80444445 
-1.7G444443  by  -397/225  « -1.7G444444 
-0.90222224  by  -203/225  « -0.90222222 
-0.033333333  by  -1/30  - -0.033333333 
0.1GGGGGGGG  by  1/G  - O.IGGGGGGGG 
-1.7G444443  by  -397/225  » -1.7G444444 
1.80444442  by  40G/225  = 1.80444445 
0.SG222224  by  194/225  » 0.8G2222224 
0.16GGGGGGG  by  1/G  = 0.16GG6GGGG 
-0.033333333  by  -1/30  - -0.033333333 
-0.90222224  by  -203/225  = -0.902.22222 
0.SG222224  by  194/225  « 0.8G2222224 
1.80444442  by  40G/225  » 1.80444445 


tiiiifix  4G9  msec. 


(D14) 

[ G. 355421 GG 
[ 

[ 3.G445783 
[ 

[ G. 77710843 
I 

t G. 52710843 
[ 

[ - 0.25 


3. G4 4 5783 
G.355421GG 

- G. 77710843 

- 7.02710843 

- 0.25 


5.77710843 
- G. 7771 0843 
4 G. 939729 
4G. 3855424 
0.55418719 


G. 5271 0843 
- 7.02710843 
46.3855424 
4G.5S5542 
0.2 


- 0.25  ] 

] 

- 0.25  ] 

] 

0.55418713  ] 

] 

0.2  ] 

] 

0.75418719  ] 
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< Deyree  of  freedom  5 is  represented  by  the  third  rou  of  the  reduced 
mntrix.  > 

< For  case  1,  unit  loads  are  applied  in  the  positive  and  negative 

direction  of  degrees  of  freedom  1 and  2 I'especti vely.  > 

(CIS)  DELTAl:RESULT(3.1]-RFSULTr3,2]; 
time»  3 msec. 

(015)  13.55421G3 

< For  case  2,  unit  loads  are  applied  in  the  positive  direction  of 

degrees  of  freedom  5 and  G and  in  the  negative  direction  of  degree  of 

freedom  7.  > 

(CIG)  DELTA2:RESULT  (3.3] -fRESULT  (3,41 -RESULT  (3.51; 
time»  4 msec. 

(DIG)  92.771085 


In  reference  (15)  Coo)^  scales  the  results  from  the  finite  element 
model  with  respect  to  a normalized  value  generated  by  applying  beam 


ILI.IISTUAIIVI-  i:XAI1l-LLS 


theory.  the  calculoted  defloctiono  of  the  free  end  of  the  conti  lever 
subjected  to  both  loading  conditions  are  tabulated  beloii.  The  values 
presented  are:  (a)  Cook's  scaled  result  from  the  finite  clement  model; 

(b)  Cook’s  normalized  value  from  beam  theory;  (c)  The  value  determined 
from  tlie  generated  stiffness  m.itrix  and  (d)  The  un-norma  I i zed  value 

caluculated  from  beam  theory. 

• 

CASE  1: 

(a)  Cook's  finite  element  model:  9.0 

(b)  Cook’s  normalized  beam  theory  value:  103.0 

(c)  Result  from  generated  stiffness  matrix;  13.55 

(d)  Un--norma I i zed  beam  theory  result;  150.0 

Scaling  the  results  form  the  generated  stiffness  matrix  produces 

13.55  * 100.0  / 150.0  = 9.03 

uhich  agiMJOS  with  Cook’s  results  for  the  finite  element  model. 

CASE  2: 

(a)  Cook’s  finite  element  model;  9.3 

(b)  Cook’s  norm,alized  beam  theory  value;  102. G 

(c)  Result  from  generated  stiffness  matrix:  92.8 

(d)  Un-norma I i zed  beam  theory  result; 

(deflection  due  to  bending  »)  1000 
+ (deflection  due  to  shear  =)  30 

- 1030 

Scaling  the  results  from  the  generated  stiffness  matrix  produces 

92.8  Vf  102.G  / 1030.0  = 9.2 

which  agrees  well  with  Cook’s  results  for  the  finite  element  model. 
The  deviation  here  is  most  likely  due  to  the  difference  in  the  values 
calculated  using  beam  theory. 

Thus,  the  values  calculated  using  the  generated  template  agree 
quite  well  with  the  tabulated  results. 


3.4.3  Isoparametric  Formulation 


* SYSTEM  INITIAI-IZATION  * 

WELCOME  TO  VnWnV*  VERSION  1.0 

It  is  now  WEDNESDAY  DECEMBER  7, 1977  10:0:7 

The  current  file  is  (SYINIT,  FCNI 

The  current  device  and  username  is  (DSK,  AKIGI 


Report  problems  to 
ALAN  R.  KORNCOFF 
DEPT.  OF  CIVIL  ENGINEERING 
CAITNEG I E-MELLON  UNIVERSITY 
CMU-lOA,  AKIG 


I LLUSTHAT  I Vfi  EXAI1I-LES 


1 (5r III i Melt e ell  I input  with  a SEMICOLON  - *■,’ 
Input  for  HELP 

I npu  t your  LOG I N NAME 
? 

AKIG; 

Is  AKIG  correct  ? (YES;  or  NO;) 

? 

YES; 

GREETINGS  AKIG 
* METHOD  SELECTION  * 


The  available  formulations  include 

(1)  THE  ISOPARAMETRIC  METHOD 

(2)  THE  GENERALIZED  COORDINATE  METHOD 


Please  enter  tlie  number  of  the  method  chosen  (1  OR  2). 
? 

1; 

I SOPARAME  TRIG  FORMULA  T I ON 
,V  ISOPARAMETRIC  FORMULATION  EXECUTIVE  >v 
,v  ISOPARARMETRIC  FORMULATION  INITIALIZATION  * 


Do  you  u i sh 
? 

YES; 

to  set  DUMP, 

BREAK  O'.  TRACE 

POINTS 

1.  ISOEIN 

2.  INPISO 

3.  SFNISO 

4.  BMOATA 

5.  MATERL 

B.  AUXTER 

7.  LIMTIS 

8.  JACOBN 

9.  BI1XISO 

10.  INTISO 

11.  DISIPP 

12.  OISPLY 

SET  DUMP  POINTS 

Designate  selected  PHASES  by  entering 

the  associated  integer  INDEX  or  ’ALL’  for  all  phases. 


Type  ’END’ 
? 

END; 

to  TERMINATE. 

1.  ISOEIN 

2.  INPISO 

3.  SFNISO 

4.  BMOATA 

5.  MATERL 

G.  AUXTER 

7.  LIMTIS 

8.  JACOBN 

9.  DUX ISO 

10.  INTISO 

11.  DISIPP 

12.  OISPLY 

SET  BREAKPOINTS 

Designate  selected  PHASES  by  entering 

the  associated  integer  lIClEX  or  ’ALL’  for  all  phases. 

Type  ’END’  to  TERMINATE. 


1 l.LU'oTRA  t I VH  L'XAIII'LtS 


1.  l&aillN  2.  INPISO  3.  SFNlliO 

5.  MATPRl.  G.  AUKTFR  7.  L HIT  IS 

3.  BHKISO  10.  INTISO  11.  DISIPP 
SET  .TRACE  POINTS 

Oesi finale  selected  PHASES  by  entering 
the  associated  integer  INDEX  or  ’ALL’  for 
Type  ’END’  to  TERMINATE. 

? 

END; 


A.  BI1DATA 
8.  JACOBN 
12.  DISPLY 


all  phases. 


* PROBLEM  PARAMETER  SPECIFICATION  - ISOPARAMETRIC  * 

Input  the  NUMBER  OF  ELEMENT  NOOES 
? 

4; 

Input  the  NUMBER  OF  DEGREES  OF  FREEDOM  PER  NODE 
? 


Input  the  NUMBER  OF  NATURAL  COORDINATES 
? 

2; 


Input  the  vector  of  the  NAMES  OF  THE  NATURAL  COORDINATES 
There  should  be  2 elements 
ELEMENT  1 - 
S; 

ELEMENT  2 - 
T; 

Input  the  vector  of  the  NAMES  OF  THE  GLOBAL  COORDINATES 
There  should  be  2 elements 
ELEMENT  1 » 

X; 

ELEMENT  2 - 
Y; 

Input  the  vector  of  the  NAMES  OF  THE  DISPLACEMENT  VARIABLES 
There  should  be  2 elements 
ELEMENT  1 - 
U; 

ELEMENT  2 - 
V; 


* SHAPE  FUICTION  PROCESSOR  - ISOPARAMETRIC  Vr 

ENTER  the  terms  of  the  SHAPE  FUNCTION  ordered  from 
node  1 through  node  4. 

The  4 elements  ui I I be  prompted  for 
ELEMENT  1 - 


ll.l.UlUnATIVt:  flKAIII-LtS 


1/A*(1hS)*(1-T) ; 

LLtiHliNT  2 = 
l/A*f.l+S).vH+7); 

CLtHtNT  3 « 

1/4*{1-S)*(1+T); 

ELEnCNT  4 = 

1/4*(1-S)*(1-T); 

yUAPE  KUNCTION  nODlFJCATJON 
The  OPTIONS  are 

(1)  DISPLAY  THE  SHAPE  FUNCTIONS 
{?)  tlOniFY  THE  SHAPE  FUICTION5  ' 

(3)  TERMINATE  THIS  FUNCTION 
Enter  the  NUMBER  ASSOCIATEO  with  the  CHOSEN  OPTION 
? 

1; 

The  terms  of  the  SHAPE  FUNCTION  ere 
Term  1 

(S  + 1)*(1  - T) 


4 


Term  2 


(S  + 1)>V(T  + 1) 
4 


Term  3 


(1  - S)»V(T  + 1) 


4 


Term  4 


(1  - S)*(l  - T) 


4 

Enter  the  NUMBER  ASSOCIATED  with  the  CHOSEN  OPTION 

7 

3; 

* B MATRIX  DATABASE  GENERATION  * 

The  OPTIONS  for  specifying  STRAIN  COMf’ONENTS  <tre 

(1)  USER-SUI>PLIED  VALUES 

(2)  LIBRARY  VALUES 

Ple^ise  ENTER  the  NUMBER  ASSOCIATEO  WITH  YOUR  SELECTION 

7 

2: 


The  LIBRARY  OPTIONS  for  SPECIFYING  STRAIN  COMPONENTS  arc 
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I1.LU15T»M7IV(-:  exauflls 

(1)  USER-SUPPLIED  VALUES 

(2)  ONE  DlflENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 

(4)  /PLANE  STRAIN 

(5)  AXISYUIIETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 30 
Please  ENTER  the  NUflBER  ASSOCIATED  UITH  YOUR  CHOICE 


* HATERIAL  PROPERTIES  SELECTION  * 

The  options  for  the  selection  of  the  MATERIAL  PROPERTIES  MATRIX  arc 

(1)  USER-SUPPLIED  MATRIX 

(2)  LIBRARY  MATRIX 

Please  enter  the  NUMBER  ASSOCIATED  UITH  YOUR  SELECTION 
? 

; 

The  LIBRARY  OPTIONS  for  SPECIFYING  MATERIAL  PROPERTIES  are 

(1)  USER-SUPPLIED  VALUES 

(2)  ONE  DIMENSIONAL  ELASTICITY 

(3)  PLANE  STRESS 
(A)  PLANE  STRAIN 
(5)  AX  I SYMMETRIC 

(G)  LINEAR  ISOTROPIC  ELASTICITY  - 3D 
Please  ENTER  the  NUMBER  ASSOCIATED  UITH  YOUR  CHOICE 


Vf  AUXILIARY  TERM  PROCESSOR  * 


ENTER  ELEMENT  VOLUME  MODIFICATION  FACTORS  OR  AUXILIARY  TERMS 
TYPE  ’ENDS’  to  TERMINATE 
? 

TH; 

? 

END; 

AUXILIARY  TERM  MODIFICATION 
The  OPTIONS  are 

(1)  DISPLAY  THE  AUXILIARY  TERMS 

(2)  MODIFY  AN  AUXILIARY  TERM 

(3)  TERMINATE  THIS  FUNCTION 

Enter  the  NUMBER  ASSOCIATED  uith  the  CHOSEN  OPTION 


? 

3; 


1 


I I.I.U5TRAT  I VH  tZXAtll-LL'S 
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i 


* 


f-NfER  thn  LU1ITS  OF  lUfEGRATION  for  NATURAL  COORDINATE.  S 
LONER  LIMIT  - 
? 

-i;  . / 

UPPER  LIMIT  = 

? 

i; 

ENTER  the  LIMITS  OF  INTEGRATION  for  NATURAL  COORDINATE.  T 
LONER  LIMIT  = 

? 

-1; 

UPPER  LIMIT  = 

? 

1; 

* JACOBIAN  GENERATION  - ISOPARAMETRIC  * 

* B MATRIX  GENERATION  - ISOPARAMETRIC  * 


BREAKPOINT  FOR  PHASE  B11XISO  ENCOUNTERED 
TYPE  ’EXIT;’  TO  RESUME 


35  msec. 

(MACSYMA-BREAK) 
^CLOSEFILE (AKIG. IQCMD) ; 


< To  avoid  exceeding  storage  space  provided  by  the  present  version  of 
MACSYMA,  processing  nas  suspended  after  the  determination  of  the  (Bl 
matrix.  The  record  file  uas  closed,  generated  VvtIucs  were  saved  on 
file  (AKIG,  IQVALS) , the  current  job  iios  terminated  and  a new  job  was 
initiated.  The  values  were  then  restored  and  selected  functions  from 
the  integration  processor  and  display  pre-processor  phases  were  then 
Sused  to  produce  the  integrand  of  equation  (2)  (section  1.2.1);  > 

(D2)  [DSK.  AKIGI 

{C3)  LOAOFILE (AKIG, IQVALS); 

AKIG  IQVALS  DSK  AKIG  being  loaded 
I oad i ng  done 

(D3)  DONE 


(C4)  LOADFILEdNTISO.FCN); 
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II.I.U^TRA1IVM  EXAIin.ES 

(C.5)  l.OAnriLE(DISIPP.FCN)! 

DI51PP  FCN  D5K  AKIG  bninrj  loodecl 
locicling  done 

(D5)  , / DOME 

(CG)  BI1ATRIX: FACTOR  (BI1ATRIX)# 

(C7)  FUNCTIONAL : QUAURAT J CF0RI1 IDI1ATR1 X. I1ATER1  Al.StlATRI X) i 


GENtlAT  FASL  DSK  I1AXOUT  beiny  loaded 
loading  done 

(C8)  FUNCTIONAL [1.1]; 

2 

(Y  T-Y  T-Y  S + Y S-Y  + Y) 

4 3 3 2 4 2 

(D8) 

G4 


2 

(1-NU)(X  T-X  T-X  S + X S-X  + X) 

4 3 3 2 4 2 
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< EKpression  (08)  represents  the  coefficient  before  it  is  multiplied 
by  any  auKiliary  terms  which  are  functions  of  the  natural 
coordinates.  > 

(C3)  I NTVARSI N5TACK  (NATURALCOOROS, STACK, STACLTOI NTER) ; 

(09)  DONE 

(CIO)  FUNCT I ONAL ; FUNCT I ONALMULT 1 NTVARS  (STACK , STACKPO I NTER, FUNCT I ONAL ) 8 


(Cll)  FUNCTIONAL [1,1]; 


(Y  T-Y  T-Y  S + Y S-Y  + Y) 
4 3 3 2 4 2 


(Dll)  - 8 (• 


G4 


(1  - NU)  (X  T-X  T-X  S + X S-X  + X) 

4 3 3 2 4 2 

+ ) 
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/(((X  -X)Y  +(Y  -Y)X  +(X  -X)Y  +(Y  -Y)X)T 
214  124  123  213 


+ ((X  - X)Y  +(Y  - Y)X  + X Y - Y X - X Y + Y X)S 

3 24  2'34  13  13  12  12 


1 


I 


II  l.USTIMTIVl-:  FXAHR.FS 


•t  (X  - X ) Y I-  (Y  - Y ) X - X Y + Y X - X Y + Y X ) 
1 3A  3 1 It  23  23  12  12 

(C12)  FACTOR;  nULTFACTOni  SO  (STACK.  STACKTOINTER) ; 

E TH 

(012)  - 

• (NU  - 1)  (NU  + 1) 

(Cl 3)  CLOSEFILEdQCMD.DISP); 


< As  thR  output  is  VO  1 utni nous,  only  the  first  coefficient  (expression 
(Dll)  ) and  the  common  multiple  (expression  (012)  ) are  displayed.  > 

< It  is  to  be  noted  that  although  integration  could  not  be  achieved, 
because  of  the  restrictions  placed  upon  storage  space  in  this 
exper i fnnnt.-j I version  of  IIACSYnA,  the  results  produced  are  in  the 
format  appropriate  for  the  hybrid  syinbo  I i c-numer  i c inetliod  dicussed  in 
section  2.1  . That  is,  a "matrix  template"  could  be  produced  for 
expressions  such  as  (Dll);  any  numerical  quadrature  program  could  call 
ouch  a template,  supplying  numerical  values  for  the  element 
coordinates  as  uel  I as  for  the  natural  coordinates  s and  t at  the 
selected  quadrature  points.  > 


1LLU5TRAT|V(-:  CXAni'LES 


3.5  Algpbrair.  Format  of  Gencrotpcl  Exprossions 

IIACSYUA  performs  no  automatic  simplification  of  nyrnliolic 
expressions.  Insteacl,  the  user  is  provided  with  . a class  of  functions, 
with  which  he  can  interactively  produce  output  formats  to  his  liking. 
Computation  time  is  generally  high  for  these  operations. 

It  would  be  a major  research  effort,  in  itself,  to  invest  VoWnVsv 
with  a set  of  heuristics  that  would  enable  it  to  make  decisions 
concerning  "optimal"  output  formats.  Current  results  are  produced  by 
applying  the  FACTOR  function  to  the  output  of  the  integration 
processors. 

Some  of  the  various  output  formats  which  may  be  generated  are 
demonstrated  in  the  run  below.  Several  equivalent  forms  are  presented 
in  algebraic  and  FORTRAN  card  image  formal.  The  example  employs 
element  11,1]  of  the  factored  stiffness  mati'ix  for  the  CST, 
generalized  coordinate  formulation. 

(D2)  [OSK,  AKIGI 

(C3)  L0A0FILE(y.3I,VALS); 

%3I  VALS  DSK  AKIG  being  loaded 
loading  done 

(03)  DONE 

< Expresions  (DA)  and  (D5)  are  those  currently  produced  by  *****.  > 

(C4)  FTLll: FUNCTIONAL [1,11; 

2 2 2 2 

(04)  -(X  NU-2X  X NU  + X NU-2Y  + 4Y  Y - X +2X  X 

3 23  2 3 233  23 

2 2 

- 2 Y - X )/4 
2 2 

(C5)  FORTRAN (%); 

FORTRA  FASL  DSK  flACSYN  being  loaded 
loading  done 

- (X  (3)  5V*2*f'HJ-2*X  (2)  >vX  (3)  >’fNU+X  (2) **2-.vNU--2*Y  (3)  **2+4*Y  (2)  *Y  (3)  -X  (3) 
1 *2+2*X  (2)  >vX  (3)  -2>vY  (2)  >v*2-X  (2)  »-nv2)  /4 

(05)  DONE 

< The  application  of  the  EXPAND  function  will  cause  products  of  sums 
and  exponentiated  sums  to  be  multiplied  out,  numerators  of  rational 
ex(5ressions  which  arc  sums  to  be  split  into  their  respective  terms, 
and  multiplication  to  be  distributed  over  addition  at  all  levels  of 
its  argument.  (101  > 
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(IJ6)  EXPAND  (Fill  1); 

2 

X NU  XX  NU 
/ 3 2 3 

(DG) + 

4 2 


9 7 7 

X NU  y"  X 

2 3 3 

+ Y Y + — 

4 2 2 3 4 


2 2 

X X Y X 

2 3 2 2 

2 2 4 


(C7)  FORTRAN (%); 

-X  (3)  ,v*2vrNU/4-. X (2)  :3)  .vNU/2-X  (2)  >v*2*NU/4+Y  (3)  >v,v?/2- Y (2)  ,vY  (3)  +X  (3) 

1 **2/4-X (2) *X (3) /2+Y (2) **2/2+X (2) **2/4 

(D7)  DONE 


< RATSIMP  "ra  t i oPci  I I y"  aiinplifies  its  nryumcnt  and  all  of  its 
subexpressions  including  arguments  to  non-rational  functions.  The 
t'esult  is  returned  as  the  quotient  of  tuo  polynomials  in  a recursive 
form,  i.e.  the  coefficients  of  the  main  variable  arc  polynomials  in 
tite  other  variables.  [11]  > 


(CS)  RATSniP(FTLll) ; 
CDS) 


2 2 2 2 2 2 
(X  - 2X  X + X)NU-2Y  + 4Y  Y - X + 2X  X ~2Y  - X 

3 232  3 233  23  22 


4 

(CS)  FORTRAN  (Tf); 

- ( (X  (3)  **2-2*X  (2)  *X  (3)  +X  (2)  **2)  *NIJ-2*Y  (3)  **2-»4*Y  <2)  *Y  (3)  -X  (3)  **2+2 
1 *X  (2)  *X  (3)  -2*Y  (2)  »-ov2-X  (2)  **2)  /4 
(09)  DONE 

< Horner  (exp)  ui  I I convert  exp  into  a rearranged  representation  as  in 

Horner’s  rule.  (12)  > 

(Cl (3)  HORNER (FTLll); 

OPT  1 11  FASL  DSK  I1AX0UT  being  loaded 
I oad i ng  done 

2 

(D10)  (((2X  - X)X  - X)NU  + Y (2Y  -4Y)+X  (X  - 2X) 
2332  33  233  2 

2 2 

+ 2 Y.  + X )/4 
2 2 

(Cll)  FORTRAN  (’,:); 

( ( (2*X  (2)  -X  (3) ) *X  (3)  -X  (2)  **2)  vrNU+Y  (3)  * (2*Y  (3)  -4*Y  (2) ) +X  (3)  * (X  (3)  -2 
1 *X  ( 2 ) ) +2*Y  ( 2 ) **2+X  ( 2 ) .v»v2 ) /4 

(Dll)  DONE 

< It  is  to  be  noted  that  the  range  of  the  number  of  multiplications 
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1 Chapter  4 


SUIIflARY  AND  CONCLUSIONS 


4.1  Summary 

A MACSYUA' based  processor  for  aiding  in  the  systhesis  of  stiffness 
matrices  for  finite  element  applications  has  been  introduced.  It  is 
hoped  that  the  system  ui  I I significantly  reduce  the  computational 
effort  in  the  generation  of  elements,  as  well  as  permit  investigations 
which  would  be  intractable  by  current  methods. 

It  is  also  hoped  that  this  project  has  suggested  direction  for 
future  research  and  demonstrated  the  potential  of  applying  flACSYflA  to 
other  engineering  problems. 

4.2  Suggestions  for  Future  work 

Some  possible  areas  of  future  research  include  the  implementation 
of  the  extensions  listed  below.  The  ramifications  of  many  of  these 
extensions  are  discussed  in  section  4.3. 

(a)  Permit  the  formulation  of  consistent  mass  matrices,  load  vectors, 
geometric  stiffness  terms  amd  thermal  effects  (see  section  2.9). 

(b)  Admit  second  and  higher  derivatives  in  the  specification  of  strain 

components.  •-  ,■ 

(c)  Remove  the  restriction  that  each  degree  of  freedom  must  be 
represented  by  the  same  shape  function. 

(d)  Remove  the  restriction  that  each  node  must  have  the  same  number  of 
degrees  of  freedom. 

(e)  Investigate  more  efficient  methods  of  performing  the  symbolic 
operations. 

(f)  Investigate  criteria  and  methods  to  produce  better  expression 
optimization  in  the  final  matrix  template. 

(g)  Investigate  extensions  to  incorporate  material  and  geometric 
non- I ineari ties. 


(h)  Investigate  alternative  formulations. 

(i)  Permit  the  specification  and  production  of  additional  output 
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formats  directly  compatible  uith  the  calling  sequences  of  a variety  of 
Invoking  analysis  programs. 

, /' 

4.3  Conclusions 


It  is  felt  that  the  current  version  of  vnViVVov  has  fulfilled  the 
problem  statement  and  design  objectives  within  the  limitations  and 
constraints  of  the  operating  environment.  This  is  not,  however  a 
statement  that  the  problem  area  has  been  exhausted.  Further  research 
in  the  areas  presented  in  the  section,  FUTURE  WORK  are  definitely 
necessary.  , - . 


Those  areas  offering  the  greatest  return  in  extending  the 
capabilities  of  the  system  for  the  least  investment  in  development 
costs  are  to:  (a)  admit  second  derivatives  in  the  specification  of  the 
stain  components;  (b)  remove  the  restriction  that  each  degree  of 
freedom  must  be  reperesented  by  the  same  shape  function  and  (c)  permit 
the  formulation  of  consistent  mass  matrices  and  load  vectors. 
Investigations  of  the  admissibility  of  material  noni  i near  1 1 i es  may 
also  prove  fruitful. 


The  implementation  of  many  algorithms  merits  re-examination. 
Operations  using  FOR  loops  should,  whenever  feasible,  be  replaced  with 
vector  operations.  An  examination  of  the  templates  generated 
indicates  that  attention  should  be  directed  towards  facilitating  the 
computation  through  the  identification  of  patterns  resulting  in  the 
repetition  of  stiffness  coefficients.  This  technique  has  been 
employed  in  Luftlll,  Anderson  and  Noor 14)  and  Anderson  and  Bowen  151 . 


Futhermore,  considerable  additional  improvements  in  execution 
speed  may  be  obtained  if  futher  research  is  carried  out  on  the 
optimization  of  the  matrix  template.  This  area  represents  a 
significant  research  effort  in  itself.  Refinement  of  these  output 
formats  could  readily  permit  execution  time  checks  on  numerical 
results.  One  such  check  could  determine  if  the  determinant  of  the 
Jacobian  is  equal  to  zero.  The  symbolic  form  of  the  determinant  is 
generated  during  processor  execution  and  would  be  Imbedded  in  an 
appropriate  boolean  expression. 


An  Important  ramification  of  this  research  is  the  demonstration  of 
the  potential  of  MACSYflA  in  the  solution  of  engineering  problems. 
Though  proficiency  in  its  use  requires  an  appreciable  overhead  to 
engineers  acquainted  only  with  FORTRAN,  the  rewards  are  enormous.  In 

addition  to  the  ease  and  variety  of  computation  possible,  the  use  of 
flACSYMA  permits  investigation  of  general  parametric  representations 
heretofore  unattainable.  Present  limitaions  in  Its  operation  will 
diminish  as  additional  storage  space  is  made  available.  There  are 
indications  ^ that,  through  the  institution  of 
processor  on  the  computer  supporting  flACSYtIA,  a 
provide  a greater  alotment  of  storage  space. 


a virtual  memory 
future  version  will 


SLJIIMARY  AND  CONCLUSIONS 


In  this  study,  It  uas  found  that  Integration  uas  proportionally 
the  moot  costly  component  of  the  symbolic  process.  The  hybrid 
symbolic-numerical  scheme  presented  could  eliminate  this  problem,  at 
the  enpense  of  a larger  amount  of  numerical  computation  at  enecution 
time.  The  tradeoffs  between  the  symbolic  and  hybrid  schemes  merit 
further  study. 

Assessing  the  present-  lack  of  tabulated  closed-form  formulations 
of  finite  element  stiffness  matrices,  it  is  apparent  that  *>v*>’nv  can 
make  a significant  contribution  to  the  present  field  of  knowledge, 
permitting  display  and  examination  of  the  closed  form  expressions. 

% 

The  impact  upon  contemporary  analysis  programs  which  employ 
numerical  quadrature  merits  investigation.  A sense  of  the  potential 
is  apparent  when  one  considers  that  these  traditional  computation 
techniques  require  an  investment  of  effort  proportional  to  the  number 
of  elements  used  in  each  problem  to  be  modeled.  Alternatively, 
need  be  used  only  once  for  each  element  type.  The  stiffness  matrix  so 
generated  is  valid  for  any  problem  requiring  that  element  type. 


p 
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FOOTNOTES 


tl3  Hathlab  Grou|>,  "IIACSYMA  Primer''  (Combrid^jc:  Mathliab  Group, 
l aboratory  for  Computer  Science,  llaar.achuaetto  Inotitutc  of 
Technology,  1975),  p.  1, 

12)  Richard  Bogen,  "riACSYflA  Reference  Manual"  (Cambridge:  Math  lab 
Group,  Laboratory  for  Computer  Science,  Massachuset to  Institute  of 
Technonogy,'  1975) , p.  2. 
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Ibid., 

P* 

2. 

14) 

I bi d. , 

P- 

2. 

15) 

Ibid. , 

p. 

3. 

16) 

Ibid. , 

P- 

125. 

[7] 

Ibid., 

p. 

40. 

[8) 

Ibid. , 

p. 

42. 

19) 

Ibid., 

p. 

81. 

[10] 

Ibid. . 

. P- 

40. 
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. P- 

42. 
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Ibid. , 
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OPERATION  DETAILS 


I.l  INPUT  PROTOCOL 

Important  bnd  unique  features  of  data  entry  are  summarized  beloui: 

(a)  All  input  is  format  free.  Input  lines  are  terminated  by  a 

semicolon  and  not  by  a return.  Input  may  be  in  upper  or  lower  case. 

(b)  Uhcn  typing  command  lines,  depressing  the  "rubout"  or  "delete"  key 
deletes  the  previous  character.  By  typing  "ctrl  K",  the  user  obtains 
a copy  of  the  current  command  line  free  of  any  echoed  erasures.  The 
two  characters  ??  delete  the  whole  command  line  and  causes  the  line 
number  to  be  redlnplnyed  [A], 

(c)  In  general  input  is  prompted  for  with  a question  mark. 


1.2  ACCESSING  THE  SVSTEfl 

The  procedure  for  gaining  access  to  vn’nvvnv  involves  three 
processes: 

(1)  Establishing  a connection,  via  the  ARPANET,  with  the  fIC  computer 
at  H. I . T. { 

(2)  Accessing  VnWnv*; 

(3)  Terminating  the  connection  to  tIC. 

1.2.1  Establishing  the  ARPANET  Connection 

This  procedure  is  installation  dependent  and  Is  outlined  os  It 
would  be  enacted  on  the  PDP-IBA  computer  in  the  Computer  Science 
Department  at  Carnegie-flel  Ion  University.  Specific  details  of  this 
particular  procedure  can  be  found  in  reference  1121. 

(1)  Login  to  the  computer  at  CflU,  henceforth  refered  to  as  the  host. 

(2)  Establish  the  ARPANET  connection  with  MC  I.T.  by  typing 

imp  telnet  mit-mc 


1.2.2  Accessing  vrvwhviv 

Details  of  logging  onto  HC  and  of  ITS,  the  operating  system  on  HC, 
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can  be  found  in  reference  1111.  The  uoer  should  read  sections  ll.B, 
and  II.G  as  an  absolute  mininmm.  The  procedure  for  initiating 
execution  of  >vvn‘nvv«  is  outlined  as: 

. / 

(1)  Login  to  I1C 

(2)  Create  a job  and  load  a copy  of  flACSYHA  by  typing 

! A 

(3)  MACSYIIA  uill  prompt  uith  line  Cl.  Respond  uith 

batch!  [fstaricmrjtdskiaklg]  .on)  i 

This  command  does  not  require  a carriage  return  since  it  is  a MACSYflA 
command.  Execution  of  is  initiated. 

(4)  The  logout  procedure  uill  depend  upon  the  option  chosen,  in  phase 
TERMIN,  to  terminate  the  system: 

If  option  3 (TERMINATE  THE  RUN.  THE  JOB  AND  LOGOUT)  is  chosen,  proceed 
to  step  (5). 

If  option  2 (TERMINATE  THE  RUN  AND  TIC  JOB)  is  selected,  control  will 
return  to  ITS  command  level  (see  Lewis  til)).  Type 

: I ogout 

and  proceed  to  step  (5). 

If  option  1 (TERMINATE  THE  RUN)  is  requested,  control  uill  return  to 
MACSYMA  command  level  and  the  user  uill  be  prompted  for  a MACSYMA 
command  with  a line  of  the  form  “ (Ci)  ".  To  exit,  type  a "ctrl  Z" 
and  follow  the  procedure  for  option  2. 

(5)  Return  to  the  host  by  entering  a "ctrl  taackarrow".  (On  some 
terminals  this  uill  be  a "ctrl  underscore"  or  "ctrl  shift  0". 

1.2.3  Terminating  the  Connection  to  MC 

This  process  Is  also  installation  dependent.  For  the  host  at  CMU,  the 
user  need  only  enter  ' 

imp  close/sel f 

and  then  logout. 


I .3  Addi t ional  Notes 


(a)  Execution  of  the  system  may  be  aborted  at  any  time  by  entering  a 
"ctrl  2".  Control  will  return  to  ITS  command  level  and  the  user  may 
restart  by  typing 
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ski  1 1 

and  then  foMou  the  (jrocedure  starting  with  stop  (2)  of  section  1.2.2 
or  exit  by  executing  the  procedure  starting  with  step  (5)  of  that  same 
sect  I on. 

(bl  The  Record  File  (section  2.B.1)  may  be  retrieved  for  printing  at 
the  host  by  using  the  file  transfer  program  (FTP)  which  is  documented 
i n reference  (121  . 

(1)  Uhlle  on  HC,  use  the  ITS  FIND  command  to  determine  the  file 
name  by  typing 

sfinci  users;  loginname  * 

where  "loginname"  is  replaced  by  your  flACSYMA  login  name  (see 
reference  (111).  This  will  display  the  directory  information  for  all 
record  files  which  the  user  has  created  by  executing  the  system.  From 
the  time  and  date  information  displayed,  the  user  can  determine  the 
desired  file.  The  file  name  will  be  your  login  name  followed  by  an 
integer.  The  file  with  the  greatest  integer  is  the  most  recently 
created. 

(2)  On  the  host  enter 

r f tp 

host  mit-mc 
user  users 

retr  users; loginname  integer 
qui  t 
qui  t 

where  "Integer"  is  the  value  displayed  by  the  FIND  command.  This  will 
create  a file  on  the  user’s  directory  on  the  host,  called  "USERS" 
which  will  contain  the  contents  of  the  record  file.  The  file  may  be 
dumped  to  the  line  printer  or  examined  using  one  of  the  editors.  The 
file  on  fIC  is  left  untouched  and  may  be  erased  using  the  ITS  DELETE 
command. 


I 


Unclassified 

'CCUfilTV  CL*»»IFtCATIOH  or  THIS  >l>8t  '»}>•»  Pit  Cnlmnd) 

f REPORT  DOCUMENTATION  PAGE 


) NUMfCA 

R-78-103^ 

4.  title  Cantf  SublltltJ 


U.  SOVT  ACCESSION  NO. 


Symbolic  Generation  of 

Finite  Element  Stiffness  Matrices 


17.  AUTMONi’#; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

S.  REClNIENT'S  CATALOG  NUMIEH 


S.  TYPE  OF  AEPOAT  A PEMOO  COVBMED 

Technical  Report 

I 

S.  PEPFOPMINO  OPS.  PEPOPT  NUMBEP 

R-78-103 

a.  contaaCT  OP  saant  uumsenctj 


Alan  R.  Korncoff 
Steven  J.  Fenves 

s.  pepfopmins  opsanization  name  and  address 

Carnegi e-Mel  Ion  University 
Pittsburgh,  Pennsylvania  15213 

tl.  CONTROLL'NG  OFFICE  NAME  AND  AOCRESS 

Office  of  Naval  Research 
Structural  Mechanics  Program 
Arlington,  VA  22217 

14.  monitoring  AG£NCY  name  ft  ADORESSf^/  di/xer*n;  /root  CcntrotUng  Offle^J 

Office  of  Naval  Research 
CMU  Campus  Representative 
5000  Forbes  Avenue 

Pittsburgh,  PA  15217 

l«.  OISTPISUTION  STATEMENT  fot  Hill  Haport) 

Distribution  unlimited  for  any  purpose  of  the 
United  State  Government 


N0014-76-C-0354^ 


PROORAM  element.  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 


NR-064-536 

12.  REPORT  DATE 

January,  1978 

13  NU^TbER  OP  PACES 

109 

IS.  security  CLASS^  re/  thTi  r«p«it; 


Unclassified 


tSa.  DECLASSIFICATION  DOWNORADING 
SCHEDULE 


17.  Distribution  statement  fot  th«  atatraci  mtmd  m Block  JO,  II  dl//»r«o(  trmm  Kmpmt) 


1 18.  Supplementary  notes 


It.  KEY  WORDS  rCenllmia  on  roaarao  alda  II  noeaaaafr  and  IdsnXiy  Ar  A»ooA  nuMAorJ 

Finite  Element  Analysis 

Symbolic  Computing  ^ 


■fr/ 


aNT^CT  fCanMmia  an  forarao  alda  II  naoaaadty  and  Mmullr  Ar  Al— A w— AooJ 

A symbolic  processor,  *****  (pronounced  HMvb  star^T^to  assist  in  the 
generation  of  stiffness  matrices  for  finite  elements,  based  on  a recently 
developed  symbolic  processor.  Is  presented.  Operations  are  performed 
upon  el eme^ characteristics  and  material  properties  in  symbolic  form  to 
produce  a ""matrix  tempi  ate  ,'***cbns  is  ting  of  the  algebraic  expressions  gener- 
ated for  the  stiffness  coefficients  as  functions  of  the  problem  parameters  4 


EDITION  OF  < NOV  At  IS  obsolete 
S/N  0102.LF4l4.SA01 


Unclassified 

security  CLASSIFICATION  OF  THIS  PAOE  rBItan  Dim  I 


Unclassified 


SeCURITV  classification  S*  this  AASe  rNTicn  Dmim  Enfr*<l) 


literal  form.  The  template  may  be  evaluated  for  a given  element  by 
binding  these  symbolic  forms  to  the  numerical  values  associated  with  a 
specific  element.  The  evaluation  process  is  further  facilitated  by  permitting 
specification  of  a variety  of  output  formats  for  the  resulting  matrix  template. 
Required  input  is  minimized  by  automatically  synthesizing  the  constituent 
matrices  of  the  formulation  from  user-supplied  specifications  of  shape 
functions,  material  properties  and  stress-strain  relationships,  all  in 
symbolic  notation. 


The  processor,  written  in  MACSYMA,  is  highly  interactive  providing 
prompts  for  user  input,  enumeration  of  available  program  options,  and  exten- 
sive on-line  assistance. f-The  user  may  input  a "?”  in  place  of  a prompted 
input  to  request  instructional  text.  The  file  handling  capabilities  of 
MACSYMA  are  utilized  to  retain  a complete  record  of  each  program  run.  These 
records  facilitate  the  handling  of  diagnostics,  assist  in  further  processing 
and  permit  the  generation  of  satisfies  valuable  for  system  development. 

Error  checking  is  accomplished  through  semantic  checks  built  into  the  program 
functions  and  syntactic  checks  pWformed  within  the  MACSYMA  operating 
environment. 

A partial  list  of  user  input  includes: 

1)  Method  Selection  - Isoparametric  or  generalized  coordinate 
formulations.  2)  Element  Parameters  - Number  of  nodes,  number  of  degrees 
of  freedom  per  node  and  related  terms.  3)  Material  Properties  - This  matrix 
may  be  selected  from  a library  of  standard  forms  (e.g.  plane  stress,  plane 
strain)  or  supplied  by  the  user.  4)  Strain  Specification  - Components  are 
entered  in  a user-oriented  calculus  notation  (e.g.  su/ax  is  input  as  D(u,x)). 
5)  Shape  Functions  - Shape  functions  may  contain  trigonometric  functions  and 
a large  class  of  intrinsic  functions  as  well  as  polynomial  terms.  6)  Output 
Control  Specification  - A description  of  the  output  format  of  the  generated 
matrix  template. 


Possible  output  forms  include  a tabular  display  of  the  matrix  coef- 
ficients in  symbolic  form  and  the  coefficients  in  FORTRAN  card  image  format. 

Background  material  includes:  The  objective  of  this  study;  The 
derivation  of  the  stiffness  matrices;  A summary  of  previous  research;  A 
brief  description  of  MACSYMA. 

Details  of  the  implementation  of  *****  cover:  The  design  objectives; 
Details  of  the  algorithms  used  and  how  they  were  implemented;  A description 
of  *****  and  its  limitations. 


Sample  runs  Include  the  formulation,  using  both  the  isoparametric  and 
generalized  coordinate  methods,  of  the  stiffness  matrices  for  a:  Bar  element 
with  constant  cross-sectional  area;  Bar  element  with  linearly  varying  cross- 
sectional  area;  Constant  Strain  Triangle  with  uniform  thickness;  Four  Node 
Quadrilateral . 

Conclusions  are  d-awn  and  recommendations  for  future  work  are  made. 
Appendix  I contains  notes  on  operating  and  accessing  the  processor. 
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